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Diaphragm-Type Micromanometer for Use on a 
Mass Spectrometer 


Vernon H. Dibeler and Fidel Cordero 


\ diaphragm-type micromanometer is described that is capable of measuring pressures 
in the range of 1 to 100 microns with a sensitivity of about 0.1 micron on the 50-micron scale. 
The displacement of the thin diaphragm is measured by the change in mutual inductance 
of two coils mounted above the center of the diaphragm. The capsule is easy to assemble, 
and the rugged construction reduces the effects of mechanical vibration. The value of the 
manometer is illustrated by a comparison of the sensitivities of the isomeric heptanes, 
octanes, and nonanes calculated from pressures measured in two ways: First by the con- 
ventional method of measurement with a mercury manometer of the pressure of a known 
volume of the vapors and expanding approximately 1,000 fold, and second, by measuring 
with the micromanometer the pressure of the expanded gases 


I. Introduction diaphragms was measured by direct coupling to a 
strain gage, the stiffness of which largely controlled 

The amount of a sample of gas admitted to a gas- | the sensitivity of the instrument. 
analysis mass spectrometer is usually measured with We have constructed a diaphragm-type manom- 
a mercury manometer on a volume of a few cubic | eter capable of measuring pressures in the range of 
centimeters. The sample is then expanded approxi- | 1 to 100g with a sensitivity of about O.lu on the 


mately 1,000-fold into a gas reservoir preparatory to | 50g scale. The high sensitivity is due to the 
entering the leak manifold. The method is quite in- | specially designed diaphragm and to the absence 
applicable if the vapor pressure of a compound is | of mechanical coupling between the diaphragm and 
less than a few millimeters, yet the mass spectrometer | the displacement-measuring device. The displace- 
can, in principal, be used with substances having | ment was measured by the change in mutual in- 
vapor pressures as low as 0.01 mm or less. Even | ductance of two fixed, concentrically wound coils 
with substances having vapor pressures of several | due to the approach or retreat of the conducting 
centimeters, large errors may be introduced by the | surface of the diaphragm. 

expansion method because of departures from the 


ideal gas laws. Taylor and Young [1]' have de- Il. The Micromanometer 
scribed a method of introducing liquid samples by 
means of a micropipet. However, it has been noted A number of capsule designs were tried. Sylphons 


{2} that for some compounds the pressure of the | and corrugated diaphragms of the nesting type with 
vaporized liquid is considerably less than expected sufficient pressure sensitivity were found to be too 
due to adsorption on the interior surfaces of the | Sensitive to ordinary building vibration to be useful 
inlet svstem. Furthermore. no suitable method has | [n the final design, a single diaphragm with a working 
been described for introducing and measuring the | diameter of 2.87 in. was made from rolled brass 
vapor pressure of compounds that are solids at room | 0-001 in. thick. The diaphragm had 26 concentric 
temperature. The above sources of error are sub- V-shaped corrugations 0.01 in. deep with 90° angles. 
stantially eliminated if a suitable micromanometer | Figure 1 is a cross section of the capsule assembly 
be attached to the reservoirs for pressure measure- showing relative positions of the diaphragm, probe 
ments on the expanded vapors. coils, and coil support, The og mbly was made 

The construction of sensitive manometers for the | in the following manner. A 0.25-in. hole was drilled 
micron pressure range has usually taken the form of | !™ the center of a 0.08-in.-thic k brass disk, B. One 
modifications of the ordinary U-tube manometer | end of a 3-in. length of 0.50-in. outside diameter 
(3, 4, 5] or MeLeod gage. Most of these require | brass tubing, C, was hard soldered in position to 
careful control of temperature and bulky, rigid match the hole: the other end was tape ‘red to fit the 
mountings on vibration-free supports. Matheson | Outer me ‘mber of a 12/30 sts andard taper glass joint. 
and Eden |6] have described a sensitive manometer The disk was undercut 0.01 in. and the diaphragm, 
by using a pair of nesting diaphragms as the pres- A, soft soldered at the periphery. The parts of the 
sure-sensitive element. Their manometer was capa- | Stainless steel yoke, D, were machined to 0.125-in. 
ble of detecting a pressure differential between the thickness from bar stock of the proper dimensions 


- * 


inside and outside of the capsule of 0.001 mm in The mating surfaces were carefully fitted and assem- 


ae . ——— 
the range | mm to 1 atm. The displacement of the M. L. Greenough and W. E. Williams, An clectronic circuit for measuring the 
— displacement of pressure-sensitive diaphragms, J. Research NBS 46, M5! 
Figures in brackets indicate the literature references at the end of this paper. RP 2168 
9138375-——50 1 l 
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Figure 1 Cross element 


showing the relative positions of the components. 


section of the pressure-sensitive 


bled with machine screws. The yoke was attached 
to the diaphragm support by set screws in the 
collar, E. The probe coils were wound on an insu- 
lating coil form, F, and mounted on the underside 
of the top cross bar. Electrical tie-points were pro- 
vided in the form of Stupakoff through-panel insula- 
tors,G. Figure 2 shows the asse smbled pressure cell. 
The entire unit was mounted inside two glass hemi- 
spheres (not shown) sealed together by a low-vapor 
pressure wax. An outer 12/30 standard taper was 
sealed to the lower he ‘misphe re at the pole to receive 
C (fig. 1). A small quantity of Apiezon grease was 
used to make the seal. The outer standard taper 





Assembled pressure-sensitive unit, 


Figure 2 


was connected to the inlet system of the mass spev- 
trometer. Connection was also made to the hem- 
sphere through an isolating stopcock to prevent 
damage to the diaphragm during venting and also 
to allow periodic checks on the zero point of th 
diaphragm. The cable connections for the elec- 
tronic micrometer were made through a four-con- 
ductor glass press in the lower he ‘misphe re 

The pressure sensitivity measured at 
of the diaphragm was about 2 10 
of mercury. Conservatively, 
have satisfactory elastic properties with a permis- 
sible central deflection of 3 percent of its diameter, 
or 0.086 in., far in excess of any deflection contem- 
plated in service. It was deemed unnecessary there- 
fore to use materials with superior elastic properties, 
such as phosphor bronze or beryllium copper. The 
capsule was easy to construct, and the 0.08-in. base 
provided mechanical strength for handling, as well 
as stiffness to reduce the effect of mechanical vi- 
bration. 

A detailed description of the displacement-measur- 
ing instrument utilizing the mutual inductance mi- 
crometer appears in the accompanying paper (see 
footnote 2). <A brief description of the operation 
can be obtained, however, by referring to figure | 
The primary winding of the probe coil is excited with 
radio-frequency current. The voltage induced in the 
secondary is a function of the proximity of the dia- 
phragm surface as a result of eddy currents generated 
in the metal. This voltage is amplified by suitable 
circuitry, rectified and applied to a direct-current 
microammeter. The gain control of the amplifier 
permits a wide range of scale adjustments for the out- 
put meter. It also provides a means of adjusting 
the calibration of the unit from the deflection of 
the diaphragm caused by a known pressure difference 
across the capsule. 


III. Operational Details 


The manometer was operated with the capsule 
connected to the inlet reservoir. The pressure in 
the sphere enclosing the capsule was maintained at 
approximately 10°° mm of mercury. Zero adjust- 
ment was made when the reservoir was evacuated 
on the same vacuum line. Calibration of the micro- 
manometer was obtained by using the mass spec- 
trometer as a pressure-measuring instrument in the 
following manner. With the reservoir open to the 
leak manifold, a pressure of approximately 50 mm 
of n-butane measured with a mercury manometer 
was metered into the gas pipet provided as an in- 
tegral part of the inlet system oi the Consolidated 
mass spectrometer. This sample was then expanded 
1,000 fold into the reservoir. The gain control of the 
micromanometer was adjusted for a reading on the 
output meter corresponding numerically to the exact 
pressure in microns calculated from the manometer 
reading and the known expansion ratio. The m/ 
13 peak of the n-butane mass spectrum was then 
scanned by the mass spectrometer in the conven- 
tional manner. A portion of the gas in the reservoir 
was pumped out, the micromanometer reading was 
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IGURE 3 Calibration curves of the micromanometer for 
full-scale sensitivities of 25 and 50 microns 
left-hand ordinate refers to the open circles The right-hand ordinate 
tot olid , 


oted, and the 43 peak again scanned. Further 
epetitions of this procedure resulted in data illus- 
rated by the lower curve in figure 3, in which the 
rdinates are computed from the initial pressure 
nd the peak heights recorded by the mass spec- 
rometer. The upper curve was obtained in a similar 
nanner when the gain was adjusted for full-scale 
leflection of the output meter with a pressure of 
5u in the reservoir. At full sensitivity a pressure 
f approximately 104 was sufficient to give full- 
cale reading on the micromanometer. The noise 
evel for the 50-4 scale was approximately 0.14 and 
vas principally the result of vibration from pumps 
und other equipment mounted in the vicinity of the 
‘abinet housing the inlet system. Temperature 
ariations of 3 to 4 deg in the ambient room tem- 
erature had no effect on the performance. Zero 
lrifts over 8-hour periods were less than Iyamp, 
ind day-to-day adjustments of the full-scale calibra- 
ion were of the same order. The capsule has been 
ound to be fully protected against a pressure differ- 
ntial of 1 atm applied to the outside of the cap- 
ile No test has been made for large excess pres- 
ire on the inside Pressures of a few tenths of a 
nillimeter have had no effect on the mechanical zero 
Except for the pressure measurements, the sensi- 
Vities of various compounds were obtained by 
sing conventional techniques on a Consolidated 
spectrometer. The approximate amount of 
mple required was introduced into the inlet reser- 
ur. The reservoir was opened to the leak, and 1 
nute later the exact pressure indicated by the 
cromanometer was noted and the sweep started. 


Liss 


IV. Experimental Results 


Table 1 summarizes the sensitivities (ion current 
of maximum peak per unit of pressure) of the octane 
isomers relative to n-butane obtained to two ways 
First, the pressure of a known volume of vapor was 
measured with a mercury manometer and the vapor 
expanded 1,000-fold into a reservoir Second, the 
pressure of the expanded vapors was measured with 
the micromanometer attached to the reservoir, the 
latter opened to the leak. Columns one and two, 
respectively, give the names and vapor pressures of 


the 18 isomeric octanes. Column three gives the 
relative sensitivities calculated from the pressures 
of the substance and n-butane measured by a mer 


cury manometer [7]. The fourth column gives the 
relative sensitivites calculated from the pressures 
measured with the micromanometer [8]. Column 
five gives the ratio, 2, of the relative sensivities ob- 
tained by the two methods. The greatest value for 
R (1.9) is found in n-octane, which also has the 
lowest vapor pressure. The values of ? apparently 
decrease with increasing vapor pressure in a rather 
uniform manner for molecules with few side chains 
Similar data [8] were obtained for all of the heptane 
isomers and the nonane isomers. As expected, the 
effect was smallest for the heptanes. For n-heptane 
(vapor pressure, 46 mm) the value of R was 1.15. 
For the other heptanes (vapor pressures up to 100 
mm) the average value of 2 was 1.02. On the other 
hand, n-nonane (4.3 mm) gave a value for PR of 2.4, 
whereas 2,2,4,4-tetramethyl pentane (20.0 mm) gave 
a value for R of 1.5. That factors other than vapor 
pressure are operative for some molecules is indicated 
in table 1 by the low value of R (1.2) for 2,2,3,3- 
tetramethyl butane as compared with compounds of 
similar vapor pressure; for example, 2-methyl hep- 
tane (R=1.5) and 4-methyl heptane (2?=1.4) 
Another paper [2] reports some measurements made 
by introducing measured volumes of liquid cyclo- 
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hexane and liquid styrene into the reservoir of the 
mass spectrometer by means of a micropipet. In 
these experiments the pressure of the vapor measured 
by the micromanometer was always less than that 
computed from the amount of liquid introduced. 
For cyclohexane and styrene the descrepancy was 
about 13 percent and 30 percent, respectively, and 
was caused at least in part by adsorption on the walls 
of the rvservoir. 

It is often noticed that the pressure of some re- 
active gases and of some highly polar molecules 
gradually decreases on standing in the reservoir. 
In such cases satisfactory measurements are only 
obtained after waiting until a constant pressure is 
attained. A great advantage of the micromanom- 


eter is that continuous monitoring of the pressure 
in the reservoir is possible, and such changes are 
immediately visible and do not escape notice. 


The authors acknowledge the valuable discussions 
with Fred L. Mohler and W. G. Brombacher through- 


out various phases of the constructionof the micro) 
nometer. We are also grateful for the cooperatio 
the Engineering Electronics Section, in particular 
L. Greenough and W. E. Williams for the desig 
and construction of the electronic micrometer. 
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n Electronic Circuit for Measuring the Displacement of 
Pressure-Sensitive Diaphragms 


Maurice L. Greenough and William E. Williams 


An instrument has been developed for indicating the position of thin pressure-sensitive 


diaphragms 


As no appreciable load is placed on the diaphragm, the system is well adapted 


to the measurement of low absolute or differential pressures in the range from zero to 100 


microns The electronic 
5x 10-*ineh 


Instrument gives 


inductances 
may be 


I. Introduction 


An electrical distance-measuring system designed 
or application to sensitive diaphragms has been suc- 
essfully employed for pressure measurements in mass 
pectrometers.' The same system is also applicable 
to many other problems requiring the remote indica- 
tion of relative position. Specifically, the instrument 
lescribed is well suited for the measurement of dis- 
tances of 5X 10~* in.as a full-scale value Motions 
if less than 10~° in. may be detected. Displacements 
if 0.05 in. may be indicated by suitable designs in- 
orporating the same method of distance measure- 
ment 

The instrument, coupled to a pressure-sensitive 
diaphragm, is shown in figure 1. With the exception 
of the distance-measuring element itself, all of the 
electronic circuit is housed in the steel cabinet. The 
pressure cell and distance-measuring element may be 
seen inside the glass bell jar. In operation the bell 
jar is evacuated, and the measured pressure is that 
inside the cell. The cell is connected to the mass 
spectrometer by an inlet through the bottom of the 
glass bell. In other applications it may be preferable 
to interchange the unknown and = the 
vacuum. Instrument performance would not be 
affected 
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cromanomele system, showing inatru- 


ment, acuated glass bell, and pressure cell, 


full-scale 
Five tubes are incorporated in a circuit based upon the principle of the mutual- 
micrometer for the detection of diaphragm motion 
utilized in other applications requiring the indication of very small displacements 


indication for motions of less than 


The same instrumentation 


The distance-measuring element used for observing 
motion of the pressure-responsive diaphragm is of the 
mutual-inductance type” Two air-core are 
wound on a single form that is mounted in proximity 
to the diaphragm The presence of the metal dia- 
phragm adjacent to the coils produces partial shield- 
ing between them, and motion of the diaphragm as a 
result of pressure changes causes a variation in the 
mutual inductance between the windings. Suitable 
electronic circuitry is then utilized to provide the 
desired pressure indication 

The actual diaphragm motion is approximately 
0.0005 in. (500 microinches) for full-scale indication 
on the most sensitive meter range. For this motion, 
the mutual inductance is changed 1.2. percent 
While at first glance it would appear that other and 
more sensitive distance-measuring systems could be 
emploved, the ruggedness and adaptability of the 
present system give it advantages for this applica- 
tion. It would not be difficult, for example, to obtain 
a 10-percent change in capacitance between two 
plates at reasonable spacing. However, this initially 
greater percentage change would be masked by the 
capacitance of the interconnecting cable, or would 
require that part of the circuitry be mounted in close 
proximity to the pressure cell. Rather critical dia- 
phragm alinement would also be required, as both 
initial clearance and tilt affeet the magnitude of 
capacitance Variation \agnetic core type distance- 
measuring elements have the disadvantage of react- 
ing on the diaphragm, as well as requiring magnetic 
material attached to the diaphragm. 

The mutual inductance micrometer has advantages 
in this application in that it possesses linear response 
as a function of displacement and does not appre 
ciably load, or alter the characteristics of, the dia- 
phragm. As the distance-measuring element ex- 
hibits low-impedance, cable length and resulting 
capacitance is not critical. The device requires a 
reasonably flat area on the diaphragm approximately 
% in. in diameter, of a nonmagnetic, low-resistance 
metal at least 0.0005 in. thick. These conditions 
were satisfied for this application by fabricating the 
entire diaphragm in a suitable shape of 0.001-in 
brass sheet. (Construction of the pressure cell is 
described in the paper referred to in footnote 1 
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Fiewurt 2 Cut-away sketch oO pressure ce mutua 


The relative locations of the cell and 
micrometer windings are indicated in figure 2. A 
stainless-steel voke provides the mounting to hold 
the micrometer coil form and the pressure cell in the 
proper relation. The flexible diaphragm of the cell 
is nearest to the face of the supporting form of the 
distance-measuring element, or probe. As the pres- 
being measured is that inside the cell, zero 
absolute pressure causes the diaphragm to assume its 
most remote position with respect. to the probe 
windings. Increased pressure then produces expan- 
sion of the cell and consequent reduction in spacing 
between diaphragm and windings 

Normally, the initial clearance for zero pressure is 
set at about 0.005 in. between the diaphragm and the 
end of the probe \ pressure of 2504 ts sufficient to 
The diaphragm moves approximately 
20 microinches per micron of mercury pressure 
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close the gap 
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The electrical characteristics of the mutual indu.+- 
ance probe are shown in figure 3. When the m 
surface is in contact with the end of the probe, 1} 
coupling between the two windings is at a minimun 
because eddy currents induced in the metallic surfac 
set up an opposing electromagnetic field, whos 
cancellation effect is then greatest. Coupling |» 
tween the windings increases with the spacing |» 
tween probe and metal surface. The relation has 
the form of an initial or minimum coupling to whic! 
is added a component proportional to separation @ 
By suitable design, the proportionality of this com 
ponent is made linear to within 1 percent 
range of separation from 0 to 0.02 in. As the motior 
of the diaphragm is restricted to a few thousandths 
of an inch, the wide range of linear response is mor 
than adequate. However, it allows variation in thy 
initial position of the diaphragm without altering tiv 
sensitivity of the instrument, as the change of mutua 
inductance for a given incremental motion is 
stant. Initial clearance is set at the value of 0.005 
in. to provide some tolerance for misalinement 3 

The cutaway view of the probe in figure 2 illustrates @ 
the approximate geometry of the probe windings 
The primary winding is usually placed in the outer 
slot, which is further away from the metal surfac 
whereas the secondary ts located in the closer on 
The secondary is wound with several times the num- 
ber of turns on the primary \s its resistance is no 
critical, fine wire is used on the secondary in order to 
accommodate the winding in as narrow a ; 
possible The narrow slot and thin supporting wal § 
are desirable in that mounting of the windings closer § 
to the diaphragm is possible. In the design of th 
probe used in this equipment, the minimum value of @ 
coupling as set by these factors is equivalent to abou 3 
0.035in. The desirability of alow value of minimum @ 
coupling is that less precision Is then demanded of the 
circuits provided for its cancellation. However, thy 
actual change in mutual inductance in microhenrys 
the slope of the relation in fig. 3), is substantiall) 
unaffected by changes in minimum coupling 

Details of the probe and its windings are g 
the appendix, part A 
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II. Circuit Details 


The circuit diagram is shown in figure 4, and com- 
ponent specifications are listed in the appends 
part B. The source of radio-frequency power is a 
6V6 oscillator arranged for two-terminal operatiol 
by means of a nonresonant feedback transforme 
L, and Ly. The frequency of oscillation is deter- 
mined by reasonance of the shunt-fed tank circuit 
formed by C, in combination with the total indu 
tance of the primary windings of the current trans 
former, Z,, and mutual-inductance micrometer, / 
The oscillator output is approximately 1 amp. at 
frequency of 2.5 Mes. 

The current transformer serves several purposes 
Its chief function is to provide a means of obtaining 
a voltage that is proportional to the current in th: 
mutual-inductance micrometer. It is also necessary) 
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for limiting the power input to the distance-measur- 
ing element to avoid inductive heating of the 
diaphragm. Only about 10 percent of the voltage 
across C, appears across the primary winding of the 
micrometer element. 

To obtain constancy of output indication for a 
given pressure, it is necessary to regulate the ampli- 
tude of the radio-frequency current in the micrometer 
primary. It is accomplished in this instrument by 
regulating the oscillator output voltage, depending 
upon the stability of the primary inductance of the 
current transformer, probe, and cable to secure 
constant current. Frequently shift, arising from 
changes in the value of the tuning capacitor for 
example, produces a corrective variation in current, 
so that the output of the mutual-inductance microm- 
eter remains unchanged 

The regular circuit commences with a 6H6 half- 
wave peak rectifier, V,,, the output of which is fed 
to the grid of a 6SL7 direct-current amplifier, V,. 
A fixed reference voltage from a gaseous regulator 
tube, V, is applied to the cathode of the 6SL7. The 
difference between these two is the grid bias on the 
tube. The resulting plate voltage on the tube is 
transmitted through a 6SN7 cathode follower, V.,, 
to control the screen potential on the 6V6 oscillator, 
completing the regulating loop. The amplitude of 
oscillation is thereby automatically held at about 
102-v. peak. The loop gain of the regulator circuit 
is approximately 70, so that the effect of varying 
tube characteristics, line voltage, ete. 1s greatly 
reduced. 

In order to obtain a stable source of voltage for 
cancellation of the coupling at normal clearance, a 
secondary winding is placed on the current trans- 
former. Through a potentiometer, /?;, an adjustable 
r-f voltage source is obtained. Voltage from this 
source bears the same proportional relation to 
oscillator voltage and frequency as that from the 
secondary of the mutual inductance micrometer. 
This potentiometer provides adjustment for equality 
of current transformer and micrometer outputs, 
under conditions of zero differential pressure on the 
diaphragm. Thereafter,’ only relative geometrical 
changes can disturb the balance between these two 
voltages. Changes in inductance of the intercon- 
necting cables for example, can cause only identical 
effects in the two secondaries but no unbalance. 
The result is an extremely stable zero indication. 

The absolute magnitude of these voltages is 
proportional to the exciting current and frequency. 
As mentioned, the factors of exciting current and 
frequency are fixed by virtue of the constant-voltage 
regulator circuit in conjunction with stable induct- 
ances. Cable inductance increases the impedance 
of the inductive branch of the tank circuit and lowers 
the frequency tending to reduce current flow in that 
branch. Change in lead inductance therefore can 
produce similar variations in indicated output. 
However, lead inductance itself is small compared to 
that of the current transformer, so that minor varia- 
tions are negligible. Had variations from this 
cause been found appreciable, a step-up type of 





current transformer feeding the 6H6 rectifier cou d 
have been employed. For this instrument, the 
additional complexity was not felt warranted, as 
scale-factor stability was found to be good with the 
regulating circuit employed. 

When the primary winding of the probe is excited 
with high-frequency current, the output of the 
secondary is a voltage of the same frequency, whose 
magnitude is proportional to the mutual inductance. 
The secondary voltage is fed to one contact of a 
single-pole-double-throw vibrating relay actuated at 
60 cycles. The other contact is connected to the 
cancellation voltage of the same frequency from 7? 
The resultant voltage at the arm of the vibrating 
relay becomes a square-wave, amplitude-modulated 
carrier, whose modulation is proportional to dia- 
phragm excursion. Zero modulation is made to 
correspond to zero differential pressure by adjust- 
ment of 2. The modulated carrier is rectified, re 
sulting in a 60-cycle square-wave signal. This is 
amplified and fed back through radio-frequency 
isolation circuits to the contacts on the vibrating 
relay. The synchronous rectification so obtained 
provides current for actuation of a direct-current 
microammeter. 

The relay used was a 60-cvele vibrating one with 
mercury-wetted contacts. A steady direct-current 
from V, through FR, is put through one of its 


windings to bring the armature to the central 
position. The addition of 60-cycle alternating- 


current in the other winding causes equal dwell time 
on each of the contacts, producing symmetrical 
square-wave modulation. In practice, the magni- 
tude of FR, is chosen to provide equal dwell time. 

The square-wave component of rectifier output: 
is fed through an attenuator, Ry», for scale factor 
adjustment during calibration. Two scale ranges 
are provided with a 5:1 ratio (2), and Ry»). 

The square-wave is next passed through an ampli- 
fier stage, whose gain is stabilized with inverse feed- 
back by omitting the cathode bypass capacitor on 
R,,. Following the amplifier is a direct-coupled 
cathode follower that feeds an output transformer 
of about 2.5:1 stepdown to the center-tapped 
secondary. The square-wave voltage at the second- 
ary is connected through the low-pass filter, 7), 
R;, Cy, and Cy, to the same relay contacts that 
handle the carrier frequency signals. The purpose 
of the filter is to prevent loading of the carrier signal 
sources by the output transformer. Blocking ca- 
pacitors C, and C, prevent loading of 7; by the radio- 
frequency circuits. As resistances are required by 
amplifier loading considerations, and should _pref- 
erably be precision components, the use of wire- 
wound elements as r-f chokes served both functions 
No interference was observed between the low 
and high frequency circuits with the component values 
shown. The only critical point is that the radio- 
frequency choke from relay arm to ground have low 
60-cycle impedance to avoid spurious feedback 

Vibrating relays in this service were sometimes found to display large switch 


ing transients. Some methods of transient reduction are discussed in the ay 
pendix, part C. The particular one used was satisfactory in this respect. 
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The reflex connection acts as its own synchronous 
rectifier, saving the use of a separate detector tube. 


Ill. Performance 


The over-all response of the system is shown in 
figure 5, where the relationship between output cur- 
rent and gas pressure is plotted. As the relation is a 
direct’ proportionality, calibration at only two 
points, for intercept and slope, is required. Inter- 
cept adjustment, at zero pressure, is controlled by 
the zero potentiometer, R,, and is set so that zero 
output current is indicated for a perfect vacuum or 
zero differential pressure. The slope adjustment for 
magnitude of indication for known gas pressures ts 
determined by the scale control, R, ). Typical 
resultant curves for output current are given in 
figure 5, reproduced from the paper referred to in 
footnote | 

The stability of the instrument in the vicinity of 
zero pressure indication is good. As a result of the 

reuilts employed for obtaining cancellation voltage, 
is little of an electrical nature that can cause 
ero shift of serious magnitude. Obviously, the po- 

ntiometer used for adjustment of cancellation volt- 
ive must be stable. Line voltage changes of +10 
produce error readings of about 0.5 uw, in the form 
| a slow drift believed due to heater voltage effects. 
he cause of the variation observed is provably the 
pendence of the effective resistance of the 6H6 
nal rectifier, Vy», upon heater voltage. As the 


there 


irees of cancellation and micrometer signal volt- 
s are not of equal impedance, there may be a 
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Ficure 6 Construction details of probe. 


differential change in rectification efficiency in the 
half-wave, peak circuit employed. For this specific 
application, the magnitude of drift encountered was 
considered negligible. Hum pickup in the amplifier 
can also give rise to zero shift with line voltage vari- 
ations, inasmuch as the pickup is usually propor- 
tional to line voltage. Proper wiring and shielding 
will reduce any effects from this cause 

The stabilitv of scale factor is also good, because 
of the high gain of the oscillator regulator circuit, 
and the use of inverse feedback and regulated plate 
potentials on the square-wave amplifier. Line volt- 
age variations of 10 v result in approximately 1-per- 
cent indicated pressure errors, 


IV. Appendix 


\. Drawings showing the 
details of the mutual-inductance micrometer are presented Ils 
figure 6.@The windings are based upon the use of feed cables 
under 6 to 10 ft in length 

B Values of circuit co npo ents in figure 4 


constructional and winding 


Part number 

a. as © ( 0.0 1-uf, 100-5 paper condense! 

Cc, 0.00 1-yuf, 400-v mica condenser. 

Go. © 0.1-uf, 400-v paper condenser 

C, 0.000 1-yf, 400-v mica condenser 

Ce, Ca, Ce, ( 0.00 1-uf, 400-\ paper condenser 

(- 0.002-uf, 400-\ paper condenser 

Cye, Crs, C ( 20-yf, 150-\ electroly tic con- 
denser. 

ei 1,000-ohm potentiometer 

R,, R 1,000-ohm, 2-percent precision, 
multiple-pie, wire-wound re- 
sistor. 

R,, R, 20,000-ohm, 's-w, 10-percent tol- 


erance resistor 

















Part number 


Vibrating relay 


200, 000-olm, lw, 
tolerance resistor. 

150,000-ohm, 4-w, 
tolerance resistor. 

100,000-ohm, Lo-w, 
tolerance resistor. 

10,000-ohm, l-w, 5-percent toler- 
ance wire-wound resistor. 

1.0-megohm, 's-w, 10-percent tol- 
erance resistor 

1.0-megohm, composition potenti- 
ometer 

600,000-ohm, 2-percent, tolerance 
wire-wound resistor. 

150,000-ohm, 2-percent tolerance 
wire-wound resistor. 


10-percent 
10-percent 


10-percent 


500,000-ohm, ‘sw, 10-percent 
tolerance resistor 
1,500-ohm, 'o-w, 10-percent toler- 


ance resistor 

12,000-ohm, l-w, 10-percent toler- 
ance resistor 

8,.000-ohm, 5-w, power resistor. 

Double-pole-double-throw toggle 
switch 

6V6 tube 

68 N7 tube 

68L7 tube 

6H6 tube 

5Y3 tube 

OC3/VR-105 tube. 

In this application a Western 
Llectrie, of the make-before- 
break type, type D168479, was 
used. The direct current and 
alternating current are applied 
to the 3,300- and 700-ohm 
windings respectively. Relays 
of other manufacturers with 
contact ratings of at least 1 ma 
could be used. 


Part number 


Bas incon Winding details given in the ap- 
pendix, part A and figure 6. 

L 12 turns No. 27 F wire, close- 
wound on 54-in. diameter Bake- 
lite tube. 

Ly 2 turns No. 38 F wire, wound 
over L3 with three layers of 
varnished cambric insulation 
between L3 and L4. 

L 68 turns No. 27 F close-wound on 
l-in. diameter Bakelite tube. 

Ly 14 turns No. 27 F, close-wound 
over +300-v end of L5, with 
three layers of varnished cam- 
bric insulation. 

L- R-f choke, 750 wh, resistance 
less than 25 ohms 

Ls Filter choke, 8 h, 40 ma. 


C. There may be high transient peaks coinciding wit! 
switchover of the relay. In many relays or choppers of this 
type, the arm may briefly touch all contacts at once, or it 
may have a completely open period. At the time when all! 
contacts are joined, the combined radio-frequency voltage 
may not lie in the range between the individual values be- 
cause of phase differences. Reversal of the polarity of either 
signal may be found to reduce the amplitude of the peaks 
If a relay is used that has a momentary open-circuit condi- 
tion, as with the break-before-make type, large transients 
occur because the carrier frequency amplitude then drops to 
zero. The effect may be eliminated by connecting a condenser 
of a few micromicrofarads between the relay arm and some 
source of r-f voltage, such as the oscillator grid. Due to its 
high reactance, the condenser can only control the potential 
of the relay arm when it is not contacting either of the low 
impedance voltage sources. Proper adjustment is indicated 
when the transient peak is reduced to zero by establishing the 
open arm potential as being equal to the cancellation voltage 
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Research Paper 2169 


Heat of Polymerization of Alpha-Methylstyrene from Heats 
of Combustion of Monomer and Four Polymer Fractions 


Donald E. Roberts and Ralph S. Jessup 


Measurements have been made with a bomb calorimeter of the heats of combustion 
Allc®) at 25° C to form gaseous carbon dioxide and liquid water, of liquid alpha-methyl- 


stvrene 


of polymerization ( 


have been calculated from the heats of combustion 


l-methyl-l-phenylethylene) and four fractions of its sold polymer having number 
average molecular weights of (1) 5,400, (II) 3,750, (111) 2,200, and (1V) 1,300. 
of heat of combustion are: monomer, 1204.87 keal/mole; fraction I, 
1195.83; fraction III, 1195.50; and fraction IV, 1194.74 keal/CyHyo unit. 
AlTp®) at 25° C of liquid alpha-methylstyrene to its solid polymers 


The values 
1196.04; fraction II, 
Values of heats 


The values of heat of polymerization 


are: fraction I, 8.83: fraction II, 9.04; fraction III, 9.37; and fraction IV, 10.13 kcal/mole 


of monomer. 


The standard deviations of the mean for heat of combustion are not more 


than 0.30 keal/mole, and for heat of polymerization, not more than 0.26 kcal/mole. The 


decrease in heat of polymerization with increasing molecular weight is discussed 


The low 


values of heats of polymerization are ascribed in part to steric interference between the 


substituent groups of the polymer chain 


I. Introduction 


This project was undertaken in connection with a 
general investigation of molecular properties of high 
polymers, and as a sequel to the work on heat of 
polymerization of styrene [1].'. The purpose was to 
determine the effect of the alpha-methyl group on 
heats of combustion and polymerization. It was 
known that alpha-methylstyrene (1-methyl-1-phenyl- 
ethylene) does not polymerize readily (in absence of 
catalyst) [2,3] and it was suspected that the steric 
interference between the substituent groups of the 
polymer chain might be partly responsible. A 
Fisher-Hirschfelder-Taylor molecular model of the 
polymer was made with the monomer units arranged 
head-to-tail, and alternating d and / units, that is, 
with Me and Ph groups interchanged, as follows: 


H Me H Ph H Me H Ph 
( ( ( ( C Cc ( C 
H Ph H Me H Ph H Me 


This model indicated that the molecule was so closely 
packed that it permitted only limited flexing or 
vibration and very little, if any, rotation. Ali 
models having other configurations (all head-to-tail) 
showed even higher degrees of interference. As a 
result of this steric effect, a low heat of polymeriza- 
lion Was anticipated [4]. A second purpose of this 
work Was to investigate the effect of molecular weight 
on the heat of polymerization. 


II. Methods and Apparatus 


lhe methods and bomb-calorimetric apparatus 
in this work were the same as those described 
iously [1, 5, 6 to 10). The observed value of 
of combustion in the bomb process, —Al’n, in 
experiment was reduced to the value of —AL’e®, 
decrease in intrinsic energy accompanying the 


1 brackets indicate the literature references at the end of this paper, 


combustion reaction 


92 
C,H,,(/ or +5 O.(q) OC (4) t 5H.A )(7/), 


with each of the reactants and products in its thermo- 
dynamic standard state. The standard state for 
C,H, in this equation is that of the liquid monomer 
or solid polymer. The procedure described by 
Washburn [11] for calculating the above reduction 
was modified so as to apply to 30° C and to the energy 
content of the gases at zero pressure. 


The values of heat of combustion are based on the 
amount of the combustion reaction, obtained from 
the mass of carbon dioxide formed (molecular weight 
44.010). No corrections were applied for unburned 
carbon in the crucible. In the few experiments in 
which carbon remained, it was found to be negligible 
in amount. Frequent tests, using a method de- 
scribed by Shepherd [12], revealed no carbon monoxide 
in the gaseous products of combustion 


The calorimetric system was calibrated with NBS 
Standard Sample 39f of benzoic acid, using the value 
for the heat of combustion of this sample under the 
conditions of the standard bomb process, QB 
26433.8 abs j/g at 25° C, obtained from the value 
reported previously [8, 13]. 

The energy unit used is the absolute joule. Con- 
version to the conventional thermochemical calorie 
was made by use of the relation [14] 1 cal=4.1840 
abs j. 


III. Materials 


The alpha-methylstyrene monomer and polymer 
were obtained from the Dow Chemical Co. through 
the courtesy of D. R. Stull. The monomer was 
stored in a closed bottle in a refrigerator at about 
11° C until it was introduced into glass bulbs pre- 
paratory to the combustion experiments. It con- 
tained no inhibitor of polymerization. 











was shut off, and observations of the pressure in the 
system were made with a McLeod gage over a period 
of about 3 weeks. The rise in pressure during this 
time corresponded to the evaporation of a few 
hundredths of a milligram per week of material, 
assumed to be methanol, from about 10 to 15 g¢ of 
polymer. 

After completion of the above treatment, the 
polymers were exposed to the atmosphere for about 
3 months. <A sample of about 3.85 g of Fraction II 
was then evacuated for 2 days at 50° to 65° C and 
weighed immediately after removal from the vacuum 
chamber. It was found to have lost about 2.0 mg 
all of which was regained upon exposure to the 
atmosphere for 6 days. It was concluded that all 
solvent and methanol probably had been removed 
from the polymer powders. Residual moisture 
would not affect the heat of combustion, as the 
results are calculated from the mass of CO, formed 
by burning the sample. 

The results of determinations of molecular weight, 
intrinsic viscosity, and oxygen content of the poly- 
alpha-methylstyrene fractions are 


Analysis of alpha-methylstvyrene by the Mass 
Spectrometry Section? of this Bureau gave the follow- 
ing results: 


99.6 mole percent 
0.4 


Alpha-methylstvrene 
Styrene 


The values given for heat of combustion of alpha- 
methylstvrene in table 2 have been corrected for the 
presence of this amount of styrene, using the value 
previously found [1] for the heat of combustion of 
styrene. 

The polymer used was formed during distillation 
of the monomer in the plant finishing still. No 
catalyst could be detected, although very minute 
amounts of some oxidation product in the iron still 
may have been present. The maximum tempera- 
ture to which the material was subjected was about 
50° to 60° C. After receipt at this Bureau, the 
sample of polymer was fractionated * and reprecipi- 
tated, the final product being four fractions con- 
sisting of fine nonecrystalline powders. In the frae- 
tionation process, the raw polymer was dissolved in 
distilled toluene to make a 1l-percent solution and 


filtered. The solution was placed in a_ 20-liter ' 
spherical flask in a constant-temperature (30° C) Fraction Molecular weight Intrinsic viscosity | Oxygen content | 
water bath. Methanol was added dropwise with ch 
stirring, and fractions were removed when the _™. Isj=dlig Avg. wt. % O: = 
methanol concentration reached 40 percent (Frac- + 3 — ©. 103 0 O59 re 
. 1 3, 750 065 . O47 pr 
tion 1), 60 percent (11), and 80 percent (111). The Ill 2 200 039 O85 du 
precipitates were allowed to settle for 1 to 3 days IV 1, 300 036 050 a 
and took the form of a viscous liquid phase at the ch 
bottom of the flask. Fraction IV was obtained by ve 
evaporation of the remaining solution. Fractions | The molecular weights' were determined by both 
I, II, and ILI were diluted to 1-percent solution with | ¢ryoscopic and isothermal distillation methods [15] 
a mixture of toluene and a small amount of methyl | for each fraction. Good agreement between the 
ethyl ketone to prevent phase separation from ab- | results of the two methods was obtained. In the 
sorbed water. The polymer fractions were repre- eryoscopic method 0.5- to 0.7-percent solutions in vi 
cipitated from these solutions, by adding each solu- spectroscopic grade cyclohexane were used, and the = 
tion dropwise to excess methanol in the approximate | Jowering of the melting point was about 0.02° to T 
ratio of 1 to 20. Fraction 1V was diluted to 4-pereent | 9.1° C. For the isothermal distillation, 5- to 15- _ 
solution and precipitated in the ratio of 1 to 10. | percent solutions in toluene were made and suspended 
The fractions were decanted, washed eight times | fy a glass dish in the vapor of the solvent at 30° C (2 
with methanol and dried at 30° C and 0.001 mm of | 4 modified method of measurement was used, in of 
mereury for about | month. The above treatment | which the rate at which the solution gains weight is 
of the polymers was such that monomer and lowest compared to the rate when a material of known 
molecular weight polymer probably were removed. | molecular weight is in solution. These methods give 
. ww 


The largest portion of the polymer yield was frac- 
tion II, lesser amounts were obtained for fractions 
Ill and LV, and fraction I was quite small. 

The vacuum pump system used for drying the 


approximate number average molecular weights 
They correspond to average degrees of polymeriza- 
tion of 45.7 to 11.0, which are not very large. Dif- 
ferences in heats of combustion per monomer unit 


polymer fractions contained a trap immersed in 
liquid air located in the line between the pump and 
the vacuum chamber, where the polymer powders 
were placed. By means of this trap, removal of 
condensible material from the polymer was con- 
tinued during periods when the pump was shut off 
(at night and on week ends). After evacuation had 
continued for about a month, the liquid air was re- 
moved from the trap, which was then pumped out 
to a pressure less than 107‘ mm of mercury. After 





for polymers having given differences in molecular 


weights are more likely to be measurable if the 
molecular weights are small. The low molecular 


weights of the polymer used here may be partly a 
result of the fact that the alpha-methylstvrene was 
polymerized at raised temperatures. Previous work- 
ers [2] have found it difficult to obtain polymers of 
alpha-methylstyrene of high molecular weight, except 
at very low temperatures. 

The measurements of intrinsic viscosity’ were 


this the stopcock between the pump and: the trap | 
? Under the direction of F. L. Mohler. 
? The fractionation was done under the™direction of 8. G 


made by means of the Ubbelohde viscometer at 


Brawley and 8. G. Weissberg 
McElwain and 8. G. Weissberg 


* Measured by D. J 


Weissberg. 5 Measured by J. W 
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mut 30° C, with concentrations in distilled toluene 
ving from about 0.5 to 2 g/100 ml. These vis- 
ty values are included only for relative com- 
ison. 
[he oxygen contents ® of the purified polymer 
tions were determined by a direct method 
developed at this Bureau [16]. In the procedure for 
vgen analysis, adsorbed water and air were re 
moved by sweeping the sample for 2 hours at 30° to 
50° C with pure helium. 

The effect of the oxygen content of the polymer 
samples on the values of heat of combustion and 
heat of polymerization will depend on whether the 
oxygen is present in some inert impurity or is chem- 
ically combined with the polymer molecules. In the 
former the effect of the oxvgen would be 
negligible. If the oxygen is chemically combined 
with the polymer, the values of heat of combustion 
will be too low, and the values of heat of polymeri- 
zation will be correspondingly too high. From the 
oxygen contents of the poly mer samples given above, 
it is estimated that if all the oxygen is combined 
with the polymer the effect would not be greater 
than 0.03 percent in heat of combustion, or about 
t percent in heat of polymerization. 

Since it is not known whether any of the oxygen is 
chemically combined with the polymer, the results 
have not been corrected for combined oxygen. The 
presence of the oxygen in the polymer samples intro- 
duces some uncertainty into the results, in addition 
to that resulting from other causes. Allowance for 
this has been made in calculating the estimated 


over-all standard deviations.’ 


IV. Results 


The results of the calibration experiments are 
given in table 1. Two series of calibration experi- 
ments were required because of repairs to the bomb. 
The results of heat of combustion measurements on 


cause 


ssured by F. W. McCulloch and W. H. Smith 
indard deviation of the mean as used in tables 1 and 2 is defined as 
z ! where d is the difference between a single observation and the 


ind nis the number of observations. In the case of the final values of heat 
istion given In table 3, the values of standard deviation include the effects 
ated systematic errors, as well as accidental errors. These values were 


ited from the relation 
Ql(se/I sq Q)?+ B+ R2 
ry » the standard deviation of the mean of the experiments with benzoic 
’ to determine FE, the energy equivalent of the calorimetric system; s@ is the 


rd deviation of the mean of the results of the experiments to determine Q, 
combustion of the sample; B is an allowance of 510-5 for the standard 
ition of the value used for the heat of combustion of benzoic acid: and R is 


wance of 2X10~ in the case of the polymer or 5X10-5 for the monomer 
ns, and takes account of systematic errors associated with the impurities 
wlytmer samples and the determination of the amount of the combustion 
rom the mass of carbon dioxide. The expression used in connection 
vaiues in table 4 was 
s’ =[(#m)?+(#,)2] 


ssion &= and s» were calculated from the data of the heat of combus 





ements on monomer and polymer, respectively, by means of expres 
Sm = Q[(ee/ F)2+(8¢@/Q)2)'/2 
Sp= Ql (sx/ F)2+ (se Q)?+ R?) 


lL here R is an allowance of 15 
e polymers. It 
com bustion 


x 10-5 for the errors associated with impurities 
is assurned that systematic errors in the value used for the 
of benzoic acid and in the method of determining the amount 

“om bustion reaction will cancel in taking the difference between the heats 
nbustion of polymer and monomer 
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the polymer and monomer for each experiment were 
calculated by using the appropriate mean of cali- 
bration series A or B, depending on whether the 
combustion experiment was performed before or after 
the alteration of the bomb. 


PaBie 1. Calibration of calorimeter with benzoic acid 
Observed energy equivalent (2) at 28.5°C 
Series A Series B 
ahs j/°C ths j/°C 
13743. 2 13780. ¢ 
13792. 6 L373. 8 
13704. 4 S782. 8 
13792. 0 IS781. ¢ 
13791. 7 13783. ¢ 
13702 1788. 0 
13780. 8 13786. 4 
13791 13785. 1 
13783. 0 
13787. 4 
LATS. 
Mean 13792. 2 13784. 4 
Standard deviation of mean *, sz 0.5 0.7 


® See footnote 7 


The results of heat of combustion measurements 
are given in table 2, where—Al’s is the observed 
heat of combustion under the conditions of the bomb 
process: the temperature (30° C) to which the reac- 
tion is referred, the mass (1g) of water placed in the 
bomb before each experiment, the internal volume 
of the bomb (377 ml), and the conditions specified in 
columns 1 and 2 of the table. The quantity —Al’e¢ 
is the reduced value, as explained in section IT. 

The values obtained for the mass of carbon dioxide 
formed in combustion are lower than the correspond- 
ing values calculated stoichiometrically from the 
masses of samples burned, assuming that the compo- 
sition of the samples is represented by (CyHjo)-. 
The average differences between observed and cal- 
culated masses of carbon dioxide are, for alpha- 
methylstvrene 0.052 percent, and for poly-alpha- 
methylstyrene fractions: 1, 0.110 percent; I1, 0.106 
percent; II], 0.188 percent; IV, 0.082 percent. 
Frequent measurements of the carbon dioxide 
formed in combustion of Standard Sample 39f of 
benzoic acid vielded results in agreement with calcu- 
lated values within about 0.016 percent. The 
observed deficiency of carbon dioxide in the case of 
alpha-methylstyrene can probably be attributed to 
air or moisture dissolved in the liquid. The defi- 
ciency in the case of the samples of poly-alpha- 
methylstvrene is probably due to the presence of 
nonhydrocarbon material in the samples as an im- 
purity. 

In table 3 are given values of —AUc° per mole of 
C,H,» (or per monomer unit) derived from the mean 
values given in table 2 by taking 9 moles of carbon 
dioxide as equivalent to 1 mole ef monomer, or the 
C,H,» unit of the polymer. There are also given in 
table 3 values of —AHc°, the decrease in heat con- 
tent in the combustion reaction, with all reactants 
and products in their standard states. 








TaBLe 2 





Romb-calori metric 


measurements on 


poly-alpha- 


methylstyrene and alpha-methylstyrene 


Mass of 
sam ple 


POLY 


g 
0. S855 
U5498 
vI254 
NOB T 
GSU 
Vsil4 


Mean 


Standard deviation of 


pres 


ure at 
sd CO 


ALPHA-METHYLSTYRENE 


wi.2 + OO141 
tt + 1vAe 
7.2 +. 418 
Se 3. 24617 
8 ‘. 14561 
2 +. 28476 


fmean *, #@ 


Mass of 


Heat of combustion at 
0) ; 


FRACTION I 


abs jig COsg ahs jig COs 

12628. 3 12621. 7 
12625 3 12618. 4 
12620. 7 12613. 7 
12620. 5 12613.9 
1216.8 12610.3 
1220. 2 12613.7 
12615. 3 

1.7 


POLY-ALPHA-METHYLSTYRENE, FRACTION I 


Oo w2i70 
“STH2 
O55 2 
USS] 
4285 
anus 


Mean 
Standard deviatior 


POLY-ALPH 


0. G2SS86 
525 
GOSI6H 
7 


“7505 


Mean 


Standard deviation of 


POLY-ALPH 


1. (2726 
0. G4ASA7 
US42S8 
WHi2 
4sn06 
M128 


Mean 
Standard deviation « 


ALPHA 


0. SY20D8 
seat 





6.4 +. OGG08 
wo 1 +, 20856 
“7 + 1080. 
AH. 7 W407 
i 15585 
4.2 24007 
of mean *, #@ 


128620. 5 12613. 7 
12620 12613. 6 
12619. 4 1212.5 
1217.5 12610. 5 
1223. 3 12616. 4 
12618. 8 12611.8 
12613. 1 

os 


A-METHYLSTYRENE, FRACTION II 


t + 101 
4 +. 1200 
5 3. 1378 
2 3. 2160 
2 + 2A540 
2.4 +. 15006 





mean *, &q@ 


2616. 5 12610. 2 
12616. 1 12609. 8 
12620. 5 12614. 2 
12613. 5 12607. 1 
1215.5 12H09. 1 
113.2 12606. 9 
12600. 6 

1.1 


A-METHYLSTYRENE, FRACTION IV 


2.9 44065 





{mean *, #@ 


METHYLSTYRENE, MONOMER 


12.7 2 GSATS 
i2. 6 2. 71183 
12.6 2. 65744 
t2 2. 4808) 
2.2 + 4100 


2.3 +. (2352 


Mean 
Standard deviation of mean *, #@ 
* See footnote 7 

Quantity 


4te? GC 
BV EL we 
SHe® (25° C 


Standard deviation of 


* See footnote 7 


mean *, # 


1M10.7 12M. 2 
12408. 1 11.7 
107.1 12800. 6 
12007. 6 ww. 2 
12808. 6 1202. 2 
12605. 9 12599. 6 
12601.6 
io 
12710. 2 12710. 4 
12708. 4 12708. 7 
12700. 4 12709. 7 
12705. 7 12706. 2 
12710.9 12711.2 
127109. 3 2700 
12709. 3 
0.7 


TABLE 3. 


Poly 


Fraction I 


abs kj kcal abs kj 
CoH CoH CoH 
unit unit unit 
4996. 80 4995. 93 
MOOS. 10 Mw2. ZB 
UMM. 24 1196. 04 MWS. 37 
1. 24 0. 30 1.10 


Fraction II 


For reducing the values at 30° to correspondinny 
values at 25°C, the heat capacity of alpha-meth!- 
styrene was derived from unpublished data obtaino<d 
at the Bureau on the heat capacity of styrene [17] 
and an estimate of the increase in ©, accompanying 
the addition of a CH, group; the heat capacity of 
poly-alpha-methylstyrene was obtained from that 
of the monomer by use of the ratio of heat capacities 
[17], polystyrene: styrene. 

Values of heat of polymerization of alpha-methy|- 
styrene to give the four polymer fractions are shown 
in table 4. These values were obtained by taking 
the difference between the value of —AHe® (25° C 
(table 3) for the monomer, and the value for each of 
the polymer fractions. 


Tarte 4. Values of heat of polymerization 
Standard 
Fraction SIT» (25° ¢ deviation of 
I mean *, & 
keal) mole 
I SSS 0. 26 
Il ¥. 04 22 
iil 9. 37 23 
I\ 10.13 22 
* See footnote 7 


V. Discussion 


The values of heat of polymerization given in 
table 4 are lower than the value (—AHp°®=16.7 
keal/monomer unit) previously reported for styrene 
[1], and are also lower than any value previously 
reported for a vinyl compound [18]. It has been 
suggested by Flory [19] and by Evans and Tyrrall 
[4] that low values of heat of polymerization of 
1,1-disubstituted ethylenes are to be attributed to 
steric interference between the substituent groups 
There are other effects that may be partly respon- 
sible for the low heat of polymerization, such as the 
effects of side groups on bond energies in the mon- 
omer as compared with those in the polymer, and 
the lack of a large “end effect’, which in ethylene 
arises from the nearness of the double bond to the 
ends of the monomer molecule [18]. It is believed, 
however, that the magnitude of the latter two 
effects is not sufficient to account for the very low 
values of —A//p° for alpha-methylstyrene, but that 
steric interference must be an important factor. 


Values of heat of combustion 


ilpha-methylstyrene (solid 


Alpha-methylst yrene 


liquid 
Fraction III Fraction I\ 

keal abs kj keal ibs kj keal . 
CoH CoH CoH CoH CyH oo a 
unit unit unit unit unit 

1904. 4 4901. 37 5034. 08 

5000. 84 4907. 67 M40. 33 
1195. 83 KOOL. OS 1195, 50 40908. 81 1194.74 M1. 17 124. 87 

0. 1.15 0. 27 1.09 0. 26 0.53 0.13 
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In general, it is found that molecular models of 
vmers of 1,1-disubstituted ethylenes indicate that 
ic interference is quite pronounced if the mon- 
er units are joined head-to-tail but is much less 
the head-to-head and tail-to-tail arrangement. 
he low values of —A//p° for alpha-methylstyrene 
therefore consistent with other evidence {20, 21) 
at polymers of 1,1-disubstituted ethylenes are 
wmed by head-to-tail addition. 


It will also be seen from table 4 that the values of 
heat of polymerization of alpha-methylstyrene de- 
rease systematically with increasing molecular 
weight of the polymer. This was not expected, and 
t seems worth while to examine other pertinent 
nformation to see whether it is consistent with the 
observed variation of heat of polymerization with 
molecular weight of polymer. 


It is well known [22] that the heats of formation of 
an homologous series of hydrocarbons, either liquid 
or gaseous, can be represented by a linear function of 
ihe number of carbon atoms in the normal alkyl 
radical, provided this number is fairly large, but that 
this relation does not hold for the first few members 
of the series. It may be inferred that the heat of 
formation of a straight-chain polymer molecule will 
be a linear function of the number of monomer groups 
in the molecule, provided this number is sufficiently 


large, but that this relation may not hold for a 
molecule containing only two or three monomer 


units Then the heat of addition of a mole of liquid 
or gaseous) monomer to a mole of straight-chain 
liquid or nonerystalline solid (or gaseous) polymer 
according to the reaction 


Mn-1tm=—mM,, (1) 


will be a constant, g, independent of n if n > r, where 
rissome small integer. Ifn S r the heat of reaction 
will, in general, differ from q. Let (q¢+A;) represent 
the heat of the reaction 


m +n m i Sr). (2) 
hen the heat of the reaction 
nim Ihe, (3) 
with n > r will be given by 
AH (n—1)q + SUA (4) 


i 


nd the heat of polymerization per mole of monomer 


ll be 


DAi— 


AH p° AH®/n 


it 


q+‘ ° (.)) 


monomer and polymer, respectively, then eq 5 may 


be written in the form 
M, [ «| 


M 


(6) 


| If M, and M represent the molecular weights of 
| 

| 

| AT p?=4 + 


AHp 


i=2 


From eq 6 it is evident that if SSA;<q, 

| 

will increase with increasing .M, as has been shown 

to be true in the case of ethylene [23]. In order that 
AHp° shall decrease with increasing M it is neces- 


sary that SA, >y. 


i=2 


The equation 


219- 
8.424+ J (7) 


was derived by the method of least squares from the 
experimental values of —A/7p° and M for the four 
fractions of alpha-methylstyrene polymer. The 
observed values of AHp® and values calculated 
from eq 7 are compared in the following tabulation. 


Alp 
V 
. Observed 
Observer Calculate ' 
bserved alculated Caleuleted 
kcal, monomer kcal monomer kcal) monomer 
unit unil uni 
5400 S83 & S83 0. 00 
3750 9.04 9. 01 03 
2200 9. 37 9, 42 . 05 
1300 - 10. 13 10. 11 02 


The maximum difference between observed and 
calculated values is seen to be 0.05 keal, which corre- 
sponds to about 0.004 percent of the heat of combus- 
tion. The standard deviation for qd is 0.05, and for 


M334, 


q)is 97. According to eq 7, q=84 keal 


r 


and (SoA q)=18.6 keal, so that >A 


27 0 keal. 


| As indicated in the derivation, eq. 5 and 6 are valid 
only for n=r, that is, for values of n such that heat 
of formation per mole of polymer is a linear function 
ofn. The good agreement between the observed and 
calculated values of —AHp° therefore indicates that 
heat of formation (or combustion) of poly-alpha- 
methylstyrene per mole of polymer is a linear func- 
tion of n, the number of monomer units in the 
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molecule,” at least for n211, and it seems probable 
that this relation would hold for n>3 or 4. 

It may be noted that in the derivation of eq 5 
and 6 it was tacitly assumed that » and V are con- 
stants for any given sample of material, that is, that 
all of the polymer molecules in a given sample are of 
the same length. When these equations are applied 
to polymers containing molecules of various lengths, 
it can be shown that the proper values of VW and n 
to be used are the number averages, which were 
actually used in the derivation of eq 7. 

The only 1,1-disubstituted ethylene for which data 
are available on which a reasonable estimate of 


> A, for that compound can be based is isobutene. 
t=2 
Chemical and X-ray evidence |[20, 21], as well as a 
low heat of polymerization, indicate that the 
monomer units in the polyisobutene molecule are 
joined head-to-tail. Evans and Tyrrall [4] have re- 
ported for the heat of polymerization of liquid iso- 
butene to liquid polymer the value 12.6 keal per 
mole of monomer, and for the heat of dimerization 
of g monomer to head-to-tail dimer 
(2,4.4-trimethyl-1-pentene) the value 19.9 keel per 
mole of dimer. These values indicate that for 
butene g=12.6 keal and A,=7.3 keal, if the difference 
between heats of vaporization of polymer and 
monomer is neglected 

The heat of dimerization of alpha-methylsty rene 
is not known. If it is assumed to be about the same 
as that of isobutene, that is, about 20 keal per mole 
of dimer, and if g is about 8 keal as indicated by eq 7, 
then A, would be about 12 keal. It would ‘be eX- 
pec ‘ed a the fact that steric interference is —_ 
small in a head-to-tail dimer, but is pronounced i 
the thd and in polymers of higher oe 
weight [4], that (A,+.A,-+ A,) would be small in 
comparison with A,, and might even be negative, 


FASCOUS rauscous 


1IsoO- 


so 


r 
value of >) A, based on the 

s=2 
above assumptions would be much more than half of 
the value 27.0 keal derived from eq 7. However, the 
fact that the above value of A,, 12 keal, is greater 
than the value 8 keal derived from eq 7 for q, suggests 
that the above assumptions are consistent with a 
decrease in heat of polymerization with increasing 
molecular weight of polymer, although they do not 
lead to quantitative agreement with the observed 
variation of heat of polymerization with molecular 
weight. In this connection, it should be emphasized 


that it is unlikely that the 


*It has been pointed out [24] 
representing the heat of combustio 


that this fact can be shown more directly by 
n per mole of polymer as a linear function of n 


ASHle® (25°C) (solid alpha-methylstyrene polymer i+Bn, 
or heat of combustion per monomer unit by the equivalent relation 
Alk 


_ ‘ 
25°C solid alpha-methylstyrene polymer)=— +B 
n 


n 
determination of the cons —_ in ee latter equation by the method of least 
hk s yielded the values A Mi, B= 1196.447 Combination of the above 
equation with the value for the bea at of combustion of alpha-methylstyrene 
monomer, designated by the symbol C, yields the relation 


i 
—Allp® = ¢ ; RB, 
r 
which is of the same form as eq 5. However, it is not immediately evident from 
this derivation how the constant A is related to g=(C— 8B) and to the end effects, 
Ay, Aa, This relation is given explicitly by the derivation in the text 
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that some of the assumptions made in comparing 
alpha-methylstyrene with isobutene are subject to 
considerable uncertainty. 

The estimated effect of the oxygen contents of thy 
polymer fractions on heat of combustion is too small 
to account for the observed decrease in heat of poly- 
merization with increasing molecular weight. 

Although the low values of heat of polymerization 
of alpha-methylstyrene indicate that the monomer 
units in the polymer molecule are joined head-to-tail 
for the most part, there is some reason for conjectur- 
ing that the first step in the formation of the polyme: 
molecule may be the tail-to-tail 


R, R, 
C—CH,—C H,--(¢ 
R R, 
combination of two monomer units, and that the 


subsequent growth of the molecule is by head-to-tail 
addition. As stated previously, a model of a polymer 
molecule with monomer units joined head-to-tail 
indicates a considerable amount steric hindrance 
A model of a polymer molecule in which two units 
are joined tail-to-tail, and the shen ‘rs are joined head- 
to-tail to one of the first two indicates only a little 
steric hindrance in the first three groups. It would 
therefore be expected that for such a molecule both 
4, and A; would have fairly large positive values 
r 


possibly large enough to bring >} A; into quantitative 


i=2 
agreement with the value derived from eq 7. 

In regard to whether the initial step in the forma- 
tion of a 1,1-disubstituted ethylene polymer molecule 
could be the tail-to-tail combination of two monomer 
units, it may be noted that Evans and Polanyi [25j 
have concluded from calculated heats of reaction that 
in the first step, tail-to-tail addition of isobutene is 
more probable than head-to-tail addition, and that 
the latter is more probable than head-to-head 
addition. 

It is believed that the heat of polymerization of 
alpha-methylstyrene decreases with increasing molec- 
ular weight of polymer, as is indicated by the experi- 
mental data reported herein. The magnitude of the 
observed decrease is greater than would be expected 
from the data on isobutene. The foregoing discus- 
sion shows, nevertheless, that the observed change in 
heat of polymerization with molecular weight is con- 
sistent with certain plausible assumptions regarding 
the structure of the polymer molecule and the heats 
of the first few steps of the polymerization reaction 


VI. Previous Work 


Auwers, Roth, and Eisenlohr [26] reported for the 
heat of combustion of alpha-methylstyrene monomer: 
the experimental values 10,183 and 10,194, mea 
10,189 cal/g. By applying appropriate corrections 
for the value used for heat of combustion of benzoi 
acid, and by applying the Washburn correction, cal 
culating for weight in vacuum, and correcting to th: 
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revction at constant pressure at 25° C, the mean of 
the above results may be reduced to a comparable 
havis with the present work. After these corrections, 
{ above value for heat of combustion becomes 
1202.8 keal/mole 

\pproximate corrections applied to the value re- 
ported by Lemoult [27] give 1209.4 kcal/mole for heat 
of combustion of the monomer. 

Previous measurements at this Bureau [28] give 
the value for heat of combustion of liquid alpha- 
methylstyrene AHe® (25° C)=1204.3 keal/mole, 
with an uncertainty of 


») 


+ (0.4 keal/mole. 
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The McKee Worker-Consistometer with Constant-Speed 
Drives 


Samuel A. McKee and Hobart S. White 


Two constant-speed drives have been developed for use with the McKee worker-consist- 
ometer for mechanically working a fluid and measuring its flow characteristics in the same 


series of operations. 


Previously, weights were used to provide various constant loads at which 
measurements were made of the rate of flow of the material under test. 


The modifications 


provide for measurements of the resultant forces when operating at constant rates of flow. 
They permit the determinatioa of the effect of wotking at various constant rates of shear and 
also provide for a much greater range of loads. 
rest data are given showing the effect of mechanical working upon the flow characteristics 
of some lubricating greases and some raw synthetic rubbers. 
The results show that the apparatus provides a tool useful in the measurement of the 
appareat viscosity and the effect of mechanical working for non-Newtonian materials having 


thixotropic properties. 


The performance of the apparatus is adequate to cover consistencies 


of materials ranging from a light lubricating grease to 100 percent of raw rubber. 


I. Introduction 


The development of a worker-consistometer for 
investigating the effect of mechanical working on 
the physical properties of lubricating greases has 
been described in a previous publication.’ With 
this apparatus a grease or similar material may be 
worked and its flow characteristics measured in the 
same series of operations. It consists essentially of 
two coaxial steel cvlinders and mating pistons with 
a capillary shearing element between them. By 
means of the closely fitting pistons the material 
being tested is forced repeatedly through the shearing 
element from one cylinder to the other. The 
shearing elements were made by drilling holes in 
steel disks 4 in. thick, using a No. 79 drill. With 
the original apparatus, weights acting on the top 
piston provide a constant load during the downward 
pass when the rate of flow is being measured; an 
air evlinder with leather-cup piston and rod acts 
against the lower piston during the upward pass. 

With materials such as greases and some rubber 
solutions, the rate of shear changes with working 
when a constant load used. Recent work has 
been directed toward increasing the usefulness of the 
worker-consistometer by modifications that (1) per- 
mit operation at constant nominal rates of shear and 
extend the range to higher viscosities. The 
modified instrument is suitable for studies of the 
rheological properties of rubber and rubber solu- 
tions, as well as greases and similar substances. 

The modifications required were concerned cheifly 
with the driving mechanism and with the method 
for measuring the forces involved. Two constant- 
speed drives have been constructed: one is a cam- 
operated mechanism with a load range up to 200 
lb; the other is a screw-operated drive with a capacity 
of 1,600 lb. The operation of each mechanism is 
automatic, and is such that readings of the forces 
may be taken at desired intervals during a run of 
continual working. 


Is 


(2) 





A. McKee and H. 8. White, A worker-consistometer for lubricating greases, 
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This paper presents details on the construction, 
operation, and range of the new driving mechanisms, 
together with examples of the application of the 
modified worker-consistometer to four lubricating 
greases and five synthetic rubbers. 


II. Apparatus 


The design of the major elements of the worker- 
consistometer has been retained. This has been 
described in the previous publication (see footnote! 
Each worker unit consists of two coaxial steel eylin- 
ders and mating pistons with a capillary type shear- 
ing element between them. Relative motion be- 
tween pistons and cylinders forces the material under 
test through the shearing element. Views of two 
worker units are shown in figure 1. A set of unas- 
sembled parts is on the left, and an assembled unit is 
on the right. Four such units are now available for 
use. Each is provided with 1-, 10-, and 50-hole 
shearing elements. These shearing elements are 
made by drilling holes in circular steel disks \ in. 
thick, using a No. 79 drill, giving capillaries \ in 
long and about 0.015 in. in diameter. 


1. Cam-Operated Machine 


With the cam-driven apparatus, the pistons are 
moved back and forth by a motor-driven cam- 
operated mechanism, and determinations of the 
pressure developed in pushing the material under 
test through the capillaries in the shearing element 
are made by measuring the force sequired to hold 
the cylinders in a fixed position. Views of this 


apparatus are given in figures 2, 3, and 4. The 
worker unit (shown in fig. 1) is enclosed in the 


duralumin housing, A, (fig. 2), which is encircled by 
the electric heater, B. A mercury-in-glass thermal 
regulator placed in a perforated tube is used to con- 
trol the test temperature, which is measured by a 
thermocouple contacting the shearing element. The 
flanges at the ends of the housing are attached to 
flexure plates, which provide “free” action along the 
vertical axis. The weight of the housing and parts 
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ported by it is counter-balanced by the weight, 
fic. 4) mounted on a lever arm. This arm is 
vided with a pointer and scale, the zero position 
which indicates when the system is in balance 

h no deflection in the flexure plates. 

lhe operating voke, £, (fig. 2) guided by bearings, 

s actuated by the cam, //7, and moves the pistons 

ck and forth, a distance of approximately 2 in., by 
eans of the push rods, A. During the upward 
ss, the force transmitted to the cylinder assembly 
by the action of the piston in pushing the material 
inder test through the shearing element is measured 
by the weighing seaic, which is connected to housing 
A by rods acting through an 8 to 1 lever arm. The 
arrangement is such that the load range is from 0 to 
120 Ib with weight DP (fig. 4) in place and from 80 
to 200 lb when J) is removed. A turnbuckle on the 
rod from the scale permits adjustment to compen- 
sate for the scale travel, so as to keep housing A in 
its zero (balanced) position when reading the scale. 
Movement of the housing is limited by adjustable 
stops, and on the downward pass the force is taken 
by the lower stop. A counter actuated by each 
downward stroke of the voke records the revolutions 
of the cam 

The pointer, ZL, mounted on voke / moves along 
a steel measuring scale and is used for measuring the 
speed of the pistons. 

Cam /7 is mounted on a shaft supported by two 
ball bearings. The shaft is connected to a positive- 
action variable-speed transmission. A 36-rpm gear 
motor supplies the power to the transmission through 
a four-step combination of V-pulleys and a 30 to 1 
reduction rear. The gear motor is mounted on a 
hinged base so that the V-belt can be easily shifted 
to a different set of grooves. This driving mecha- 
nism in conjunction with the use of the 1-, 10-, or 
\0-hole disks provides a continuous range of nominal 
rates of shear from 200 to 100,000 reciprocal seconds, 
the latter value being limited, not by the drive, but 
by the length of time required for making an accu- 
rate load observation. The range of rates of shear 
covered with any given material under test is, of 
course, also limited by the load range of 1 to 200 Ib. 
For a broad range of rate of shear this load range 
corresponds to a viscosity range from 1 to 100 
and for a relatively low rate of shear of 200 
reciprocal seconds, the viscosity range is from 20 
to 4,000 poises. , 


Poses ; 


2. Screw-Operated Machine 


With the serew-operated apparatus the evlinders 
of the worker unit are moved back and forth, and 
the force measured is that required to hold the pis- 
ons ina fixed position. This apparatus is shown in 
. 6, and 7. The worker unit is enclosed in 


vures 5 


steel housing with a thermo couple and tempera- 
ire control as described for the cam-driven appa- 
The housing flanges are attached to a steel 
perating-voke, M (fig. 5), which is mounted in three 
eve bearings attached to a heavy steel plate 
ounted in a vertical position. 


itus 


This mounting is 
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such that the yoke may be moved up and down, 
but rotation is prevented. 

Fastened to the lower end of the yoke is a bronze 
screw, N, which engages a steel nut, ?, mounted 
on two ball bearings housed in a steel box, R. One 
of the interchangeable gear trains that connect the 
nut to a reversible gear motor is shown in figure 5. 
Adjustable screws on the voke contact push buttons 
that operate magnetic switches that automatically 
reverse the motor at the end of each upward and 
downward stroke of the voke. These screws and 
push buttons are adjusted for a stroke of about 
2in. A counter is actuated at the end of each down- 
stroke of the voke. 

During the reciprocating motion of the cylinders 
the pistons are restrained from moving by loading 
pins, S, acting against calibrated steel springs, 7. 
A dial indicator measures the deflection of the canti- 
lever arm of the spring and hence the force neces- 
sary to push the material under test through the 
capillaries in the shearing element. Five pairs of 
springs of different capacity cover the range of loads 
up to 1,600 lbs. The gear trains provide a number 
of rates of shear ranging from less than 1 to 100,000 
reciprocal seconds. The range of viscosity is eight 
times greater than that of the cam mechanism for 
its range of rate of shear (200 to 100,000 sec™') 
With the screw-operated apparatus the range of 
viscosity determinations is extended to 600,000 
poises at 10 recriprocal seconds or 6,000,000 poises 
at 1 reciprocal second. 


3. Accessories 


In testing samples that cannot be readily poured, 
the filling apparatus shown at the left in figure 8 was 
used. It consists of a split evlinder, piston clamping 
ring, and positioning ring. After the two halves of 
the evlinder are filled with the sample, they are as- 
sembled and positioned against the end of a worker 
cylinder, and by means of the filling piston the 
sample is transferred to the worker evlinder. When 
the worker is assembled the worker pistons are coated 
with pew sample to insure that the clearance 
space is filled, and care is taken to prevent trapping 
air with the sample. 

For the removal of air from spongy materials such 
as rubber, steel tubes with loose-fitting steel pistons 
are used. After small pieces of the sample are forced 
into the steel tube, a piston is started into each end, 
and a load of a few hundred pounds is placed on the 
pistons for several hours. The plug-shaped sample 
thus formed is pressed out of the tube, assembled in 
the split evlinder, and transferred to the worker. 
When using rubber samples, a piece of the rubber 
moistened with benzol is rubbed against the worker 
pistons until a coating is formed, and the coated 
pistons are pushed through the cvlinders a few times 
until both are coated. The benzol is then allowed 
to evaporate before assembly. 

In cases where it is desirable to store a sample in- 
tact for future test, a special cylinder is used. A 
storage cvlinder with its piston, cap, positioning ring, 








and push rods is shown on the right in figure 8. A 
storage-cylinder assembly containing a stored sample 
is shown also. The push rod having two steps is 
used to push the worker piston when transferring a 
sample to a storage evlinder, and the push rod with 
one step is used to push on the storage piston when 
transferring the sample back to the worker. 


III. Calibration 


The shearing elements were calibrated in their re- 
spective cylinder assemblies with oils of known vis- 
cosity in the manner described in the previous paper 
(see footnote 1). As these oils are Newtonian fluids, 
the flow may be represented by the equation, 


rP 
S “2L K,W 
1 "R 4Vi K,Vie 
Ti 


where 
n —absolute viscosity, poises 
S=shearing stress at the wall of the capillary, 
dynes/em? 
?=rate of shear at wall, sec 
P=pressure drop across the capillary, 
cn’ 
L=length of capillary, «em 
r=radius of capillary, em 
V/t=volume rate of flow per capillary, em*/sec 
W’=load on piston, |b 
K.,=shearing stress constant 
K,=rate of shear constant. 


-1 


dy hes 


The piston friction for each assembly was also de- 
termined by operating with an oil of known viscosity 
with the shearing element removed from the unit. 
The values obtained were used as a correction to As, 
the shearing stress constant and also for determining 
the average diametral clearance between the cvlin- 
ders and mating pistons. 

The values obtained in the calibrations of the re- 
spective worker-units and shearing elemens are given 
in table 


Taste l. Calibration constants of worker units and shearing 
elements 

Work Average Number Average 

oe ~‘ diametral Disk o diameter Kr Ks 
m clearance holes of holes 

m m 

0. 00015 1 1 0. 01508 183, 000 5, O55 

| 41 1 O1L506 182, O00 4, ORD 

4 O00g8 42? w O146y8 106, 000 4,710 

| 43 ) Oo1no 200, 000 4,135 

| | l OL487 189, 000 4, 0 

0008 | 2 10 O1453 2, 000 4,700 

| NM Mw 01477 193, 000 4, 530 


For convenience in measuring a very wide range 
of consistencies, the different units have different 
clearances between pistons and cylinders. In gen- 
eral, unit 3 with disk 31 is used with low-viscosity 
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samples, as the small clearance minimizes leakag» 
and the higher piston friction is still a relatively smal] 
correction “when using the I-hole disk. Unit 5 with 
disk 53 is used for the very high viscosity samples 
When a 50-hole disk is used, the larger clearance of 
unit 5 is desirable to reduce the piston-friction cor- 
rection. Unit 4, which is a compromise between 
leakage and piston friction, is usually used for samples 
of intermediate consiste nev. W ye this unit the 
piston-friction corrections are 0.4, 3.4, and 14.7 per- 
cent for disks 41, 42, and 43, respectively. 

Some indic ation of the accuracy obtainable with 
the apparatus is given in figure 9, where viscosity 
data are plotted against rate of shear. Sample N-15 
is a Newtonian oil, and the horizontal lines represent 
the viscosities at the temperatures indicated obtained 
with a modified Ostwald viscometer, whereas the 
points represent the viscosities obtained at the same 
temperatures, using worker-unit 4. When using this 
unit for viscosities below poise, a greater spread of 
the data would be expected because of the low loads 
involved. 

Data are also given in figure 9 for a non-Newtonian 
solution of GR-I rubber in a light mineral oil ob- 
tained with worker-unit 4 with all three disks. In 
order to cover the same range of rates of shear with 
disks 43 (50 holes) and 41 (1 hole), some runs with 
the latter were made by applying small dead-weight 
loads to the upper piston. These are indicated by 
the crosses. The good general agreement of the data 
obtained with all three disks is of particular interest 
The apparatus was calibrated with Newtonian fluids, 
and as the rates of shear at capillary and piston Ww a 
are not the same, some error would be expected in 
applying piston-friction corrections when sethinn 
with a non-Newtonian material. However, the data 
indicate that in this particular case the error is quite 


small. 
IV. Test Results 


The test data are given in this paper primarily to 
illustrate the application of the apparatus to certain 
investigations. They are not of sufficient scope to 
provide specific relations between the rheological 
properties of the materials and other properties or 
factors involved. 

In the previous publication (see footnote 1) results 
were given on the effect of mechanical working on the 


flow characteristics of various lubricating greases 
when operating in the worker-consistometer at 


Under such conditions the 
the grease was subjected 


various constant loads. 
rates of shear to which 


varied throughout the test. Flow measurements 
were made with two of these greases, grease A 


(lithium soap) and grease D (sodium soap), when 
operating at 100° F at various constant rates of 
shear. The results are given in figures 10 and 11, 
where the apparent viscosity (plotted on a log scale 
for convenience) is plotted against the number of 
passes through the shearing element. In each test 
at a given rate of shear, flow measurements were 
taken at intervals during a period of mechanical 
working of 200 passes, after which the sample was 
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ferred to a storage cylinder. After 4 weeks of 
we at room temperature, flow measurements 
made with each sample during a working period 

) passes. ‘They were then subjected to an addi- 
‘| storage period of 35 weeks and another 26- 
working period 

With both of these vreases the effect of mechanical 

king in lowering the viscosity was quite marked 
during the first 200 passes. With each grease this 
change in viscosity is roughly proportional at all 
rates of shear over the range covered. The ratio of 
the viscosity at the 200th pass to the viscosity at the 
first pass was approximately 0.6 for grease A and 0.5 
for grease D. 

The data taken at the first pass after a storage 
period provides an indication of the so-called healing 
process. In general there was a measurable increase 
in viscosity after storage. Apparently grease A had 
become relatively stable (with respect to mechanical 
working) after the first working period, as subsequent 
working after storage periods did not lower the 
viscosity much below that at the 200th pass, whereas 
with grease D there was a definite lowering in vis- 
cosity bevond the 200 pass 

Figure 12 shows the relationship between apparent 
viscosity and nominal rate of shear for both greases 
at various degrees of working. The points for 1, 10, 
and 100 passes were obtained from the curves in 
figures 10 and 11. Over the range covered with a 
given grease and a given pass, a linear relationship 
was obtained on the logarithmic plet. 

Flow tests were made also on four other greases: 
(soda-lime soap) and greases G, H, and J 
soda soap At a rate of shear of 18,100 reciprocal 
seconds and test temperatures of 100° and 245° F, 
apparent were determined at intervals 
during 226 passes with a 4-week storage period after 
the 200th pass. The data obtained are given in figures 
13 and 14, where the viscosity (on a logarithmic 
plotted against the number of 
Figure 13 indicates that with grease G the increase 
in temperature from 100° to 245° F caused a very 
marked reduction in the viscosity of the grease in 
both the unworked and worked condition. On the 
other hand, when grease F was worked at the higher 
temperature its viscosity increased to such an extent 
that at the 200th pass there was not much difference 
v at the two test temperatures. A similar 
trend of lesser degree is shown also by grease H in 
igure 14. a. 

lhe serew-operated mechanism was designed pri- 
irily for an investigation of the rheological proper- 

s of various synthetic rubbers and solutions con- 
ning high percentages of rubber. The range of 
peration of the apparatus is such that it is capable 
measuring viscosities as high as several million 
ses at low rates of shear and as low as a poise at 

h rates of shear. Accordingly, it may be used for 
termining the viscosities of rubber solutions 
ving from 10 to 100 percent of rubber. 


grease F 


viscosities 


scale Is passes 


Ih Viscosity 


nvestigation is being conducted in cooperation with the Office of Rubber 


Reconstruction Finance Corp 
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Data showing the effect. of working on the ISCOsity 
of five GR-S rubbers are given in figure 15. Samples 
X-418 and X-—478 are rubbers in commercial pro- 
duction, and 1, 2, and 3 are special experimental 
samples. In these tests, 100 passes were made with 
each sample at a rate of shear of 99 reciprocal 
seconds and at a temperature of 212° F. In each 
case there was an increase in viscosity during the 
first few passes, after which the behavior was different 
for different samples. It is of interest that whereas 
samples 2 and 3 have the same Mooney numbers, 
their apparent viscosities differed considerably under 
the present test conditions. The deviations of the 
points from the curves for samples 1 and X—418 at 
30 and 50 passes, respectively, were caused by ap- 
parent thixotropic effects after periods of standing 
in the apparatus. 

Data showing the effects of rate of shear on the 
apparent viscosities at 212° F of X-418 and X-478 
rubbers are given in figure 16. The data for each 
curve shown were obtained with a given sample, 
hence the individual values obtained are affected to 
some extent by the work done previously on the 
sample. These curves deviate slightly from the linear 
when plotted on a log-log scale 

Curves showing the change in viscosity with tem- 
perature for these two rubbers are shown in figure 17 
These tests were run at a rate of shear of 24.3 recip- 
rocal seconds on samples that had been previously 
worked over 100 Over the range covered, 
the relation between log Viscosity and temperature 
for the X-478 rubber was found to be approximately 
linear, while with the X—418 rubber the increase in 
log viscosity with decrease in temperature was pro- 
portionally greater at the lower temperatures 


passes. 


V. Summary 


One of the chief advantages of the worker-con- 
sistometer is that samples of various materials mas 
be worked and their flow characteristics measured in 
the same series of operations. The multihole disks 
give flexibility in covering wide ranges of rates of 
shear and consistency. 

The cam-operated mechanism is suitable for mak- 
ing flow measurements at various constant speeds, 
which give a range of 200 to 100,000 reciprocal sec- 
onds for the nominal rate of shear. The load range 
is from 1 to 200 Ibs. For a broad range of rate of 
shear, this load range corresponds to a VISCOSILY 
range from 1 to 100 poises; and for a relatively low 
rate of shear of 200 reciprocal seconds, the viscosity 
range is from 20 to 4,000 poises. 

The screw-operated mechanism is provided with 
a number of gear trains to give constant speeds cor- 
responding to rates of shear ranging from 1 to 100,000 
reciprocal seconds. Five pairs of calibrated springs 
cover a range of load from 1 to 1,600 Ibs, which 
gives a range of viscosity eight times greater than 
that of the cam mechanism for its range of rate of 
shear (200 to 100,000 sec™'). With the serew- 


M. Mooney, A shearing-disk plastometer for unvuleanized rubber, Ind. Eng 
Chem., Anal. Ed. 6, 147 (1934 











Newtonian materials having thixotropic properties, 


to 


operated apparatus, the range of viscosity deter- 

minations is extended to 600,000 poises at 10 recipro- | The performance of the apparatus was adequate 

cal seconds or 6,000,000 poises at 1 reciprocal second. | cover consistencies of materials ranging from a light 
The test data presented indicate the usefulness of 

the modified worker-consistometer in the study of 


the effect of temperature, rate of shear, and mechani- 
viscosity of non- Wasnineton, July 25, 1950. 


lubricating grease up to and including 100 percent 


of raw rubber. 
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FiGgurReE 1 Major elements of worker-consistometer. 
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Figure 5 Cross section of screw-operated worke 


Fieure 4 Side view of cam-operated worker-consistometer. 
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FIGURE 8 Views of filling and storage 
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Oxide Cathode Base Metal Studies 


Ralph Forman and Glenn F. Rouse 


\ method for comparing the properties of two oxide cathodes, one with a high- 
purity-nickel base and the other with a nickel base containing a small percentage of mag- 


nesium, is discussed 
electrolytic ‘activation and life of an oxide 


Data are presented showing the effect of magnesium on the rate of 
cathode The influence of the cathode and 


anode on pulsed current decay at high current densities is also discussed. 


I. Introduction 


During the construction of an oxide cathode and 
the tube of which it becomes a part, there exists a 
multiplicity of factors that may enter to affect the 
subsequent performance of the cathode. Solutions 
to the problem of obtaining unambiguous data on 
any of these factors would lead to a much better 
understanding of the physical nature of the emissive 
properties of an oxide cathode 

The composition of the metal that serves as a base 
for the oxide coating is widely recognized as a source 
of uncertainty in cathode behaviour. In a majority 
of cases the base is made of nickel that contains 
minute amounts of certain other metals. Studies of 
the role played by these other metals in connection 
with the emission properties of the cathode have been 
made by several investigators [1, 2, 3, 4]' and are 
still a source of considerable current interest. The 
inadequacy of the information regarding the role 
played by the base metal is brought out by the fact 
that a tube manufacturer still feels compelled to ac- 
cept or reject a new melt of cathode nickel, not on 
the basis of available information about the composi- 
tion of the material, but on the basis of test results 
obtained with one or more production tube types. 

The work described in this paper is part of an in- 
vestigation undertaken to obtain additional informa- 
tion about the nickel-base metal. More specifically, 
this paper deals with the study of the influence of 
magnesium in cathode nickel on the thermionic per- 
formance of an oxide cathode. The choice of mag- 
nesium was prompted by its relatively high reducing 
action on an oxide coating [5]. 

A description is given of two distinct phases of the 
work. One is the development of suitable test ap- 
paratus and proper testing technique. The other is 
the analysis of the results obtained from a comparison 
of a cathode on high-purity nickel with one on a base 
of magnesium nickel. 


II. Description of Apparatus and Methods 
1. Test Diode 


A logical method for evaluating the effects that 
result when a small amount of magnesium is added to 
nickel is to compare the performance of two cathodes, 
one on high-purity nickel, the other on nickel con- 
taining magnesium, both of which have been exposed 


Figures in brackets indicate the literature references at the end of this paper 
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equally to all other factors that could influence thei 
behavior. The best practical means of doing this 
would seem to be the use of a properly designed 
double diode. Particular features of the double- 
diode structure that was used are illustrated in fig- 
ure 1. 

A channel made from high-purity nickel is bent 
and welded to form a square frame to which the 
cathode can be fastened (see fig. 1, a). A flat, pure 
tungsten heater, uncoated, is positioned parallel to 
the plane of the frame and equidistant from the two 
faces. Squares of the sheet metals that serve as 
the bases of the two cathodes to be compared (one 
is illustrated in fig. 1, b) are fastened to the cathode 
frame so as to enclose the heater filament com- 
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Fieure 1 Component parts of double diode 


1, Cathode mount; b, cathode: c, sheet anode; d, water-cooled anode 
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This is done, of course, after application of 
The diode structure is com- 
In one case 


ssion coating 
»vy the addition of two anodes 
wiles consist of flat sheets with edges turned to 
rigidity This type of anode 
easy to make, but it has the disadvantage of 
ermal dissipation. Another type of anode is 
ited in figure 1, d. This consists of a body, 
ned from a ‘Grade A” nickel rod, to which a 
evlinder is brazed so as to make possible a 
etal seal. An advantage inherent in this 
is that the part of the anode nearest the 
ithode can be very effectively heated or cooled as 
esired. To accomplish the heating, a coiled tung- 
ien filament, mounted on a press that permits rough 
vacuation of the space within the anode, is placed at 
he position marked A. By controlling the current 
hrough the filament, the temperature of the anode 
wad can be held at any desired value as high as 
j00° C for as long a time as seems necessary. This 
nakes it possible to degas the anode very thoroughly 
to remove any contaminants that can be re- 
noved by heating. The cooling is accomplished by 
placing the filament with a tube through which 
made to flow against the inner surface 
earest the anode face. This method of water 
rooling was used when tests necessitating high anode 
Jissipation were run on the tube. 


see fig. 1, ¢) 
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Figure 2 shows the completed diodes, illustrating 
the two types of anodes It should be that 
the getters are mounted in side tubes 

Some of the advantages to be cained from the use 
of a double diode deserve mention. The tempera- 
tures of the two cathodes being compared are equal, 
every instant during life 
It is realized, of that unequal thermal 
emissivities will result in temperature differences 
However, as the heat radiated by the coated areas is 
a small fraction of the total heat radiated by the 
cathode, it seems reasonable to suppose that, if the 
uncoated areas of the cathode have equal emissivi- 
ties, variations in the emissivities of the coatings 
themselves will effect the coating temperatures only 
slightly During the processing and the subsequent 
life of the tube, the two cathodes are, in general, 
exposed equally to conditions that might affect their 
initial quality and life performance. A_ possible 
difficulty that may arise in the use of a double diode 
is that material present in one cathode may somehow 
migrate to the other. This point must be given 
proper consideration when an interpretation of a 
given set of data is being made 

In the interest of completeness, another fact, not 
dependent on tube design. can be mentioned here 
It is that the two cathode .re subjected to identical 
such as cleaning, degassing, hydrogen 


noted 


or very nearly equal, at 


course, 


processes, 


Figure 2. Doub 





firing, coating application, ete., previous to insertion 
into the finished tube. 

The following comments cover the more important 
details in the treatment of tube components and in 
the processing of a tube 

All nickel parts were chemically cleaned and fired 
at about 700° C in hydrogen that had been passed 
through an electrodryer. Particular care was taken 
to make sure that both cathodes, and both anodes as 
well, were treated exactly alike The two cathodes 
were sprayed simultaneously and in such a way that 
the resulting coatings were essentially identical in all 
respects. The coating weight was approximately 5 
mg/cm’, the coated area 4 em*. The emission 
coating used contained 57.2 +2 percent of barium 
carbonate, 38.8 +2 percent of strontium carbonate 
and 4 +0.5 percent of calcium carbonate. These 
values conform to a specification of the American 
Society for Testing Materials 

The completed tube was evacuated on a system of 
the conventional type containing a three-stage oil- 
diffusion pump with octoil S.A cold trap was 
located between the diffusion pump and the tube. 

The tubes were processed according to a schedule 
that included an oven bake, degassing of the anodes, 
breakdown of the coating by temperature treatment, 
another degassing of the anodes, and an ageing period. 
Details pertinent to any particular experiment will 
be given in the description of the experiment. 

Temperatures were measured with an optical 
pyvrometer. All values quoted are brightness tem- 
peratures 

Data are reported in the literature that clearly 
indicate that the magnitude of the direct-current 
emission obtainable from an oxide cathode differs 
markedly from that of the peak emission observed 
microsecond current pulses are used. For 

Dillinger [6] has shown that a direct- 
current density of 19 amp/cm’*, can be 
drawn from a _ well-activated cathode on a pure 
nickel base (1,150° K), whereas Coomes [7] has 
reported peak emission densities as high as 90 amp 
em’, from a cathode on a pure-nickel base (1,075° K). 
In the latter case l-usec current pulses were used. 

It was believed, however, that a better correlation 
between direct-current and pulsed emission data 
might be found if pulses of relatively long duration 
were used. If this were found to be true, one would 
then be in a position to learn much about direct- 
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current 
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Pulse-measurement circuit, block diagram 

















Ficure 3. 


current emission while taking advantage of the pulse 
technique, the advantage being the ability to measur 
temperature-limited emission without excessive | 
ing of tube parts. With this idea in mind, pulsed 
measurements were made on the experimental dic 
Square wave voltage pulses of 250 usec duration 
used. Setting the repetition rate at 40 pps gave g 
duty evele of 1 percent. A block diagram of th 
circuit used for the pulse measurements is show) in 
figure 3. / is the pulsing circuit. The peak read- 
ing voltmeters, VV, and Vy, are type GR-1S800A 
Meter \, is provided with a capacity-compensated 
voltage divider that increases its range by 10. Mete 
\, is calibrated to read the voltage developed by th 
cathode current across the 100-2 noninductive re- 
The CRO in parallel with meter V, is 4 
DuMont 248 oscillograph. Connections are mad 
directly to the vertical plates of the latter instrv- 
ment, whose sensitivity is known in terms of volts 
inch 


sistor 


III. Discussion of Experimental Results 
1. Cathode Performance Studies 


The magnesium-nickel alloy referred to in thus 
paper is INCO T-17652, for which the manufa 
turer specifies the additives as 0.18 percent of Mg 
and 0.01 percent for both C and Si. 

The high-purity nickel used in this experiment 
INCO 1001 electrolytic nickel. Spectrographic analy- 
sis of this nickel shows approximately 0.01 percent o! 
Cu and Pb. Traces of Mg, Mn, and Si, appear ir 
quantities less than 0.01 percent. 

For ease of description, the magnesium-nickel allo 
and the high-purity electrolytic nickel will be referred 
toas M g-Niand 1001, respectively. Likewise, cath- 
odes made on these materials will be referred to as 
Meg-Ni and 1001 cathodes. 

The first tube for which experimental results ar 
presented was one in which one of the cathodes was 
on 1001 nickel and the other on Mg-Ni. The parts 
of the tube not of 1001 nickel were the tungster 
heater, short lengths of tungsten extending insi« 
the glass press, the emission coatings, and an RCA 
batalum getter in a side arm (see fig. 2,a). Immedi- 
ately following breakdown of the coating, during 
which the maximum temperature reached was 950° 
the cathode temperature was set at 800° C, and 
direct-current potential of 100v. was applied to eac! 
anode. Data showing the variation of current wit! 
time were taken. Curves of these data are show 
in figure 4, A being for the 1001 cathode, B for th: 
Mg-Ni cathode. The rate of increase of curren! 
exhibited by the Mg-Ni cathode is markedly lower 
than that for the 1001 cathode. 

Before any significance can be attributed to this 
difference in behavior, it is necessary to know whethe! 
the curve for the reference cathode, that is, the 100 
cathode, truly represents its behavior. That is to 
say, would the activation curve for the 1001 cathod 
have been quite similar to the one shown if the Mg 
Ni cathode had not been present? Also, what statis 
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ictivation 


base; B, Mg-Nit 


should one expect 


Information bearing on these questions was ob- 


tained from a tube that differed structurally from the 


ne described in just one respect, namely, the Mg-Ni 
Thus both 
Care was taken to 
to which these 


athodes were subjected, identical. After the tube 


had been processed according to a schedule as nearly 
is possible identical to that for the previous tube, 


were taken. Curves of these 
ata are shown in figure 5. 

It will be observed that curve A of figure 4 Is very 
One may infer from 
his that magnesium did not migrate from the one 
athode to the other, or, if it did, that it had no 
fleet. The separation of the curves of figure 5 may 
¢ taken as a measure of the experimental uncertainty 

be expected. It seems, therefore, that the curve 
the 1001 cathode in figure 4 is a typical one, and 

tnay be concluded that the behaviour of the Mg-Ni 

ode does differ significantly from that of the 1001 
ode, 

Che low initial current values in each of these tests | 
an indication that there was no appreciable 
mal activation of the cathodes. 

‘he lower rate of increase of current exhibited by 
\Mg-Ni cathode may be due to the presence of a | 


high resistance layer of some magnesium compound 
at the metal coating interface. D. A. Wright [8] has 
reported finding a MgO interface laver on a Mg-Ni 
cathode. 

It should be pointed out here that the absence of 
thermal activation and the slow rate of electrolytic 
activation have been observed only when the maxi- 
mum temperature used in the conversion of the 
coating was not greater than 950° C. Evidence 
obtained from some tubes indicates that high initial 
currents are obtained from a Mg-Ni cathode for 
which the temperature during breakdown reaches a 
maximum of 1,050° C, and any growth of current 
occurring after breakdown is very rapid. 

Immediately following these tests that were made 
to obtain the data plotted in figures 4 and 5, the tubes 
involved were aged for several hours, and then pulse 
tests were made. It was found that both the 1001 
cathodes and the Mg-Ni cathode could supply peak 
current densities of about 4 amp/cm*. This result 
has been found to be rather generally true for cath- 
odes in other similar tubes. 

The tubes described above were of the type that is 
unsuited for direct-current life testing at normal 
cathode temperature because of the low thermal 
dissipation of the anodes. <A direct-current life test 
has been conducted on a tube having water-cooled 
anodes (see fig. 2, b). One cathode has a base of 
1001 nickel, the other a base of Mg-Ni The anodes 
were made of grade A nickel. After the tube had 
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Figure 5. Activation curves. 


Each cathode has a 1001 nickel base 
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been processed and certain preliminary meas 
ments were made, it was put on a life test sched);); 
like that used by Kern and Lynch [9] of Bell Te‘e- 
phone Laboratories. The test was started by sett 
the cathode temperature at 800° C and adjusting 
anode potentials until the current from each eathod 
was 125 ma (current density 500 ma/cm As 1 
passed, the anode potentials were adjusted so as to 





keep the currents at 125 ma. Curves showing tly 
values of potential plotted as a function of ti 
appear in figures 6, a, and 7, a, the former being 
the 1001 cathode, the latter for the Mg-Ni eathod: 
At intervals the direct-current tests were stopped 
temporarily, and pulse tests were made. The data 
for each of three different times, 361, 878, and 1,142 
hours, respectively, appear in figures 6, b, and 7, } 
The voltage seale is a three-halves power plot, as this 
gives a linear relation between current and voltag 
when space-charge current is being drawn. Points 
that fall on the space-charge line are not plotted i: 
figures 6, b, and 7, b, as they are duplicated for eacl 
curve, 

A study of figure 6, b, reveals that the position of a 
temperature-limited curve drops progressively as 
the hours of life increase, a fact that applies also to 
the point of departure from the space-charge ling 
However, the break from the space-charge lin 
occurs at a point well above the direct-current oper- 
ating point (indicated by a filled cirele) even after 
1,142 hours of life. This behaviour correlates wel 
with that to be deduced from the curve in figure 6, a 
which shows that the voltage required to maintain a 
space-charge current density of 500 ma/em*, had not 
changed during the 1,142 hours of operation. — In 
fact, this cathode was able to maintain a space- 
charge current density of 500 ma/em?, for the duratior 


of the life test period of 3.750 hours 


CURRENT DENSITY, omp/cm* 








ANODE VOLTS x1i072 


le-test data for the 1001 cathode 


A study of figure 7, b, also reveals a progressiv: 
drop in the position of the temperature-limited curvy: 
In this case, however, the break from the space- 
charge line reaches and drops below the direct 
current operating point. Although it is impossib| 
to locate the break exactly, it seems that it is nea! 
the direct-current operating point after about 361 
hours. One will note that in figure 7, a, the upward 
trend in the voltage value starts at about 361 hours 
It can be seen also that the 878-hours curve in figure 
7, b, crosses the 500 ma/em? line at a point corr 
sponding to a voltage of 600 or slightly higher 
This is a fairly good check of the value corresponding 
to 878 hours in figure 7, a 





(3) 


A direct-current density of 500 ma/em?*, is bi 
lieved to be a fairly severe life test condition. As 
the conditions that define end of life for a cathod 
are arbitrary, it may not be correct to say that tly 
Mege-Ni cathode has reached end of life at 1.50% 
hours. However, it seems quite obvious that tli 
1001 cathode should have the longer life. 


CURRENT DENSITY, omp/cm® 
o 
ow 


0 : . A difficulty that arose while the cathode structur 
"ee ey A, ea was being sealed in resulted in a noticeable differen 
ANODE VOLTS x10"2 between the anode-cathode spacings of the two diodes 
Ficure 7. Life-test data for the Mg-Ni cathode As a result, the anode voltage for the Mg-Ni catho«: 
., Direct-current date; b, pulse dat is considerably higher than that for the 1001 catho 














effect, if any, which the higher voltage may 
had is not known. 
ests somewhat similar to the ones presented here 
Kisenstein [10] on cathodes having 
on nickel pure nickel bases He showed 
the emission values obtained for the silicon 
el cathode, using pulses of the order of 1 to 10 
osec, declined more rapidly with time than did 
values for the pure-nickel cathode. He attrib- 
d this behaviour to an interface laver. It seems 
ly that the behaviour of the Mg-Ni cathode can 
wise be attributed to an interface laver 


{ made by 


and 


2. Decay Phenomena 


n the process of pulse testing the tubes to investi- 

te cathode performance, considerable evidence of 
pulse-current decay appeared. This effect has been 
diseussed by many authors, Sproull [11], 58. Wagener 
2), G. R. Feaster [13], D. A. Wright [14], and 
\letson [15]. Particular efforts were made to dis- 
) cover Whether cathodes 1001 and Mg-Ni behave in a 
different manner with respect to the 
decay effect. The conclusion was that there is no 
consistent difference. In general, the appearance of 
a decay effect in one cathode was accompanied by a 
similar effect in the other cathode. Although there 
were sometimes differences in degree, they were not 
consistent enough to justify the formulation of any 
definite conclusions 


: I; 


oticeably 


Some discussions of the decay effect have been 
based on the tacit assumption that the controlling 
factors exist in the cathode only [11, 12]. More 
recently, views have been expressed to the effect 
that anode contamination may also play a major 
role [13, 14, 15] Conclusive evidence as to exact 
location of the cause in any given situation is difficult 
to obtain, because in a simple diode one can presume 
that any treatment given the anode may affect the 
cathode and vice versa. Results will be presented 
now that would seem to indicate that both the cathode 
inode can be involved 


\ tube like that shown in figure 
d while still on the pumps was aged at 130 ma/cm 
12 hours. Pulse measurements were then made 
The results obtained for the 1001 cathode and the 
\ig-Ni cathode were quite similar. Those for the 
cathode are plotted in figure 8, a. It may be 
point out just how the decay effect: was 
ling observed and measured At setting of the 
voltage, readings were taken on the peak- 
ling meters, and the shape of the current pulse 
> noted on the oscilloscope When a decay effect 
ved on the oscilloscope, values of the current at 
eading and trailing edges of the pulse were read 
( oscilloscope These readings are recorded on 
curves open and filled respectively 
re 8, a, reveals that for potentials at and above 
vy there is pulse-current decay accompanied by a 
irture from the space-charge line 
mediately following the test just referred to, 
anode was held at a temperature of 800° C for 
n, the cathode being kept at normal operating 
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Any evolved gas was removed by the 


temperature 
vielded 


data 


pump Repetition of the pulse test 
that are shown in figure 8, b. It will be observed 
that the decay effect and departure from the space- 


charge line now begins at about 900 The change 
brought about by heating the anode is interpreted 
to mean that conditions at the anode do influence 
decay phenomena 

The pulse data in figure 9 were obtained from the 
Mey-Ni cathode of another tube identical to the one 
discussed The open and filled circles have 
the same. significance in the previous figure 
Following this test the cathode was pulse aged at 
300-v peak for 67 hours, care being taken to keep 
the average power into the anode low enough to 
avoid appreciable heating. Another pulse test then 
gave the data that are plotted as half-filled circles 
The evidence of decay has completely 


above 


as 


in figure 9 
disappeared, and the current values satisfy the space 
Obviously the 67-hour period was 


charge equation 
These results 


sufficient, but perhaps not necessary 
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may be considered to be fairly strong evidence that 
the causes of decay were in the cathode itself. They 
are not conclusive, however, because one cannot be 
sure that. bombardment of the anode with 300-yv 
electrons did not indirectly affect the decay. 
Mention has already been made of the fact that 
at certain times during the life test of the water- 
cooled tube pulse measurements were made (see figs. 
6,b and 7,b). A very interesting and significant fact 
in connection with these tests is that there never was 
any evidence whatever of pulse-current decay. It 
was somewhat surprising that such evidence did not 
appear in the later tests, as it was thought that 
changes occurring in the cathodes and at the sur- 
faces of the anodes during several hundred hours of 
operation might bring about conditions that would 
cause decay. It should be pointed out that early in 
the life of this tube, before the life test was started, 
it was possible to induce decay by heating the 


cathodes to a high temperature. This decay was 
removed by heating the anodes and ageing the 


cathode. 


| 


The amount of information concerning  pul.e- 
current decay presented here is hardly enough to 
justify any very conclusive statement regarding «t 
However, one may say that the information suppo: ts 
the view that decay, of the kind described here, is 
associated with conditions that exist during the 
activation period and early life but may not appear 
during the normal operation of a properly processed 
tube. 

The need for an accelerated direct-current life test 
is often pointed out by those who use tubes in large 
quantities, for instance defense agencies. A 
difficulty encountered in setting up an accelerated 
test for emission is the establishment of conditions 
that permit correlation of results with those found 
under actual operating conditions. This being so 
one may be justified in a little speculation concerning 
a possible application of the decay effect to this 
problem. 

Suppose a tube when delivered to the user is sub- 
jected to a pulse test, the pulse being square and of 
relatively long duration. At sufficiently low values 
of voltage the current pulse will doubtless be square, 
and its magnitude will satisfy the space-charge 
equation. As the peak voltage is raised either of two 
things may happen; decay may appear in the current 
pulse, or the current may reach the temperature- 
limited value and deviate from the space-charge line 
without evidence of decay. Performance of the firs’ 
tvpe would indicate inadequate or improper process- 
ing of the tube, not necessarily of the cathode alon 
Performance of the second type would be indicativ: 
of a properly processed tube, and if the break from 
the space-charge line occurs at a value of current 
considerably above the direct-current operating 
value, one might expect long life. Tubes tested in 
this w ay could be life tested later under actual con- 
ditions to establish a correlation. It is realized, of 
course, that the greater structural complexity of the 
usual vacuum tube may make it difficult or impossible 
to apply such a test. One essential condition to be 
met in the test would be that of controlling the duty 
cycle so that the average power input to the tube 
structures outside the cathode would never go above 
that in the specified direct-current operation. 


as 


IV. Summary 


Details of an experimental method for determin- 
ing the relative merits of two oxide cathodes, one on 
a base of high purity nickel, the other on a base of 
nickel containing a small amount of magnesium, 
have been presented. A guiding purpose in the 
design of the apparatus and in the establishment of 
processing and test procedures was to remove, 
actually or effectively, a maximum number of those 
factors that might introduce uncertainty into thy 
interpretation of results. 

The study of conditions associated with the ap 
pearance and disappearance of pulse-current deca) 
leads one to conclude that this effect may prove to 
be a useful tool in evaluating the initial quality of a 
tube. Uncertainty exists as to the location of th: 
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Effect of Changes in Crystalline Structure on the Infrared 


Absorption Spectrum of Cellulose’ 


Florence H. Forziati and John W. Rowen 


The absorption spectra of native, mercerized, and amorphous celluloses are recorded 
for the region from 2 to 16 microns. Cellulose having the X-ray diffraction pattern of the 
native, or type I, crystal lattice is represented in this study by bacterial cellulose, cotton 
fibers, Valonia cell wall, and the residue from the hydrolysis of cotton linters. Cellulose 


having the X-ray diffraction pattern of the mercerized, or type II, crystal lattice is repre- 
sented by regenerated cellulose, the residue from the of mercerized cotton 
linters, and mercerized bacterial cellulose Amorphous cellulose was obtained by milling 
the residue from the hydrolysis of cotton linters in a vibratory ball mill 

Che spectra of the type I celluloses showed sharper and more clearly defined absorption 


methanolysis 


bands than the spectra of the type II celluloses. The most striking differences were found 
in the 7- to S- and 9- to 10-micron regions In the former region of the spectrum of 
type I celluloses, definite absorption maxima of nearly equal intensity were found at approx- 


imately 7.0, 7.3, 7.4, and 7.5 microns. In the spectrum of type II celluloses, the maximum 
at 7.4 microns was lacking, while the absorption at 7.3 microns was considerably more in- 
tense than that at 7.0 and 7.5 microns. Observations in the 7- to 8-micron region of the 
spectrum of the amorphous cellulose were prevented by the use of a suspending medium 
that absorbed in that region In the 9- to 10-micron region, definite absorption maxima were 
found at approximately 9.0, 9.4, and 9.6 microns with the type I celluloses. When cellulose I 
was converted into cellulose II or amorphous cellulose, these maxima merged into a single 
broad band. Conversion of cellulose I into cellulose Il or amorphous cellulose also resulted 
in an increase in the intensity of the absorption at 11.2 microns. Although factors other 
than decrease in type I erystallinity may contribute to these changes in absorption char- 








acteristics, 
cellulose 


A re port of work done under cooperative 
zed by the Research and Marketing Act 
{the Bureau of Agricultural and Industrial Chemistry 


I. Introduction 


Absorption bands associated with crystal structure 
are frequently found in the infrared spectra of erys- 
talline substances {1, 2, 3].2 Accordingly, the possi- 
bility that native cellulose (tvpe | crystal lattice), 
mercerized cellulose (twpe IL erystal lattice), and 
amorphous cellulose might differ in spectral charac- 
teristics was considered during preliminary work on 
the infrared absorption of cellulose. A comparison of 
the spectra of regenerated cellulose films and mineral- 
oil slurries of cotton ground in a Wiley mill to pass a 
200-mesh screen, cotton that had been completely 
converted to the amorphous form by milling in the 
vibratory ball mill [4], and mercerized cotton ground 
in a Wiley mill to pass a 200-mesh screen failed to 
reveal any significant differences. ‘Therefore, it was 
assumed that the crystalline and amorphous forms 
of cellulose were identical with respect: to infrared 
absorption. Later observations on the infrared ab- 
sorption of a variety of celluloses, however, led to 
further consideration of this point and to the conclu- 
sion that cellulose | shows certain spectral character- 
istics not shown by cellulose I] and amorphous cellu- 


lose. These characteristics and the observations 
that led to their recognition are discussed in this 
paper. 

Figures in brackets indicate the literature references at the end of this paper 


the changes may prove useful in evaluating changes in crystalline structure of 


wreement with the United States 
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II. Materials and Methods 


The cellulosic materials used in this study were re- 
venerated cellulose films, bacterial cellulose films 
Empire cotton fibers, the cell wall of the unicellula: 
alga, and the residues from the 
hydrolysis and methanolysis [5] of purified and mer- 
cerized cotton linters,* respectively. The 
erated films were obtained by deacetylation of cellu- 
lose acetate [6]. The films of bacterial cellulose wer 
produced by the growth of Acetobacter rylinum ° on a 
medium, recommended by Aschner [7], that con- 
tained 4 percent of sucrose. The films were har- 
vested at intervals of 2 or 3 days, washed in wate! 
and freed from nitrogenous material by storage in 
5-percent sodium hydroxide for 3 to 4 days. In 
order to prevent oxidation of the cellulose, this treat- 
ment was carried out in an atmosphere of nitrogen 
The purified films were washed, allowed to stand in 
0.5-percent acetic acid for 30 min, and washed until 
neutral to litmus. They were dried on chromium- 
plated or glass plates. then loosened at the edge with 
a razor blade and carefully stripped from the plat: 


Valonia ventricosa, 





regen- 
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rhe authors are indebted to Wanda K 
for the Valonia rentricosa 
‘ Richard E. Reeves of the Southern Regional Research Laboratory of the ! 
Department of Agriculture kindly supplied the residues from the hydrolysis 
methanolysis of purified and mercerized cotton linters 
Cultures of Acetobacter rylinum were given us by G 
Shirk of the National Research Council 


Farr of the Celanese Corp. of An 
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e cotton was purified by the procedure recom- 
ed by the Southern Regional Research Labora- 
of the U.S. Department of Agriculture [8]. A 
if parallel fibers several fibers in thickness was 
ared, cut to the appropriate size, and mounted 
nha liquid-type absorption cell equipped with 
salt windows separated by a silver spacer 50 yu 
and having two screw cap closures. In order 
duce seattering of radiant energy, infrared trans- 
ting liquids such as carbon disulfide, carbon tetra- 
oride, and mineral ot! were placed in the cell when 
sorption measurements were made. Volatile 
ds were removed from the cell with a stream of 
while liquids such as mineral oil were washed out 
ha volatile liquid. By following this procedure, 
vas possible to use the same mat of fibers with a 
es of liquids 
Large were cut from the cell walls 
Vulonia ventricosa that had been preserved for sev- 
il vears in 70-percent alcohol containing 1.4 per- 
cent of formalin. The fragments of wall were washed 
th water, then mounted and dried on the silver 
chloride window of an absorption cell. <A 
window was used to hold the dried fragments in posi- 
tion during the absorption measurements 
Samples of bacterial cellulose and of the 


sections 


second 


residue 


from the hydrolysis of purified cotton fibers were 
mercerized for 10 min. with 25-percent sodium hy- 
droxide at 25° to 28° C, conditions within the range 


shown by Sisson and Saner [9] to give complete mer- 
cerization of cotton 

\bsorption measurements were made on mineral- 
oilslurries of the cellulose powders, that is, the residue 
from the methanolysis of mercerized cotton linters 
and the residue from the hydrolysis of purified cotton 
linters both before and after mercerization and after 
milling in the vibratory ball mill [4]. The slurries 

re prepared by thoroughly mulling 50 mg of sample 
n 0.25 ml of mineral oil between weighted ground 
glass plates. Two absorption cells were used with 
wh slurry. One contained a laver of slurry of 
capillary thickness; the other a laver approximately 
Jy thick 

The spectra obtained with a double beam 
pet rear cing. Baird spectrophotometer {10} 

upped with a sodium choloride prism. The per- 
lormance of the instrument was checked frequently 
against a polystyrene film of known spectral char- 
ristics. No attempt was made to compensate 

lor the absorption of suspending media. Conse- 
quently, in those cases in which a suspending medium 
the regions of the spectrum in which the 
Dotted lines are used 
he figures to represent these regions. 


einaye ie « 


were 


used, 
lium absorbs are obscured. 
} diffraction patterns obtained with a Geiger 
spectrometer with copper radiation, 
54050 A, were to identify the crystal 
e of the samples 


7 ( iter 


used 


III. Results and Discussion 


are shown the infrared absorption 
four samples of cellulose showing the 


figure | 


of 


diffraction pattern of cellulose | (native 
and two samples showing the X-ray diffrac- 
tion pattern of cellulose Il (mercerized cellulose 
The former group consists film of bacterial 
cellulose, lavers of parallel cotton fibers, fragments 
of the cell wall of the unicellular alga, Va/onia 
rentricosa, and the residue from the hydrolysis of 
cotton linters. The samples showing the X-ray 
diffraction pattern of cellulose Il are a film of 
regenerated cellulose and the residue from the meth 
anolysis of mercerized cotton linters. 

A comparison of the spectra of bacterial 
and regenerated (F591) cellulose films shows that 
the absorption bands of bacterial cellulose are, in 
general, sharper and more clearly defined than the 
corresponding bands of regenerated cellulose. The 
difference between the two spectra is most striking 
in the 7- to 8- and 9- to 10-y regions. In the former 
region of the spectrum of bacterial cellulose, definite 
absorption maxima of nearly equal Intensity are 
found at approximately 7.0, 7.3, 7.4, and 7.5 wu. In 
the spectrum of regenerated cellulose, the maximum 
while the absorption at 7.5 yu Is 


X-ray 


cellulose 


of a 


(F621 


at 7.4 uw Is missing, 
considerably more intense than that at 7.0 and 
16 M In the 9- to 10-4 region, definite absorption 


maxima occur in the spectrum of bacterial cellulose 
at approximately 9.0, 9.4, and 9.6 «A single broad 
band showing poorly defined shoulders at essentially 
same wavelengths is found in the spectrum of 
regenerated cellulose. The spectra also differ with 
respect to the intensity of the absorption band at 
11.2 uw. Examination of a large number of spectro- 
grams obtained with the two types of films showed 
this band be intense in the spectrum of 
bacterial cellulose 

The spectrograms obtained with layers of parallel 
cotton fibers (F354, F357) and with fragments of the 
cell wall of Valonia (F309) are inferior to those 
obtained with bacterial cellulose. Nevertheless, 
these spectra show clearly the same characteristics 
between 7.0 and 7.5 w as the spectrum of bacterial 
cellulose. A definite band at approximately 9.0 u 
also evident. Because of excessive scatternig, 
observations at longer wavelengths were of no value. 


these 


less 


to 


Is 


Both hydrolysis and alcoholysis attack the amor- 
phous regions of cellulose, thereby lowering the 
degree of polymerization and increasing the degree 
of crystallinity of the residual cellulose [11,5] 
Such decreases degree of polymerization as are 
indicated by the cuprammonium fluidities of these 
residues would not be expected to result in changes 
in infrared absorption. Accordingly, the spectra of 


the residues obtained by such treatment of cotton 
linters (F444A) and mercerized cotton linters (F257 
may be considered to be essentially the sponte ol 
cellulose I and = cellulose II, respectively. Un- 
fortunately, the absorption of the mineral oil in 
which these residues were suspended for infrared 
absorption measurements obscures the 7- to S-u 
region of the spectrum. In the 9- to 10-u region, 


three absorption maxima are found in the spectrum 
of the residue from the hydrolysis of cotton linters 
The location of these maxima is essentially the same 
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F591, regenerated cellulose film; F621, bacterial cellulose film; F354, F357 


layers of parallel] cotton fibers in carbon disulfide from 2 to 16 microns, in carbon tetrachloride from 6 to 16 microns; F309, fragments of cell wall of Valonia ventrisoca 


F 444A, residue from hydrolysis of cotton linters in mineral-oil slurry; 


as that of the maxima in the spectrum of bacterial 
cellulose. The spectrum of the residue from the 
methanolysis of mercerized linters shows a single 
broad band in this region. The spectra of the 
residues from cotton linters (F444A) and mercerized 
cotton linters (F257) resemble, respectively, those of 
bacterial and regenerated cellulose with respect to 
the intensity of the absorption band at 11.2y. 
Cellulose I is converted by mercerization to 


cellulose II, and by vibratory ball milling to amor- | 


phous cellulose. Hence, the effect of these treat- 
ments on the infrared absorption of bacterial cellu- 
lose and the residue from the hydrolysis of cotton 
linters is of interest. In figure 2, spectra obtained 
from mercerized bacterial cellulose and from the 
mercerized residue from the hydrolysis of cotton 
linters are presented. Also included in this figure 
are spectra that show the effect of milling in the 
vibratory ball mill for 15 min. and 1 hour on the 
residue from the hydrolysis of cotton linters. As 
expected, the spectrum of the mercerized bacterial 
cellulose (F456) closely resembles that of regenerated 


40 


F257, residue from methanolysis of mercerized cotton linters in mineral-oil slurry 


cellulose (fig. 1, F591). The spectra of the mineral- 
oil slurries of the mercerized (F462) and vibratory 
ball-milled samples (F536, F559) of the residue 
from the hydroloysis of cotton linters, like those of 
mercerized bacterial cellulose, regenerated cellulose, 
and the residue from the methanolysis of mercer- 
ized linters show a single broad band in the 9- to 10-4 
region. 

The foregoing observations suggest that resolution 
into absorption maxima at approximately 9.0, 9.4, 
and 9.64 is peculiar to celluloses that show the 
X-ray diffraction pattern of cellulose 1. Conversion 
of cellulose I into cellulose II or into amorphous 
cellulose spreads these maxima into a region of gen- 
eral absorption with poorly defined shoulders at these 
wavelengths.° The presence of a single broad band 
in the 9- to 10-u region, however, is not proof of the 
absence of cellulose 1. Thus, cotton that has been 
milled in the Wiley mill to pass a 200-mesh screen 

* During the course of this work, it was learned that Stanley C. Burket |12] 5 
Richard M. Badger of the California Institute of Technology had observed t! 


same maxima in the spectrum of cellulose from the marine anima), Ciena int 
nalis. They also found that mercerization resulted in decreased resolution 
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gives the X-ray diffraction pattern of cellulose I but 
fails to show the three distinct absorption maxima 
in the 9- to 10-u region. Whether the failure to secure 
resolution is attributable to the amorphous material 
in the milled sample [13] or to some other factor is 
known. The fact that better resolution is ob- 
tained with bacterial cellulose films than with min- 
eral-oil slurries of the residue of the hvdrolysis of 
cotton linters indicates that the sharpness of the 
bands in the 9- to 10-u region is not determined 
solely by the proportion of cellulose I present. 
Nevertheless, the possibility of using changes in the 
9- to 10-4 region of the absorption spectrum of 
cellulose as a measure of change in the degree of 
crystallinity (type I) merits further investigation. 
Because of the absorption of the mineral oil used in 
the infrared absorption measurements on cellulose 
powders, observations in the 7- to 8-4 region were 
limited to films and cotton fibers. These observa- 
indicate, however, that mercerization and 
egeneration of cellulose result in a change in the 
tensity of the absorption at approximately 7.3 and 
7.4u, as well as in a general broadening of absorption 


not 


Some of the absorption bands observed with native 
lose have not been previously reported in the 
ature. Accordingly, a complete list of the bands 
d in infrared absorption spectra of bacterial 

cellulose films is shown in table 1. Suggestions as to 

dentity of the bands are also included in the 

These suggestions are based on published 
pretations of spectra of other compounds. 


TaBLe 1. Infrared absorption spectrum of bacterial cellulose 
Seah ame Intensity rentative assignments Related reference 
m 

3 01 3,320 | Strong OH stretching [15] 

1 48 2 870 Medium CH stretching 115, 16] 

6.09 1. 640 do HOH bending ad [17] 

sorbed water 
6.95 1. 440 lo HCH deformation [15, 18, 19, 20 
0) 1, 370 d 

7.45 1, 340 do CH bending [15, 16 

7. 55 1,325 do 

7.78 1, 285 Weak 

ag 1. 250 do CH rocking in ring {21 

8. 10 1, 235 do 

8. 30 1. 205 do CO stretching {15] 

s 1,170 Strong Ring frequency 16, 14, 21, 22, 23 
x9 1, 115 de coc [19] 

9. 38 1, 065 d OH bending 115, 18, 19, 2 

9. 60 1,04 de C-OH stretchin {19 

0. 80 1.040 W eak 

9 1, 00 do 

11.20 SU4 d 


IV. Conclusions 


The results indicate that conversion of cellulose I 
into cellulose II or amorphous cellulose produces a 
decrease in the sharpness and definition of the in- 
frared absorption spectcum, as well as changes in the 
intensity of absorption at certain wavelengths. Cel- 
lulose I is characterized by absorption maxima of 
nearly equal intensity at approximately 7.0, 7.3, 7.4, 
and 7.5 uw, whereas cellulose II shows maxima at 
7.0, 7.3, and 7.5 uw. Of these maxima, that at 7.5 u 
is considerably more intense than those at 7.0 and 
7.5 wu. No observations were made in this region of 


113375—50 i 4l 


the spectrum of amorphous cellulose. With cellu- 
lose I, clearly defined absorption maxima are found 
at approximately 9.0, 9.4, and 9.6 ». On conver- 
sion to cellulose Il, these maxima merge into a 
single broad band with poorly defined shoulders at 
these wavelengths. Conversion to cellulose Il or 


amorphous cellulose also results in an increase in the 
intensity of the absorption at 11.2 4 


The authors express their appreciation to Earle 
K. Plyler for many helpful suggestions, Howard E. 
Swanson for the X-ray diffraction measurements, 
Sanford B. Newman and Emil Borysko for prepara- 
tion of the bacterial cellulose, and Walter K. Stone 
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Analysis of Mixtures of Olefin Hydrocarbons Produced by 
Codimerization of Butenes'’ 


Augustus R. Glasgow, Jr. 


his paper describes the analysis, 
mixtures of olefin hydrocarbons 
reactions are called ‘‘codimers”’ 
made by 


of such 
were 


with 
These were produced by the combination to form dimers 
of mixed butenes in a tubular type reactor with phosphoric acid as catalyst 
and 


respect: to 


the process 
utilizing distillations performed at high efficiency 


individual components, of three 
Products 
The analyses 
reflux ratio 


codimerization”’ 
with high 


with accurately measured values of the boiling point and refractive index of the fractions 


of distillate 
In the 
mixt 


case of one codimer 


I. Introduction 


\s part of its war research program, the American 
Petroleum Institute Research Project 6 at the Na- 
tional Bureau of Standards was called upon, by the 
Fechnical Advisory Committee of the Petroleum 
Industry War Council and the National Advisory 
Committee for Aeronautics, to determine the com- 
ponents of mixtures of hydrocarbons produced by 
commercial processes important in the preparation 
aviation fuels. Included in_ the 

ogram were three samples of mixtures of olefin 
hydrocarbons produced by the codimerization of 
butenes and, in addition, a sample of the product 
obtained by hydrogenation of one of the 
mixtures Since, in- the petroleum industry, a 
mixture of olefins such as that produced by the 
ombination of butenes is called ‘ta codimer’’, and 

product formed by hydrogenation of the codimer 
s called “hydrocodimer’’, these terms will be 
the text 

The analysis of various mixtures of alkylates* and 
hvdrocodimers, including the hydrocodimer referred 

n this papel by analytical distillation and the 

asurement of physical properties of the distillates 

. previously deseribed This report 
method and results of the analy ses of the mixtures 

monoolefin hydrocarbons that comprise the 
The analysis of a mixture of olefins is 

complicated than the analysis of a mixture 
the corresponding paraffin hydrocarbons because 
a) the many more 
ers, (b) the lack of accurate data on the physi- 
properties of all the pure components, and (c 
vreater susceptibility of the olefins to oxidation 


II. Material Analyzed 


he material analy zed was the product of the codi- 
ation of material, which was mainly a mixture 


of high-octane 


above 


used 


vives the 


mers 


presence of close-boiling 


and with separation by adsorption of the paraffin portion of selected distillate 
use was made of the previous analysis of the corresponding 
ire of paraffins produced by hydrogenation of the codimer 


of butenes, in a tubular type reactor with phosphoric 
acid as the catalyst (the process of the Universal 
Oil Products Co.) The samples are listed accord- 
ing to the American Petroleum Institute Research 
Project 6 numbers. Samples numbers 38 and 39 
were prepared at the Bayway Laboratories of the 
Standard Oil Co. of New Jersey as Bayway test runs 
11 and 12, respectively. Codimer sample number 
15 and hydrocodimer sample number 44 were pre- 
pared at the Heysham Refinery of Trimpel Limited, 
London, Eng., and had the designation HSO 1685 
and HSO 1684, respectively The analysis of hydro- 
codimer sample number 44, which was the hydro- 
Parte |. Ope 


ating conditions and product inspections 


samples anatiyred 





genated product of codimer sample number 45, was 
previously reported (see footnote 3). 

The pertinent operating conditions and inspection 
data on the samples, submitted by the above com- 
panies, are given in table 1. The analyses of the 
olefins in the fresh feed to the reactor and in the de- 
butanized overhead from the reactor, expressed as 
the disappearance of olefin gas in table 1, give a 
measure of the components reacting in the process. 
The large decrease in the butene content shows that 
the predominant reaction is between the butenes. 


III. Method of Analysis and Results 


The analyses were made by utilizing distillations 
performed at high efficiency with high reflux ratio, 
with accurately measured values of the boiling point 
(to +0.01 deg C) and of the refractive index (to 
+ 0.0001) of the fractions of distillate, and with 


TABLE 2 


Sample 


auesher Material 


Distilling 
column 
number 


plates at 


(approximat 


Codimer, TAC, Bayway 11 
Codimer, TAC, Bayway 12 
Cedimer, HSO-1685 


Changes in composition of the distillate from the 
analytical distillations caused by oxidation, polymer- 
ization, ete. of the component olefins, were mini- 
mized by storage of the distillate at 8° to 10° C in 
amber glass bottles with screw caps. To determine 
the magnitude of any changes in composition, meas- 
urements of the refractive indices of the distillate 
from codimer sample 39 were made as follows: 
(a) Within 1 to 2 days after collection from the dis- 
tilling column with storage at 8° to 10° C; (b) after 6 
weeks of storage at 8° to 10° C; and (c) after 6 weeks of 
additional storage at room temperature (near 25°C). 
These data are given in figure 1. The results given 
in the figure as the two solid lines of refractive index, 
for conditions (a) and (b) above, show that under 
refrigerated conditions the refractive index remained 
substantially constant. However, after 6 weeks of 
storage at room temperature, changes in the refrac- 
tive imdex of the material occurred, as shown by 
the broken line. These data on refractive index 
show that in the analysis of material of this kind 
it is necessary to refrigerate the fractions of distil- 
late and to determine the properties 
without delay. 

Reduction of the data from a large-scale plot of the 
boiling point and the refractive index of the distillate 
as a function of the volume of the distillate, to obtain 
the amounts of the individual components, was made 
in the same manner as that previously described for 
the analysis of the paraffin mixtures of alkylates and 


hecessary 


| 


total reflux 


ely 


44 


separation by adsorption of the paraffin portion | f 
selected distillate. 

In preparation for the analytical distillations, t 
bulk of the peroxides was removed from the samp]. s 
of codimer by simple distillation in vacuum, with tly 
distilling temperature near or slightly below room 
temperature and the temperature of the receive: 
about —78° C. Hydroquinone was added to t! 
pot and the receiver to serve as an antioxidant 
The bulk of the peroxides was contained in the un- 
distilled portions, which were withheld from further 
analysis. These residues, expressed as percentage 
of the volume distilled, were for samples 38, 39, and 
45 about 10, 10, and 1 percent, respectively. The 
distillates, after the addition of 25 g of hydroquinone 
to each, were subjected to analytical distillations, 
the details of which are given in table 2. Complete 
details regarding the distilling operations are given 
in a paper by Willingham and Rossini.‘ 


‘C. B, Willingham and F. D. Rossini, J. Research N BS 37, 15 (1946) RP17%4 


Information on the analytical distillation of three codimers 


Distillation 


Number of equiv 
alent theoretical 


Volume of Results Analysi 
each plotted in | given in 
fraction figure table 


Rate of 
collection 
of distillate 


Reflux ratio 
(approxi 
mately 


Volume of 


charge 


hvdrocodimers (see footnote 3). In calculations for 
sample number 45, use was made of the results of 
the analysis of the corresponding hy drocodimer undet 
the assumption that the process of hydrogenating 
the codimer to form the hvdrocodimer did not alter 
the carbon skeleton of any significant number of 
molecules. This assumption appears to be substan- 
tially correct since it served to produce a logical 
correlation among the amounts of the components of 
the respective mixtures, as shown in table 4. Fur- 
thermore, to simplify the analysis of codimer sample 
number 45, the presence of 2,2-dimethylhexane 
(0.9 percent) in the hydrocodimer and corresponding 
olefins in the codimer was ignored. 

Upon inspection of the data, it was evident that 
paraffin impurities were present in the codimers 
For example, in codimer number 38, as shown in 
figure 2, a flat occurs in the boiling point curve at 
the boiling point of n-pentane. Moreover, the cor- 
responding index of refraction is lower (Np at 25 °C 
1.3623) than would be expected for any olefins boiling 
in this region, although higher than the refractiv: 
index of pure n-pentane (1.3547). Therefore, re 
course was made to filtration through silica gel, to 
separate any n-pentane from the olefins. After 
passage of a 1-ml portion of this distillate through 
silica gel ° the first few drops of filtrate had refractive 
indices increasing successively from 1.3531 to 1.3546 


R. Glasgow, Jr., C. B. Willingham, and F. D. Rossini, Ind. Eng. Chet 


41, 2202 (1949 
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of the refractive index and boiling point for the 
two properties lhe relative amounts by 


values 


ertainties in the 
er ng these 


ind circles d 
The adsorbed portion, after displacement by ethanol, | 
was found to have the refractive index 1.3757, indi- 
cating that the material boiling in the flat was largely 
n-pentane with some isopentane (Np at 25° C, 1.3509) 
and some olefin. 
Similarly, the refractive index curve for codimer 
mber 45 indicated the presence of paraffin hydro- 
arbons. Therefore, 10-ml portions of fractions of 
distillate that appeared to contain paraffin 
lrocarbons (A, B, C, and D in fig. 3) were filtered 
ugh a column of silica gel (28 to 200 mesh, about 
m in diameter and 50 cm in length). One milli- 
of the first portion of filtrate was then passed 
ugh a small column footnote 5). The 
itive indices of the paraffin portions and the 


(see 


IN 


ve 


2500 


3000 
ML 
stillation of codimer sample number 39 

rents are indicated by the limits of the arrowed line drawn through 


pure comyx 
din the upper portion of the figure 


“ume of the various components are indicate 


components of the mixtures, as deduced from the 


refractive index 
1.3798, 2,2- and 
2,3-dimethylpen- 


boiling point together with the 
(Ny at 25°C), were as follows: A, 
2,4-dimethylpentane; B, 1.3896, 
tane; C, 1.3890, 2,2,4-trimethylpentane; and D, 
1.4025, 2,3,4-trimethylpentane plus some 2,3,3-tri- 
methylpentane. The amount of 2,2,4-trimethyl- 
pentane in sample number 45 and its apparent 
absence in samples 38 and 39, as indicated by the 
boiling point and refractive index curves, are shown 
in figure 4. The circles and crosses in the figure 
refer to the boiling point and refractive index, 
respectively, of the pure compounds. The volume 
of each of the components in the distillates, based on 
the reduction of the boiling point and refractive 








PERCENTAGE 








4435 


3NINON 





» IN3LN3d+ 
TAHA SARL OF 2, 


INI 1N3d-2 


VAHL IL -p''2 


S3IN3L90 


INALN 3d 2 


HL BAIO-6°2_ 


Md WOML-P OZ 


} 


S3N31d34 ONY SONVid3H . 


> INIX ZH ONY 


> INV x 3H 


INILNG-Z 
IN 31 N3d-2-5 


INFLNG 


AHL IN-2 a =~ 


ONY -SNVuL | 
INVINGd- 


is 
VAHLIN-2 ONY JNDINGS 


BAYWAY 


TAC CODIMER, 


38 


SAMPLE NO 


» « 
a 
a 
x 
ra 
4 
a 
<4; 
a 
z 


ree 


BF 


INV LN 3d0S 


4500 


3500 


2000 


1500 


O0r 








numbe 


sam ple 


mie 


) 
t 


oO 


on 


sf [ 


analytical 





= vs 














BY VOLUME 


PERCENTAGE 





100 


90 


70 


50 


40 


30 


20 











) a A i) ~ 
¢ ay s +: ? $ m © 
222 
~ 

t o J ee 
- Or 
SNId Sve < 

SNI43 c 

© 

t ° 7 

SNAsveNd ° 4 

SNIZZW | - 14 z 

SNAsVeNd ONY SNIZ1 St o 

ow 
IN3X3H-2-UHL3NO-¢'2 4 5 - ™ 
NW X3H AHL Nd-'e TF ‘ a6 Q 
x “ ¢ 
VO 
a <a 
q w 
w « 
W” 2 
w @ 
INILN3d-2 at x 
VAHL SLL ES 4 ° 
4 4 - oO 
ae ” 
3 
<a r 
= 
t ! 

INV LN3d UHL IL-2 , =y S 
S3N3190-1 1 i o 4 
53N3190-| 2 ©} 

INFIN3d : ow 
2-H NLL vE-(H-9) - a > — ¢ 
IN3X3H -| -WHLawo-re-t , —w—< -< oO Oo 
swwwi007 fhe . nz i 
IN3LN3d-1-TAHL INL E'S’ + xr S 
3N3X3H-£- VAHL IG-b'2-0h°2)" 4 Ww “ 
ee 
INBLN3d-4- AHL IL -b'S'Z - om fe ~- us 
2a 
aw 
t oO 4 
INBIN3d-2-TAHLIWL-b'b'Z F o ° 
INBIN3d+ AHL ms -'e'e | _—_*—- > 
INSIN3d4- AHL IL -b'b't > ee -o 
a TT 
°o 
° oO 
< i 
c a 
(=) @ . Qa 
. = 
oO 
INF LN3d-+- AHL IL-2 * us 
a ° 
> ° 
a 6 
a ! 7 ; 
INVINSd VAHL INL 22 | | . 
Ff oO 
INZLIN3d-Z- AHL BNNO-€'2 _—* 2 
INVINSd AHL INO-2'2 3} i 
INFLNG+- AHL 3-2- UAH wee - =): * : a 
IN3BLN3d-+- AHL NO-€°2_ ONY SNBINGd-Z- TAH JN 2 si <- ’ . 
W3aie-2 INVIN3d TAHLNIO-b'¢ ONV-2'29) ) ~~ 1 — 1. ~» - 
5» ct ANGUNGd 2-H NI -t-(h°9) TF +0 + + S N S3N3.1N3d 
‘WiLared €< WNGIN3Gd + UHL NO -b'p ts = x INVINGd ONY aaete i. 
: “ : > & z z 4 a 
v x ’ N : vw 


OLUME 


Vv 


n imbe y 


‘ 


samp 


of codimer 


distillation 


analytica 


3 Results of the 


URE 


Fi 


d by the 


ipounds are indicate 
ym pone 


cou 


uw 


n 


ind boil 
The relative 


eX 


portion of the 


upper 


smounts 


two propertics 


1 LD we 


vy adsorption 


it 


re process 





47 














os ae BP at 725 ww 


TEMPERAT 





r2 = 2 zy zy 
we = = 2 " 
2: 22 bf 2° bf t Fo } 
Fs + t> FS 1.45 
: tt * 
7? ? 
+ aa 
; ; —- 
> 41.43 & 
DIMER SODIMER SODIMER ms 
AME 10. 38 APL 40.39 SAMPLE NO. 45 s 
py a2 & 
* * . z 
H , 
es} 4 
X an 
— = = = = = * . 
= 
, 
80r 4 
| 
75} i 39 
| 
| 
Se ee = — | rt se 
aX 4 25 4 
“ vv 
Figure 4. Comparison of distillate from codimers in the 


region of the boiling point of 2,2,4-trimethylpentane. 


index data, is indicated in the upper portion of the 
figure. The broken lines accompanying the refrac- 
tive index curves indicate the division into compo- 
nents, as deduced from the above properties. 

The results for codimer samples 38, 39, and 45 
are given in figures 1, 2, and 3 and in table 3. The 
material boiling below pentanes and pentenes, which 
constituted 7.5 and 5.0 percent of the total volume 
distilled for samples 38 and 39, respectively, had the 
following composition expressed as percentage of the 
total volume distilled: 

Sample number 38, isopentane plus lower-boiling 
material, 1.7 +0.4; l-pentene plus 2-methyl-1-butene, 
0.6 +0.4; n-pentane, 3.1 +0.8; cis-plus trans-2- 
pentene, 0.8 +0.3; 2-methyl-2-butene, 1.3 +0.3. 

Sample number 39, isopentane plus lower-boiling 
hydrocarbons, 4.0 +0.3; pentenes plus n-pentane, 
1.0 +0.5. 

From the analyses given in table 3, the relative 
amounts by volume of the trimethylpentenes in the 
codimers that would yield the same trimethylpentane 
on hydrogenation without alteration in the carbon 
skeleton were as follows: 

(a) 2,4,4-trimethyl-1l-pentene, 2 to 4 times that of 
2,4,4-trimethyl-2-pentene; (b) 2,3,4-trimethyl-2-pen- 
tene, about two times that of 2,3,4-trimethyl-1- 
pentene; (c) 3,4,4-trimethyl-2-pentene, two to three 
times that of 3,3,4-trimethyl-l-pentene; and (d) 
2,3,3-trimethyl-l-pentene, slightly greater than 3,3,4- 
trimethyl-1-pentene. 

Figure 5 shows the distribution of olefins in the C, 
portion of hydrocodimer, sample number 44. 

The open and solid circles give the boiling points 
of the pure paraffins and olefins, respectively, and 
the crosses refer to the values of the refractive indices 
of the pure compounds. The refractive index of the 
original distillate containing olefins is given by the 
broken line and that of the paraffin portion, by the 
solid line. The difference in the refractive index of 





the paraffin portion and that of the original distill; 
is shown in the lower portion of the figure. Tie 
areas obtained by plotting the increase in refractive 
index due to olefins as a function of the volume of 
distillate, together with the corresponding plot of 
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TABLE 3. 


Codimer 
: 
Boiling point Sample No 
at | atm 


4,4-Dimethy! 
2,3-Dimethy1-2-butene 


pentene 








Components of material boiling above pentanes 


Percentage by volume 


and pentenes 


| 
Codimer 
ample No. 45 


Codimer 
Sample No 


tS 


above Cs componer 


4.4-Dimethyl-2 pentene * 76 0.7 +0 4+0.4 9 20 +0 
Hexenes 41 to 70 | 0.7 0.4 
j Hexenes »+hexanes 49 to 0.7 0 
' 2,2-Dimethylpentane 79. 2 ) : 
: 2,4-Dimethylpentane sO) p 0.4 40.2 
2, ‘ Dimethyl-2-pentene 82.4 - ss 
2,3-Dimethy]-1-pentene aS es _ if ¢ u 
3. Methyl-2-ethy!-butene 89 £0 ( 0 +0 O+0.5 | 8.9 0 
2,3-Dimethylpentane 80.8 | O08 0.4 | 
Heptenes 6 to 0 +0 | 20 40.7 | 
2,3- Dimethy]-2-pentene v7 as & .3 i 1.¢ ? 14 } 
s ; 2,2,4-Trimethylpentane 09.2 L o- ) , 7 2 ~~. 
. 2,4,4-Trimethyl-1-pentene o1.4 W6241.3 J ~ 19+0 1.9 0.4 0.4 41.0 | : 
t f 
F ( ctenes 02 to 104 yu +0.9 | 
‘ +,4,4-Trimethyl-!-pentene 04 — ‘ . 226 0 ‘ 
: }.4-Trimethyl-!-pentene 10 | ‘ | 2 29 2 40.6 9 2 
= ; 2,4,4-Trimethyl]-2-pentene 04.9 7.8 +1. 2 0.9 ri - 
° 4 .4-Trimethyl-!-pentene OS S +3.0 0 +! 7.8 +20 
: 4 Dienethyho-hosen . ON 7 40.8 
i 2,3,3-Trimethyl-1-pentene 108. 4 4.9 423 9409 
| Octenes OS to 112 0.5 +0.3 a8 . 
P , 2h. f 2.0 26 2 yo 4 } 
Octanes 09 to 110 0 0.2 - 
3.4-Dimethyl-1-hexene 112 { { 1 OR ( 
£ 3.4.4-Trimethyl-2-pentene * 112 .7 ¢ 7.2 6 4 es 
z Octenes 110 to 112 20 0 
4 Octenes 1I2to 115 } . 7.2 +2.0 
| Octenes to1ie jf °*9 = | P 2.4 | - . | : 
; 2,3,4-Trimethyl-2-pentene ! 16.3 16. ¢ 0.6 | . re nm 0 r 27.6 +2.0 1a 0 - 21.1 ! 
J 2,3,4-Trimethylpentane 11 | 22 0 | 
é 3. 4-Dimethylhexane 17 } } 0 +0. OF 
7 2,3-Dimethyl-2-hexene * 22 9 +0.4 6.4 40.4 6 +0.7 7 +40 
Octenes @ 117 to 12 YW 0.4 | 6.4 ().¢ | | 
Nonenes 24 to 13 +f 0 
Nonenes 24 to 149 & +0 
Olefins+ paraffins 23 to 14! 29 40.7 
Nonenes+ higher olefins 136 6 +0) ) 1 +0 24.8 +0 2407 
Olefins+ paraffins 149 0 +0.4 { 
Olefins+paraflins 141 to l¢ 940.7 | 
Olefins+ paraffins > 165 j 4 +( 
3 rotal 100.0 100.0 on 
* (ts and trans isomers 
May contain some 4,4-dimethyl-!-pentene (a heptene 
May include other octenes 
May include-some nonenes 
May include some octanes 
Includes 2,2,3-trimethylpentane, 2,5-dimethylhexane, and 2,4-dimethylhexane 
« May include some 3-methy1-2-isopropyl-!-butene near 104° C 
hese octenes may include 3,3-dimethyl-2-ethyl-1-butene, 2,4-dimethy]-2-hexene, 4-methy]l-2-ethyl-1-pentene, 2,4-dimeth hexene, 2.5-dimethvl-!-hexene 
~dimethyl-2-hexene, 2,5-dimethyl-2-hexene, and 3-methyl-2-ethyl-!-pentene : 
Includes small amount of 2,3,3-trimethylpentane and 2,3-dimethylhexane 
May include some +t 4-dimet hy]-2-hexene 
* May include some 4-dimethy1]-3-hexene 
boiling point versus volume, show the distribution of | used to calculate the shapes of the broken curves 
the olefins in the distillate of the Cs paraffins. The lowering in refractive index due to paraffin was 
Figure 6 shows the distribution of major C, par- | then calculated as the difference between the broken 
affins in the Cy portion of codimer sample number 45 | curve and the observed refractive index. The areas 
ind the method used in deducing this distribution. | A, B, and C obtained by plotting the decrease in 
In figure 6 the open and solid circles give the boiling | refractive index due paraffin as a function of the 
1) nts of the pure paraffins and olefins, respectively, volume of distillate, together with the corresponding 
id the crosses refer to the values of the refractive | plot of boiling point — volume, show the dis- 
lices of the pure compounds. The uncertainty in | tribution of the major Cy paraffins in the distill ate 
h of these values is indicated by the limits of the | to be as follows: A, > 2 4-trime ‘thylpentane; B, 
rowed line drawn through the property. The | trimethylpentane; C, 2,3 4-trimethylpentane 
‘ken curves accompanying the observed refractive | 2,3,3-trimethylpentane. 
lices (solid lines) are the curves calculated for the Figures 5 and 6 also show that no significant 
n ae 
ration that would have been obtained in the | number of molecules underwent alteration in their 





tilling column if no paraffins had been present in 
material. 
condition 


The b locked broken lines were for the 
with perfect separation, and were 


ne 


carbon skeletons during hydrogenation of the 
codimer to form the hydrocodimer. That is, the 
predominant C, olefins present as impurity in the 


49 














hydrocodimer have the same carbon skeleton as the 
predominant paraffins, and similarly the predominant 
C, paraffins present as impurity in the codimer have 
the same carbon skeleton as the predominant olefins. 

Table 4 is the correlation of the olefins found in 
codimer sample number 45 with the paraffins of the 
same carbon skeleton in the corresponding hydro- 
codimer (sample number 44) for the material boiling 
above pentanes and pentenes. The material boiling 
below pentanes and pentenes was 0.8 and 0.2 percent 


of the total 


(see footnote 3). 


volume distilled for the codimer 
hvdrocodimer, respectively. 
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Fred: 


Tapie 4 Correlation of the olefins in the codimer with paraffins of the same carbon skeleton in the corre sponding h ydrocot 


for the material boiling above 


Cy and higher 


livd wir 
: N 4 

( npone Volume percent 
Hexa 0.840 
22-1) t} 1M 
2 Din I " > 440 a 
2.4-Dimethylpe ul ~ v 
Olefir te wt ( “1 ‘ | 
»3-D tl went 7.4 i |. 
- ir 
Olef i] nw f 1 a 
Olef wer an | vf +t lon o 9 
2.2.4-Trimethyly 27.4+1.0 {* 
2,2-Dimethylhex 0. +4 
2.5-Dimethythex » Ons 
2. 4-Dimethy lhe t { oo 
2 Ir hyip 7.9 be 
Olef vt ” ¢ 240.1 ] 
2.3,4-1 ‘ | i 28. 3+2.0 low g >, 
Olefins 114° t wie 0.640.3 ] 
2 t ul 10 
2 3-Dimethylhex +1.0 
i-Dimett hexa 11+. ¢ 
Isononanes 122° to 144° ¢ 2. 640.5 | 5 sig 8 
Olefins 120° to 144° ¢ u.2+.1 f 
Isodecanes 145° to If ( 1.3+.6 
Olefins 144° to 178° ¢ 0.4+.2 

> 2+). f 
lsoundecane t ae ¢ 204.9 
Paraffins 17s” 7.32.6 
Olefins 178° 0.34.1 
Total 100.0 


* Cis and franea ner 
Mav include smaller amount of 2,5-dimethylhexane and 2,4 dimethylhexane 
The predominant olefin in this boiling range 


The five olefin isomers of the same carbon skeleton as 2,2-dimethythexane 


},4,4-trimeth yl-2-pentene, which 
The predominant olefin in this boiling range Is 2,5 4-trimethyl-2-pentene, which is grouped with 2,3,4-trimethy Ipentane 
normally boiling in the region 101 


Components 


dimethyl-l-pentenc 
bimethyl-2-pentene 
Jimethyipentane 
dimethyl pentane 
Jimeth yl-2-pentene 


2,3-Dimethyl-1-pentene 
- Methyl-2-ethyl 
2,3-Dimethyl pentane 

$- Dimeth y1-2-pentene 


rimethyvipentane 


l 
rrimethyl-l-pentene 
rrimethyl-2-pentene 


lrrimethyl-l-penter 


rrimethy!lpentame 


rrimethyl-2-penter 
me lt to 1124°C 


lrimethyl-l-penter 


rrimethylpentane 


}.4-Trimeth yl-2-pentene 


3. 4-Trimethyl-l-pentene 
rrimethyl-l-penter 


}- Dimethyl-2-hexene 


3,4-Dimethyl-1-hexene 
$.4-Dimethyvlhexane 


Olefins 123° to 141° C 
Paraffins 123° to 141° C 


Olefins 141° to 165° C 
*araffins 141° to 165° C 


Olefins 1s? C 
‘araffins 15° C 


> 


is grouped with 2,5,3 


somers made impracticable ar lirect assignment of the 0.9 percent of this isomer in the codimer 


W asninaton, October 2, 1950 
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1-tuteme 


pe ntanes and pe ntenes 


Codimer 
Sample 


trimethyl pentane 


No 


2040 


to 107 





percent 


7+0 

6+0 

t 0 
) 0 
7+0.8 
O+1.4 
++1.8 
1+). S 
» s 


100.0 


to large concentrations of 
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Surface-Area Determination by Adsorption of Nitrogen 


from Nitrogen-Helium Mixtures’ 
William V. Loebenstein’® and Victor R. Deitz 


A new apparatus and technique, in which vacuum pumps and associated equipment are 
not required, have been developed for the determination of the adsorption of nitrogen from 


mixtures with helium for use in the determination of surface area 
The adsorption of helium is negligible at the temperatures 
The gas mixture is cycled repeatedly over the adsorbent until no further change 
Temperature equilibration is attained rapidly 


initially in a stream of helium 
concerned 


in pressure is observed 


The sample is heated 


It is shown that 


thermomolecular pressure differences are negligible under the conditions of the experiment 
Surface-area determinations by the simplified technique are in good agreement with results 


obtained by 


various residual helium pressures 


I. Introduction 


he volumetric method for the determination of 
s adsorption is still the most popular method in 
ise today. Since the investigations of Titoff [1] 

and of Homfray [2] in 1910, the apparatus and the 
technique employed for gas adsorption have remained 
essentially unchanged. Certain refinements have 
been introduced, such as the substitution of mercury 
traps for stopeocks by Coolidge in 1924 [3] and also 
by Goldmann and Polanyi in 1928 [4], but such re- 
finements did not change the basic technique and 
procedure and tended to make the apparatus more 
complicated. 

Current interest in the measurement of the ad- 
sorption of nitrogen at low temperatures is in the use 
of the data to determine surface areas. The latter 
are based upon determinations of the volumes of 

trogen that are } cover the surfaces 
with a monolayer of the gas. Relationships to do 
this have been developed by Brunauer, Emmett, 
and Teller [5]. W.th certain assumptions as to the 
packing of the nitrogen molecules, the extent of the 

lsorbing surface may then be estimated 

The conventional gas-adsorption apparatus in use 
today consists of (1) an adsorption cell to hold the 
adsorbent sample, (2) a gas buret (usually of the 
multibulb variety), (3) a manometer, and (4) a 
sumping system that includes a mechanical pump, 
diffusion pump, liquid-air trap, and McLeod gage. 
The system is first evacuated to remove air, water, 
d other interfering impurities from the adsorbent 
he latter is maintamed at a temperature suffi- 
utly high to assist in the expulsion of desorbable 
purities, vet not so high as to bring about an 

versible change in its structure. Following the 

issing treatment, the sample is cooled to the 
perature at which the adsorption is to take place, 

a determination of the “dead space” is made 

i purified helium. This consists of the 


necessary to 


space 
ilted from a joint research project undertaken by the 
States Cane Sugar Refiners & Bone Char Manufacturers, a greater part of 
g industry of the British Commonwealth and Belgium, and by the 

| Bureau of Standards 
irch Associate at the Natic 
iting manufacturer 
es in brackets indicate the 


vestigation res 


mal Bureau of Standards, representing the 


literature references at the end of this paper 








the conventional volumetric method when the same 
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adsorbents are used with 


volume within the adsorption tube not occupied by 
the sample, that is, the space between the particles, 
and the volume of the connective tubing that joins 
the adsorption cell with the buret manifold. It is 
customary to admit purified helium to the buret 
system and to measure its volume and_ pressure 
before and after exposure to the adsorption cell con- 
taining the sample. The dead space is computed 
from these data by application of the gas laws 

It. is necessary next to exhaust the helium from the 
system before the adsorbate (usually nitrogen) can 
be admitted. Customarily, the sample is allowed 
to warm up during the evacuation to assist In the 
removal of helium. The procedure for introducing 
the gas is similar to that outlined for the introduction 
of helium. When a nitrogen-adsorption isotherm is 
Lo be used for determining surface urea according to 
the “free-surface” equation of Brunauer, Emmett, 
and Teller [5], it is important to confine all equilib- 
rium pressures to the region where the equation is 
applicable. This is usually valid below relative 
pressures of about 0.3, where the relative pressure ts 
the ratio of the pressure to the saturation pressure 
at the temperature of the experiment. A minimum 
of two points is required to determine the straight 
line resulting from a plot of p/|V(po—p)] against 
where p is the equilibrium 
pressure, p» the saturation pressure, and VV the 
volume of nitrogen adsorbed. However, it is good 
practice to obtain four points for check purposes 
The reciprocal of the sum of slope and intercept of 
this straight line is directly proportional to the sur- 
face area. In many cases the multibulb buret can 
be used to obtain successive points on the adsorption 
isotherm with only a single initial filling of the buret 

The simplified technique described in this paper 
[6] is based upon the selective adsorption of nitrogen 
from nitrogen-helium mixtures. Because of this 
phenomenon, the following advantages are possible 
(1) the entire evacuation system can be eliminated, 
(2) equilibria are attained rapidly by repeatedly 
eveling the gas mixture through the sample, (3) time 
is saved by not having to wait to remove the helium 
or to recool the adsorption tube before admitting 


p/po (relative pressure), 


nitrogen. 











The description of the new apparatus and the 
account of the new procedure will be given in some 
detail. The instructions are intended to be sufficient 
for use by others not too expert in this field of phys- 
ical chemical measurements. 


II. Design of Apparatus 


A diagram of the apparatus is shown in figure 1. 
The adsorption tube, 7, consists of a glass tube 
30 cm in length and 12 mm in diameter fitted at the 
top with a ground-glass joint. The inner member 
of the joint extends downward to within 3 or 4 em 
of the bottom of the tube, where it is rounded off 
like a test tube. The annular space between the 
outer tube and the extension of the inner member 
constitutes part of the dead space, which should be 
kept as small as possible [7]. For this reason the fit 
should be extremely snug. The sample is contained 
in the bottom of the adsorption tube. If the sample 
is finely grourd, a small pad of glass cotton placed 
both above and below the sample will hold it intact. 
The gas is brought in contact with both top and bot- 
tom of the sample by capillary tubing, which termi- 
nates in ball and socket joints. One lead connects 
with the annular space just below the ground-glass 
joint of the adsorption tube, and the other lead ex- 
tends from the bottom of the tube where it forms a 
U-bend, as shown in figure 1. The tube is purposely 
made long in order that the extremes in temperature 
to which the sample is exposed do not affect the 
lubricant of the joints and stopcocks at the top. 

The other half of each ball and socket joint is sealed 
to capillary tubing that communicates to the tops of 
two equal-volume cylindrical burets, P and S. 
These are made of 30-mm tubing and are calibrated 
to contain a volume of about 200 ml. The volume 
of P is bounded at the top by stopcocks J and E 
and at the bottom by the graduation mark A. The 
volume of S includes the space between graduation 
marks B and C. Leveling bulbs, Z and M, are 
joined to the bottoms of S and P, respectively, by 
means of rubber tubing. Air traps, not shown in 
figure 1, are located at 
and /? for the purpose of catching any stray bubbles 
of air that might be introduced through the rubber 
connections from the leveling bulbs or from gases 
entrapped in the mercury. A stopcock, D, located 
below graduation B serves as a positive cut-off for 
holding the mercury level either at mark B or C 
while a pressure measurement is being made. All 
pressure measurements are made relative to atmos- 
pheric. A sidearm made of 10-mm tubing is con- 
nected to bulb M and is used as one arm of the 
manometer with a vernier-adjusted slide on the milli- 
meter scale, ?. The other mercury level is in the 
buret system, and readings are only made when the 
mercury level is at mat’:ss A and B or C. Atmos- 
pheric pressure is determined at frequent intervals. 
The purified gases, helium or nitrogen, are introduced 
through stopcock -. Stopcock /, located between 
the primary buret ? and the adsorption tube, enables 
pressure measurements to be made of the gases in 


the lower extremities of S | 
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the primary buret before exposure to the sam |» 
Stopcock X serves as an exit to the atmosphere. 


Ill. Experimental Procedure 


The sample is first freed from interfering impurit os, 
such as air and moisture, by heating it in a con- 
tinuous stream of helium. The helium in = tanks 
(obtained from Bureau of Mines, Amarillo, Ty 
was of a high grade of purity and might be used as 
received. As an added precaution, it was furthe 
treated by passage through a trap containing ac- 
tivated charcoal surrounded by liquid air. The 
burets are kept completely filled with mercury so 
that all of the evolved gases are carried along with 
the helium through stopcock YX (fig. 1), where they 
escape into the atmosphere. An electric sleeve 
furnace fits around the lower portion of the adsorp- 
tion tube, and aluminum foil is placed in position to 
serve as a heat shield for the greased joints above, 
as well as for the burets at either side. 

After degassing, stopcock X is closed and helium is 
permitted to enter P by lowering M slowly untii the 
primary buret is filled. Stopcock £ is then closed, 
thus isolating a quantity of helium within the system 
The sleeve furnace is removed and replaced with a 
Dewar flask containing liquid oxygen.‘ Liquid 
oxygen was prepared according to the method de- 
scribed by Lane and Watson [13]. Care was taken 
to maintain the level of liquid oxygen by frequent 
additions. The mercury level in S is adjusted to 
mark C and stopcock D closed. The mercury level 
in P is brought to mark A, and the level in the side 
arm of .M is read against the vernier slide on R. This 
scale reading, When added to the atmospheric pres- 
sure and diminished by the null-point value of the 
scale, gives the pressure of helium in the system 
The level of mercury is next dropped from C to B by 
opening stopcock D and lowering bulb Z, thus allow- 
ing helium to fill the secondary buret. Stopeock D) 
is again closed before attempting to bring back the 
mercury level of the primary buret to mark A 
When this is done, the pressure is again measured 
These data suffice to determine the dead space 
However, it is recommended that the process be 
repeated for check purposes after first expelling a 
portion of the helium from the system. The dead- 


space volume may be calculated according to the 
equation: 
V/T=[pe/(pi—p2)](V2/ To) — Vi/ To. (1 


where )’=geometric volume of dead space 
\’, geometric volume of buret P 
V,—geometric volume of buret S 
7=temperature of refrigerant bath (deg K) 
T,=room temperature (deg K) 
p:=pressure of system when helium is con- 
fined to primary buret and adsorption 
tube 
pP2=pressure of system when helium is expande: 
to include secondary buret. 


* Liquid nitrogen may equally be used if it is desired to perform the measu 
ments at the lower temperature 
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The pressure of helium may be further reduced, if 
d, before continuing with the determination of 
the nitrogen adsorption by raising the leveling bulb 
\/ until the helium in the system is compressed 


above atmospheric, then opening momentarily the 
exit stopcock. The residual pressure of helium is 


finally measured with the mercury levels at marks 


Aand C. 


quation 1 assumes that all the volume between C 


and /isat the temperature 7. However, the volume 
not at 7’ G. e. not immersed in the bath) is deliber- 
ately made as small as practical. Moveover, the 
small error resulting from this assumption is almost 
compensated by the use of the identical procedure 
when nitrogen is introduced. 

Durine all the preceding steps, stopcock / is left 
in the open position. It is closed immediately after 
the residual helium pressure is measured. The 
leveling bulb .M is raised so as to bring the helium in 
buret ? nearly to atmospheric pressure. With some 
experience this enables a better control of the volume 
of nitrogen to be admitted. Nitrogen gas for the 
adsorption measurements was purified by passing 
tank nitrogen first over heated copper (450° C) 
supported vertically in a column, and second by 
through tubes containing Ascarite and 
Anhydrone. The purified nitrogen at atmospheric 
pressure is admitted to the primary buret through 
stopcock £, care being taken that the mercury level in 
P remains above that in .V/ so that the pressure in P is 
below atmospheric, and no helium is lost in the process, 
Sufficient nitrogen is introduced to yield a measurable 
increase in pressure over that of the residual helium 
after equilibrium is established with the adsorbent. 
Stopcock FE is then closed, the mereury level brought 
to mark A, and the pressure measured as before. 
The difference between this total pressure and the 
residual or partial pressure of the helium previously 
measured gives the pressure of nitrogen introduced 
into the buret system. Stopcocks / and D are now 
opened. By raising bulb VV and lowertiig 1, the gas 
mixture is caused to pass from buret ? through the 
adsorbent 7 and into buret S such that the latter is 
filled with gas at approximately the same rate thai 
Pisemptied. When P is completely emptied of gas 
and S filled, the positions of the mereury levelling 
bulbs are reversed, which causes the gas mixture to 
pass from S through the sample and back into P. 
This constitutes one complete cycle of exposure. 
\fter several such cycles (four or five) it is usually 


passage 


found that the pressure (measured with mercury 
levels at marks A and (’) no longer decreases. This 
s recorded as the equilibrium total pressure. The 


librium pressure of nitrogen is obtained by sub- 
ting from this value the residual pressure of 
im (corrected, if necessary, for changes in room 
perature) The isolating stopcock / is’ then 
d; a second increment of nitrogen is drawn into 
the primary buret and the resultant pressure meas- 

las before. The cycling exposure to the sample 
s repeated according to the technique already 

ribed, and the second equilibrium total pressure 

‘tained. Although two points on the isotherm 
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are all that are necessary to determine the straight 
line used to compute the surface area, it is recom- 
mended that three or preferably four points be 
obtained to make evident an unusual error of 
measurement. 

The determination of the liquid oxygen or nitrogen 
temperatures and the necessary vapor pressure of 
oxygen or nitrogen were computed from published 
results.” The vapor pressure of nitrogen in milli- 
meters as a function of the absolute temperature 
(7) in the neighberhood of its boiling point is given 
by eq 2 [8] 


_ 
~ 
to 


334.64/74-7.5777 


log Pr, 
The vapor pressure of oxygen is given by eq 3, 
which is valid over the same range of temperature [9] 

log po, 419.31/7+8.1173—6.487 x 10-5. (3 
Utilizing eq 2 and 3 for adsorption at liquid oxygen 
temperatures, it was convenient to plot px, as a 
function of po,, with the corresponding temperatures 
marked along the resultant curve itself. 
pressure of the oxygen was taken as 
pressure; hence, it was necessary only to 
atmospheric pressure along the abscissa (jo, axis 


om. 
The vapor 
barometric 
locate 


to obtain the temperature of the bath and the corre- 
sponding vapor pressure of nitrogen. As saturation 
pressure is of the order of magnitude of 2,700 mm, 
equilibrium partial pressures of nitrogen up to atmos- 
pheric during an adsorption determination all fall 
within the range where the free-surface equation of 
Brunauer, Emmett, and Teller [5] is valid (relative 
pressure is less than 0.3 

The pressure readings made from day to day were 
corrected for ambient temperature differences by 
taking into account the difference in thermal ex- 
pansion between mercury and the graduated scales 

Correction must be made for the nonideal behavior 
of nitrogen above the sample in the adsorption tube 
An approximation of sufficient accuracy for adsorp- 
tion measurements of this type is obtained by multi- 
plying ihe observed volume of dead space \ by the 
correction factor [1+ ap], where p is the pressure 1 
millimeters, and @ is a constant for a given gas at a 
particular temperature [10] 

This same factor may also be derived from the 
the second virial coefficient or fiom 
viscosity data [11]. This treatment has been sum- 
marized in an unpublished report by William 5S 
Benedict [12] These results are compared in table | 
with results quoted by Emmett and Brunauer [10] 
The correction must be made, since its magnitude 
may be as large as 5 percent, but the relative dis- 
crepancy among the correction factors from the 
various sources is small even at atmospheric pressure 
For example, the disagreement between the two 
extreme values for nitrogen at 77.3° K and 760-mm 
pressure amounts to only 1.5 percent. 


evaluation of 


More recent data on the vapor pressures of oxygen and nitrogen have been 


in tabular form compiled from a critical study of all « 


Tarie | Volume correction factor a for nonideal gases 


Viscosity 
data [11 


Emin 
B 


tt and 


unauer 


IV. Results and Discussion 


The basic assumptions of the proposed method are 
that nitrogen is quantitatively adsorbed from nitro- 
en-helium mixtures, and the helium remains in the 

as phase. The most direct test for the validity 
of these assumptions is to compare the results of 
surface areas determined by means of the simplified 
technique with results for the same samples by the 
conventional method (where no helium is used dur- 
ing the nitrogen adsorption Samples of materials 
varving widely both in composition and surface area 
were for these Furthermore, the 
residual pressure of helium in the nitrogen-helium 
mixture was varied in repeated runs with the same 
sample and in one case (Bone Char 32) exceeded 
half an atmosphere. Although such practice is not 
recommended in routine measurements, it was done 
here merely to emphasize the complete lack of 
dependence of nitrogen adsorption upon the pres- 
ence of helium. The results are shown in table 2. 
Liquid oxygen was used as refrigerant except where 
indicated. The comparison of the values for surface 
area given in the last two columns of table 2 shows 
the good agreement between the two methods. 
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The possibility of introducing an appreciable error 
dlue thermomolecular diffusion effects was con- 
sidered during the development of the simplified 
technique. As the manometer is connected to the 
buret system, which is at room temperature, it is 
assumed that the same pressure exists above the 


to 
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sample in the refrigerated adsorption tube. 

assumption is not strictly valid, and the error 1 
conceivably reach appreciable proportions as 
total pressure is diminished. This phenomenon 
been investigated for individual gases by Wy 
(et al.) [14 to 16] from theoretical considerati: 
The general equation derived in Weber's investi 
tion was simplified and a relationship deduced 
estimating the lowest pressure that would still 
compatible with a pressure difference amount 
to no more than a predetermined percentage of 
pressure. This relationship is given as follows 

T 4 


(7) -@) 


thermomolecular pressure difference: 

manometer pressure in microbars 

radius of the connective tubing 
centimeters 

temperature of the warm end of the 
connecting tube (room temperature, 

KK) 
temperature of the cold end of th 

connecting tube (adsorbent tem- 

perature, K) 

273.1 °K 
parameters depending upon the gas 
employed. 

The values of &,, n, and ,\, for nitrogen are given 
by Weber as 1.40, 0.41, and 5.91 and for helium 1.32 
0.147, and 17.85, respectively. By setting Ap/p 
0.01 in eq 4 and using 7;=300°, 7,—90°, and 
R=0.05, it is possible to solve for p (converting 
pressure to mm Hg). For these conditions py, is 


found to be 1.37 


2+2n 


2+2n7) 


} 


6k ,()A_)” 
(1-4 n) ah? p? 


Ap 
pP 


where 
Ap 
p 
R 


n 


> 
>. 


mm Hg and py. 4.27 mm Hg. 

It is seen that an error no greater than 1 percent 
of the total pressure would result by virtue of the 
thermomolecular pressure difference (thermal diffu- 
sion) as long as the total pressure of nitrogen by itself 
or helium by itself remained greater than the values 
shown above. Although no theory bas yet been 
advanced to predict the thermomolecular pressure 
difference in mixtures of gases, it is reasonable to 
expect that for nitrogen-helium mixtures the erro! 
should be of the same order of magnitude as for the 
separate Furthermore, as the free surface 
equation of Brunauer, Emmett, and Teller is usually 
valid within the relative pressure region between 0.| 
and 0.3, and since po for nitrogen at liquid oxygen 
temperatures is of the order of magnitude of 2730 
mm, the partial pressures of nitrogen need not by 
less than 270 mm. At this pressure the error due to 
thermomolecular pressure difference must certain!) 
be negligible. Even at liquid nitrogen temperatures 
it can be shown from eq 4 that the pressure limit for 
l-percent error is virtually the same as at liquid 
oxygen temperatures. A partial pressure of nitroge! 
of 76 mm corresponds to a relative pressure of 0 
Here, too, the correction is negligible. 

The final step in completing each cycle of exposu 
prior to reading the pressure consists in passing t! 


vases. 





yas mixture from S through 7 into P 
flicient 
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HELIUM 
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Simplified apparatus for adsorption of nitrogen 


from helium mixtures. 


(see fig. 1 


heat must be conducted away from the 


gas mixture on its way to the sample to prevent a 
rise in temperature of the latter during the half-cycle 


pre ceding 


requirements may 


measurement. The heat 
be estimated from the volume of 


each pressure 


vas mixture to be cooled, its heat capacity, and the 


temperature change 


Irom 


ni 


nly about 50 cal 


If the buret system were com- 
tely filled with gas at atmospheric pressure, 200 
would correspond to about 9 millimoles. Its heat 


ipacity would be of the order of 6 cal/mole (7 eal 


helium 
cooling 
K would require a heat transfer of 

As the latent heat of vaporiza- 
of liquid oxygen at its boiling point is 51 eal, 
it 1 g of oxygen is boiled away in each half-cycle 
xposure. The level of liquid oxygen was main- 
d by frequent additions 


if all nitrogen; 5 cal/mole, if all 
a temperature change of 210° in 
300° to 90 


irthermore, the gas mixture must pass a distance 
0 em below the surface of the liquid oxygen 
igh the annular before it reaches the 
le. The annular ring is 1 em in diameter, and 
If the time 
red to complete one half-cyele is 3 mins, it can 
hown from heat tranfer considerations that the 
is It reaches the sample, is virtually at the same 
erature as the bath. 


space 


vall thickness of the glass is 1 mm. 


feasibility of using nitrogen-helium mixtures 
een amply demonstrated here, but by the same 
other adsorbates mixed with helium might 
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[18] The 


also be used. It is believed that the greatest single 
achievement of the simplified technique is the com- 
plete elimination of the vacuum pumps and_ the 
associated equipment. It is now both feasible and 
practicable to house the entire assembly in a com- 
pact unit and thus make portable 
apparatus 


possible a 


V. Summary 


1. A simple apparatus requiring no vacuum system 
has been developed for the determination of surface 
area by the volumetric adsorption of nitrogen at low 
temperatures. 

2. The presence of helium, even at a relatively 
high partial pressure, does not affect the adsorption 
of nitrogen insofar as surface area determinations 
are concerned 

3. Comparison of the results obtained for surface 
area by the simplified technique with those obtained 
by the conventional method shows excellent agree- 
ment over a wide range of adsorbents 

4. The effect of thermomolecular pressure dif- 
ference was investigated and found not to be appre- 
ciable under the experimental conditions of the 
measurements 
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Transformations to Speed the Convergence of Series' 
J. Barkley Rosser * 


Numerical instances are given of the speeding of the convergence of series by the Euler 
transformation. This is even applied advantageously to certain divergent series, and a 
rigorous justification is given. An example is given of a series for which use of the Euler 
transformation is not useful. Instances are given of several less widely known methods. 
Finally, the method of summation by transformation into a continued fraction is illustrated 
successfully in the case of certain divergent series. The possibility of applying two different 
methods in succession to a given series is exploited throughout the paper, in spite of the fact 
that this often requires summing a divergent series. 


A remarkably useful such transformation is based Applications of the Euler Transformation with 
on a formula due to Euler: r=1 give 
¥ 1 r r n =. % 1)” ! 7 | 
>" a zr" Ss (A"a )( ) D (1) log 9 “pe Ss. , 
n= l—z £=0 “Al—z ” n=l n nai n2" 
We shall refer to this as Euler’s transformation. ee, ym 1S n! 
A purely formal derivation of this is as follows. ee fp 2n+1 2A=91-3-5---(2n+1 


Recalling that 


a,=(1-+A)"a,, 

One recognizes the right side of eq 2 as the expan- 
sion of —log (1—1/2). 

An interesting variation of the Euler Transforma- 


we write 





x . x l . ° . 
2, Gat" = >. 2%(1+4)"a, ~ a, tion is obtained as follows. 
n=) n=() | r(l T A) 
N—-1 D 
| | | Syagz*= 3 ant" +2*% Do Onin" 
le le n=() n=(0 m= 
I—z—zA i—z 7 
j— , A V1 x r m 
b 4 S° nm S* m ? 
a,xr"+ (A ay)( ) (3 
n=0 l1—z m=0 l—r 
l r P 
l—zraow \l—2z oe Applying this gives 
° ‘ . ° . > (—_,y)" N 1 / r)" 
A rigorous derivation of eq 1 is given in [1], pp. | Jog (4  % beam. = 
> 2e “a : 7 ; | log (li+27)—2 2, IZ 
62-66. It is there stated that if the a’s are all nao n+) nao n+ 1 
positive, then the right-hand side of eq 1 converges ( \V a . \m ont 
for every negative value of z for which the left-hand | tose. > ( = ) : ode 
side converges. That this is not necessarily true I+z2 mo\l+2/ (N+1)-(N+m-+! 
for positive z is easily seen by taking r=1 and 
Taking N=8 and r=1 gives 
1 
a» : 
(n+ 1)* P08. S20. 2. eg J 
log 2=1— — — z_— 
a. : | = 913 475 677 bd 
Actually, it is clear that the right side of eq 1 is | 
not likely to be an improvement over the left side | ee 6. . p.-. 1 1-2.3 ) 


+5)qt + + sates f 


unless ¢ is negative and the a’s are positive. In 919°99-10' 49-10-11 8 9-10-11- 
other words, the Euler Transformation might be sian Sas , 
expected to speed convergence only in the case of | We vet 

alternating series. However, for alternating series, 
the Euler Transformation is usually very helpful. | eae _l_9 63452 38] 


The Euler Transformation is a special case of a 2 
more general transformation given by Markoff (see 
[9], pp. 178 “194 ). However, we confine our atten- Also (see table 1), the terms inside the curly bracke ts 
tion to the simpler Euler Transformation. add up to 0.11724 674, and we get 


The preparation of this paper was sponsored (in part) by the Office of Naval 
Research 


? Help in carrying out the computations was furnished by Nancy Mann and ¢ = -9°@ l —< a yt - 
Shirley Marks under the direction of Gertrude Blanct log 2 0.63452381+, (.011724 674) =0.6931471s 
Figures in brackets indicate the literature references at the end of this paper - 
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TABLE 1 


O.LLLLE 


0.11724674 


{ more accurate value is log 2=0.69314 71806. 

One of the great advantages of the Euler trans- 
formation is that it can easily be used for numerical 
computations even in cases where one has only 
numerical values of the coefficients to work with: 
indeed, the computations are very direct even in 
such cases, because only a simple differencing opera- 
tion is called for. 

\ rather sensational use of the Euler transforma- 
tion is in connection with asymptotic series. We 
give an instance 

We have for y>0 


t 
€ 


‘ dt 
t 


y °o 14 


If the series within the square brackets were con- 
tinued indefinitely, it would diverge strongly. Thus 
it does not appear that the Euler transformation is 
applicable in the present case. However, let us pro- 

| formally to try to compute 


ise of the series alone. 


We get 


We have stopped just before the least term 24/625 
and if we add up the four terms shown the 
less than the next term. That is, we get 0.83200 
with an error less than 0.03840. Comparing with 
a more accurate value, 0.85211, we see that our 
error is indeed less than 0.03840 

Formally, we have neglected the terms 


ciror 1s 


and properly so, since they begin to diverge rapidly 
Nevertheless, let us now try to sum them by use of 
the Euler transformation, taking + l, 


24 


ly ; 
O20 


ete. 
We readily form the differences of the a’s numeri- 
cally (see the appendix for some tricks for shortening 
the computation) and get 
d)—= 0.03840 00000 
Aad, 0.00000 QOOOO 
A‘ay=. 0.00768 00000 
A*ag= 0.00307 20000 
Atay = 0.00645 12000 
Ady = 0.00761 85600 
Aa, = 0.01406 97600 
0.02602 59840 
0.05613 40416 
0.13145 60410 
0.33766 21486 
0.93823 63791 
2.80697 67610 
S.9SS51 15364 
30.66826 95734 


111.03898 71074 








Kuler 


vives the SCTICS 


So the transformation, if applied formally, 


10) 1920 OO000— 96 OOOOD 1Y 2POOOO +20 16000 


1} 90400 10 99200 10 16640 10 96368 


12 8$3750+-16 48741 —22 90616+4-34 26485 


jd SOHL52 OS 5OY 1 4 160 4 POD 


This is actually a divergent series, and the last 
terms listed are already increasing. Nevertheless we 
treat it like an asymptotic series, and add up all 
terms preceding the one that is numerically least 
This 0.02016 04800. Adding this to the 
0.83200 that we had already. we get 0.85216 04800. 
which compares favorably with the more accurate 
value 0.85211 
We neglected the terms 


vives 


107". 10 16640— 10 96368 + 12 83750—16 


tS741 


22 90616— 34 26485 + 54 86152—938 59213 


169 43205 
We now apply the Euler transformation to these 


terms, taking sr l 10 16640, 4 LO 96368, 
This gives 


-anda 


107") 5 08320— 19932 + 13457 —4372+-916— 1383 


50 5704-6 


This adds up to —0.00004 96486 if we stop before 
the least term. Adding this to our previous value 
gives 0.85211 08314, which compares very favorably 
with the more accurate value, 0.85211 08814 

The agreement with the more accurate value is 
experimental evidence that use of the Euler transfor- 
mation is justifiable even in the present case where 
the series diverges Actually, we can justify use of 
the Euler transformation rigorously for the present 
series. Looking at eq 4, we see that the remainder 
can be written 


l 


>|, Gy) 


which can be expanded as 


’ \ if *. PF \ 
») J. (y) ’ 


if (x) ("-1). 
ted (yw) CG) 


We then see that all differences of a, of even or 
are positive, and that if we stop with a term involy 
a difference of odd order, the error involved Is 
than the next term (which involves a positive dif 
ence of even order 

To justify a second application of the Euler tra 
formation, we modify the remainder term gi 
above as follows 


LG@ & 


1) 1 )e~‘dt 


We can apply the Euler transformation to oth: 
asymtotic series with equal success. For insta 
we have 


(we) E-) GG 





n justify applying the Euler transformation to | We find 


3\ mptot mw series 


applied the Euler transformation to this series 
r=6 to compute a 20 decimal value of 


2-day, 


h was then checked by a different means of 
putation (see [5], p. SO 
e can use the Euler transformation to derive 
v's “Converging Factor” see 2] In those 
es (as above) in which we can justify using the 
transformation on an asymptotic series, we 
then afforded a means of justifying the use of 
ev's “Converging Factor.”’ 
To see how this is done, let us write 


So use of the Euler transformation gives us 


If we rearrange this in powers of l/yv (such a rear- 
rangement needs justification), we find that it ts 
identical with the formula on p. 526 of [2] 


If we choose 1 N, y V, so that A=0, B 
“Converging Factor’ reduces to 


Lirey refers to the series in the curly brackets as 
he “Converging Factor” and gives a purely formal 
evice for estimating a value of it (Airey’s estimate | with C,=1, ©, »_ Gy , Cy=13, ete 
s justified in 112] We will now show how to derive values of the first 22 C’s are given on p. 529 of 
\irey’s formulas by use of the Euler transformation, | Using the first 17 of these with N=5 we get the 
Which we have already shown to be justified for the | following series for the “Converging Factor” 


The 


») 


present series We write 
| 


») 


10-41 00000 GO000 00000 5000 00000 OO0O00 


250 00000 00000— 12 50000 00000 +8 1250000000 


en the “Converging Factor’ takes the form te 
1 46875 OOOOO 11406 25000 19117 IS750 


6566 7O6BSS 503 T3047 1686 04590 — 927 48286 


+26 SOSS5S +364 18569—-295 53896 +41 19945 
154 50012— 177 02724 


This series 1s probably divergent if earried fat 
enough, but we sum it as if it were an asymptoti 
series and stop just before the least term shown 
This gives the value 0.52372 086 Airey gets0.52372 OST 
[2] p. 529), but the difference of 1 unit in thi 
last place is likely due to differences in rounding pro- 
cedure. This should be multiplied by 0.0384 and 
added to 0.832, which gives 0.85211 OS810 as com- 
pared with the more accurate value 0.85211 OSS14 

In [2], Airey applies the “Converging Factor’ to 


sec 


Various asvmptotic series without rigorous yustifi ii- 
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tion, but with great numerical success. By deriving 
the “Converging Factor’ by means of the Euler trans- 
formation, one could justify its use in some of the 
cases given in [2], though there would still remain 
the problem of justifying various rearrangements 
made by Airey of the seriesin the “Converging Fac- 
tor.” Some of these can be justified by the pro- 
cedure in {12}. 

Although we have seen that the Euler transforma- 
tion can be applied in some rather remarkable situa- 
tions, nevertheless there are cases where it cannot 
he applied. In particular, for nonalternating series 
(for instance, when z and the a@’s are all positive in 
(1)), it is of almost no use. However even for alter- 
nating series it occasionally fails, and we will now 
give an alternating series for which it is of no value 
whatever. 

Let us expand 


by the Euler-Maclaurin sum formula (see [4], p. 128) 
using an interval of length unity. We first need the 
values of the derivatives of 

1 +2" 
of odd order evaluated at Oand 1. At 0, the deriva- 


tives of odd order are zero. To get the values at 1 
we write 


Putting r=1 and recalling 


we get 


n+l)q@i n+I)q@t 
d" l (—1)"n! ) ‘ 4 d 
dz" 1+27 Jens 24 2)" (‘ ) 


Substituting in the formula 


a © ws Bl ses 
F(r)\dx=5 F(0)+5 FO) 


. 


(—1)"B,, 


(2m)! 


F2=-)(1)—F2"-(0)} +2, 


(see [4], p. 128, with w=1, r= 1) we get 


"1 dz 3 #-! (—1)"B, . xm 

a m R, 

Jol+x? 4 <4 m2”?! sin 5 +h 
3,  B,.(—1)Y 

Ss 2 I R 

at 44ers parti - 
5 l l l l 
4*24~ 20164224 1536 


ease 29608. abe 
435 867 = ‘Ss 657081 _ nie 
7354368 565248 © 98304 

This is certainly an alternating series. However, 
the partial sums do not lie alternately above and 
below the true value. In table 2 we have given the 
terms of the series with the partial sums to the right 
If we try to improve the value by use of the Euler 
transformation, we get 107~°(79117+- {12—10+61 


758+-- - -}), which is certainly no better. 
TABLE 2. 

Series Partial sums 
+0). TA000 0. 75000 
+4167 79167 
oU 79117 
+24 79141 
65 7OO76 
+506 79672 
0. 13743 65929 
+. 69282 7. 38211 
True value 0. 7854 


As indicated, the Euler transformation is valueless 
for speeding convergence of a nonalternating series 
Occasionally one can use a trick to transform a non- 
alternating series to an alternating series. Thus the 
formula 


n+l 
en l 9 ~ 1) 
amd 2 hm 2 
n=} / n=l hn 


by (indirect) use of the Euler transformation. How 
ever, in general such tricks are not available. 
Fortunately the Euler-Maclaurin sum formula ts o! 
great value for summing many nonalternating series 
For instance, by use of the Euler-MacLaurin sur 
formula (see [4], p. 128) we prove that for m>1, s 
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— ] m m 
\ > > + 
—_ i vi s l 2 
= s(s 1) (s PA 2) ( 1) B, 
r=1 (2r)! "deli 
/ lhe series on the right is divergent, but it is an 


asymptotic series, and by taking m large compared to 

msiderable accuracy can be obtained, which can 
be further increased by applying the Euler trans- 
mation to the asymptotic series. Thus, consider 


{ 
0 


the ease Where s=2. We have the well known result 
oe Se 9402 ANA 
ant 7 1.64493 40668 


») 


Even with m as small as 2 we get 


Aer I 
1'2 §'48 960 5376 


] 5 
~ + —— — 
15360 1 3516S 
691] 7 
+ 
64160 © 1 9660S 


3617 $3867 
668 46720 © 4183 81824 


»~Y 


el 1 74611 8 54513 
+ 
nd 6920 60160 © 11576 27904 
he 
: This is divergent, but if we sum all terms before 
t . . 
the least, which is 
691 
223 64160° 
we get 1.64495 as compared to the more accurate 
value 1.64493. If we now apply the Euler trans- 
formation to the terms neglected, we get the series 
LO 154488 11765 7249 11523 
14389—20493 
Doubtless this diverges, but if we sum all terms 
before the least, we get 1.64493 36 as compared with 
the more accurate value 1.64493 41 Another use of 
the Euler transformation on the terms neglected 
SS above gives further improvement. 
‘S In [5], Bickley and Miller give a method for deal- 
n- with series Sa, in which 
hie 
a, A A Ay : 
] + = + 4 7) 
a» i nn iN 
In those cases in which a, is a differentiable fune- 
tion of n, one can treat such series by the Euler- 
Maclaurin sum formula, and usually one gets the 
same approximation as would be given by the 
B ley-Miller method. 
the Bickley-Miller method, one writes 
\ 
’ \ 
S > > a» 
() n=) 
. , 
’ x 
Sm > ay 
n 0 


Then one tries to determine constants 
so that 


a 


\ a 
a@_)M-+ a+ + 
( m 


a , / 


m* ) 


Om» 


In general, this is impossible, since the series 


ay 
+ 
m 


will be divergent for each m. However, usually this 
series is an asymptotic series, so that for large m 
great accuracy is attainable. 


To determine the a’s, we proceed formally as 
follows. By eq 7, 
’ . \ a 
S S - , 6. a@_,\m 1) Qt 
mn l 
; a / 
4 4 
(7A 1)- ) 
Subtracting eq 7 from this gives 
@a=Sa—Sa-1=—Anu-1 A» 
However, by eq 6, 
A, , As 
@n=Gn-1( 1——'+—3+ --- ) 
Mi me 
Substituting this and factoring out a,,_, and dividing 
by m gives 
l Ay, La - a 12 4 a 
mm m m m m(m ] 
A, 1 \ ¥ 
pes 1 — Atay ...) fa ,44 I. 
m mm ( m ) 
If we now expand 
l l r rr 4 
> T 
m(m 1) m I'm’? 2!'m 
and equate powers of »we can determine the a’s. 
Mn 
Indeed we have 
a (A, l l 
\a l Ay a Ag, 
a(A,;+1 a+1)A a » Ais, 
ete 
| (Note error in sign in formula (12c) on p. 757 of [5)). 
If we apply the Bickley -Miller method to the series 
| wi 
—— 
ne} / 


61 














we get the same expansion that is given by the Euler- 
Maclaurin sum formula Bickley and 
Miller do not apply their method to any series of 
interest that cannot be handled by the Euler-Mac- 
laurin sum formula, but they indicate how such series 
might arise. Bickley and Miller do not give a rigor- 
ous justification of their method, but Szasz (in [10] 
has discussed the conditions under which one can use 
at least the term a@_,m in eq 7 for a first approxima- 
tion to S. It would be very helpful to have a general 
treatment of the conditions under which the Bickley- 
Miller method can be used 

In [7] and [8], Stieltjes gives some interesting ways 
for estimating the error in certain well known 
asymptotic series. However, his methods are rather 
special, and can hardly be applied to series at random 

We turn now to what is perhaps the most widely 
applicable method that we know. This the 
method of transforming the series into a continued 


sec abov ¢ 


Is 


fraction 
Let us return to the example given above in which 
we got an approximation for 


by the Euler-Maclaurin sum formula. We had a 
divergent of which we neglected all terms 
beginning with the least. Let us now try to find a 
sum for these neglected terms by transforming them 
into a continued fraction. To write them as a 
continued fraction, we first consider the series 


series, 


\ 691 l 7 | 3617 l 
223 64160 2 1 96608 668 46720 
{3867 | 1 74611 l 
~ 4183 81824 -°' 6920 60160 
§ 54513 l / 

11576 27904 . \ 


Transforming this into a continued fraction by the 
procedure given on pp. 196 to 202 of 16] we get 


a a > a 


with the coefficients 


30897 65053 


d,—= 0.0000 
a 1.1523 15485 
,=0.3674 32414 
a 1.7288 91714 
a,=0.8615 52894 
a;—=2.4273 31997 


(Our a’s after a) are taken to be the negatives of 
defined in [6]. Due to a gradual loss of signi! 
figures, the final digits of our a’s are increas 
inaccurate 

We then put 1, and use the resulting conti: 
fraction to compute an estimate for our remai 
Which, by subtraction from the known sum of 
series should come out to be —0.00001 54985 


successive convergents are 
0.00003 OS9O77 


$3555 


0.00001 


0.00001 53294 
0.00001 57006 
0.00001 54620. 
We note the characteristic continued fraction be- 
havior of the convergents in that they are successively 
above and below the true value, and monotonel) 
approach the true value from above and_ below 
Already the sixth convergent gives a value correct 
within less than 4 in the eighth decimal place, and 
seems very plausible to conjecture that later converg 
ents will continue to come ever closer, in spite of the 
fact that it is a divergent series that we are trans- 
forming into a continued fraction. This is a not 
uncommon phenomenon with continued fractions: in 
fact this is the basis of the Stieltjes method for sum- 
ming asymptotic series. In [6] on pp. 349 to 361, p 
p. 367, and pp. 372 to 373 are given man) 
examples of convergent continued fractions, which 
are the formal equivalents of divergent asymptoti 
series, including some that resemble very closely th 
series with which we are dealing. 
Note that we used only odd powers of 


365. 


>in our sertes 


\ 691 l } 
' 2a 64160 = \ 


> 
» 


This has the effect of greatly reducing the labor o 
finding the corresponding continued fraction (in thi 
notation of pp. 196 to 202 of [6], all the 6’s will b 
zero 

For an additional instance of the use of continued 
fractions, let us recall the series for the “Converging 
Factor” which we had earlier. We used all terms 
preceding the least, leaving severai terms unused 
Let now trv to attach a value to these unuse¢ 
terms by transforming them to a continued fractio! 
To avoid difficulties with the continued fractior 
expansion, we start with the term after the leas 
We formally transform the series 


us 


] (364 18569 295 


” ] 


53896 41 19945 
+ “ 


154 50012 
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» continued fraction where the \,, are increasing positive numbers, is 
sometimes helpful in estimating the sum of an alter- 
a a nating series of which S,, are the partial sums 


Appendix. Methods of shortening the 
computation of A””,. 


a 0.00000 OOTS2 O9O284 
In an earlier draft of this appendix, we set forth a 
ae 
a 0.81150 624 scheme for shortening the computation of A”a 
Upon reading this earlier draft, Prof. J. C. P. Miller 
sais - ; 
a 0.670210 Les proposed the following very) ingentous scheme for 
me computing A”“a 
S 0.80673 455 In the text, we put 
(4 0.02419 195. 24 
ad = —y 
. ; 625 
Because one of the a’s Is negative, we ho longer vet 
he familiar performance of the convergents. Never- 120 
theless, they are apparently converging to the true a,=., 52? 
value 0.00000 00078 69 since they are 3125 
> 720 
0.00000 00182 09 - 
) De- 2 7; vo 
15625 
Vvely ~ 
0.00000 00100 52 
] ely 
l ele 
TOW - 
ae 0.00000 00052 40 
: and wished lo compute A™a To do this more 
id " 
0.00000 00079 44 quickly, we write 
r thy * : V . 
he 0.00000 00079 16 : U 1)! 
itis i» = Gy? 
ne ' N)* a 
i respectively. Using the last, we get a value of 
um- )S521 l OSS] } 25 for je 5 ey | 5 as compared “Ww ith W here N 5 and 
the more accurate value 0.85211 O8814 237. 
on. One other remarkable use of continued fractions 
hicl s given in [11], where the asymptotic series a=) 
ot . . _ | Oo --= 
‘ . PY “odo > rd Le 
he é ey dy + + t - es | we | be 
aT , 2 tw sw low 320 ar 
converted into the continued fraction 
a 1+ 
l 2 } ) Ss N 
Dor Oy ou ow Du 
2 2 ; ; 
ro s : Qs (14 -)(14 -) 
| his continued fraction was used to compute por- N \ 
1} ns of a table (see Section 26 of [11]) in spite of the 
that the continued fraction ts likely divergent » 3 
to the fact that occasional convergents will have Oey ( 1+ \ )( l+ \ )( 1+ \) 
remely small denominators. However, as shown 
itn the estimates in Section 8 of [11], whenever we oie 
lave a convergent with a large denominator, this 
3 vergent is a close approximation to the value of | pen 
9 lunction 
ln closing, we might remark that some of the A" N—1)! A 
. . 1 xs 
isformations studied in theories of summability ; N _ 
occasionally useful to speed convergence Thus 
use of Riesz sums, So we seek a formula for A"a 
ee Mets Actual subtraction gives 
o . . 
‘ hi 
\ + > Xr, Aa > Aa. 
- \ 
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Then we prove by induction on m that ae 10320 , 64224 26432, 2520 
7 N° N’ N’ N 
/ yn mw . 2 200OnR 2499972 ‘ [22 - 
A" tla. , y A™"a, 4 * A™—'a,. Ae 3 62880 _6 23376 ‘ 3 03660 , 44100 945° 
j N® N® \ N‘ N* 
Pabtinn ou gives Ate, 36 28800 , 66 36960 , 36 78840 , 7 05320 , 34650 


Ne * N® V3 V Ne’ 


Moan ” , 399 16800 , 76998240 | 473 24376 . 110 98780 


8 66250 . 10395 


Then rN Ne 
- 48 4790 01600 | 9675 24480 | 6475 36032 
A™*'a, y A” a a a _ N® Vu Nw 
: 1773 31440 | 188 58840 5 40540 
N*® N® Nz’ 


Reealling that it is actually 
62270 20800 | 1 30967 36640 , 94189 45536 

















A™a , NB NB yu 
9" 29205 25608 , 3894 49060 | 182 88270 , 1 35135 
Nw NY V3 . 
that we seek, we write the simple recursion 8 71782 91200 19 00603 35360 1454105 80224 
| Aa vu t Nu + NR 
aera m (2° doy 4 a ). 499528 62960 , 79349 27000 | 5200 59540 . 94 59450 { 
gn+2 4N \~ 2" 3° yu ro Nw vv Te - 
of 
2 . the 
Using this, we readily compute numerical values of References - : 
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Random Walks and the Eigenvalues of Elliptic 
Difference Equations 


Wolfgang Wasow 


1 


Chis paper is concerned with the application of Monte Carlo methods to differential 


equations of the form 2u+g(zr,y)u=0. 


is shown that under suitable conditions the 


Monte Carlo solution converges and gives the solution of the difference analoz of the above 


equation 
I. Introduction 


(he aim of this paper is to describe and investi 
rate a random walk procedure that can be used to 
appreaimate the solutions of elliptic partial differen- 

equations containing the unknown function 
itself and not only its derivatives. The method also 
leads to a scheme for the numerical determination 
of the lowest eigenvalue of such differential equations. 

\s a computational technique, our method is 
somewhat similar to the one used by Donsker and 
Kae [3]? for the calculation of the lowest eigenvalue 
of Schroedinger's equation, but the underlying 
theory is more elementary than the theorem on 
Wiener integrals of [4] used in [3] 

Like all computational methods based on random 
sampling, those described in the present paper re- 
quire the use of a high-speed calculating machine. 
Numerical tests are in progress. 

The random walks considered lead to difference 
equations. By virtue of known results [1], [6], the 
solutions thus obtained are, for small step length 
approximations to the corresponding solutions of the 
limiting differential equations. 

Let the symbol A denote the finite difference ana- 

vue of the Laplace operator, 1. e 


h Y h Y 


h)+u(a,y—h fu(r,y)]. 


Then we shall be concerned with the difference 
lation 
Au+q(r,y)u=0, (] 
re g(r,y) is to be sufficiently regular to guarantee 
vergence of the solutions considered to those of 
corresponding differential equation. 

Everything that follows can be easily extended to 
e than Extensions to other 
tie differential equations are also possible. 

the sequel, the word “point’’, without further 
ification, is meant to refer to the points of a 


ire lattice with mesh length / 


two dimensions. 


parat 


II. Random Walk Procedure 


Let B be the interior of a finite domain whose 
boundary points form the set C. Consider the 
problem of finding the solution of eq 1 in B which 
assumes prescribed values /(S) at the points S of C 
This problem will be shown to be related to the fol- 
lowing random walk procedure. 

Let k(P) be a positive function defined in B A 
particle of mass 1 starts from an inner point P. 
Before moving to one of the four neighboring points 
P;, (g=1, ... , 4), its mass is multiplied by the 
value of k(P) at P. Then it moves to a neighbor- 
ing point, all 4 points having the same probability 
4 of being chosen. After a certain number of steps 
the particle arrives for the first time at the boundary, 
say at the point S. Consider now the random 
variable which is equal to the product of the mass 
upon arrival at the boundary by the value of f(P) 
at that point. We claim that, for appropriate 
choice of k(P), the expected value of this random 
variable—if it is finite—is the solution of eq 1 as- 
suming the values f(S) on C 

In the description above, the infinite random walks 
which never reach C have been ignored. This is per- 
missible, since the total probability associated with 
such walks is zero, (see [1], p. 44) 

Without loss of generality we may restrict the dis 
cussion to the special boundary values 

ro * for S=R 
RS 
(0 for S¥R 


where R& is some fixed arbitrary point of C. For, if 
we denote the expected value upon arrival corre- 
sponding to these special boundary values by /(P, R), 
then the expected value for any boundary function 
f(R)is >) E(P,R) A(R). If ECP.R) is indeed the solu- 

Ret 
tion of eq 1 assuming the values 6(R,S) on the 
boundary, then it follows immediately — that 
SCE(P.R) F(R) is the solution with boundary values 
Ret 
fUS) 

Theorem 1 Net 
h 
gP 
(1-79) 














and let h be so small that k(P)>0 in B, then the fune- 
tion ke(P.R) described above satistes if it is finite 
with respect to P the difference eq 1 and the boundary 
condition 2 

Proof: Let {L,(P,R); be the set of all possible 
paths starting at P and reaching C at R on the nth 
step. The probability of the particle moving along 
any such path is4~". Let m,(P,R) be the mass upon 
arrival corresponding to a given such path and write 


¢(P.R) i = 


i it Pe 


m, (PR) (4) 


R) | 


the summation extended over all paths of {L,(P,R) 
This quantity ¢,(7,2) is the expected mass upon ar- 
rival on C, if only arrivals at R and on the nth step 
exactly are counted with their actual mass, and all 
other random walks ending on C are considered to 
contribute the mass zero. If we set 


ee (P,R) (5) 
i ” 


BAP R 


then 


k(P.R)= lim E,(P, R) (6) 


provided the limit is finite. 

Let P,, (j=1,- ~~ ,4), be the particular neighboring 
point of P reached by the particle at the first step, 
when following the particular path L,(P,R): then 

m,(P, R)=k( Pym, —(P,, R), for n>0, (7) 
if m,—,(P,,R) is the mass upon arrival at R corre- 
sponding to that path of the set |Z, —,(P;,R)} which 
is part of L,(P,R). From eq 4 and eq 7 we have 


l a , 
e,(P, R) ym » eP)m,_(FP,,f%), a>6, ©) 
j=l il PR 
and, therefore, 
bP) 4. 
e,(P, R) - » P,,R), n>O0 (9) 


By means of definition 5 this yields 


BP) 4. =» 
E,(P. R) F P E, iCP;, R) (10) 
or, upon passage to the limit, 
: OP) 8. a 
E(P,R)="~" DS E(P,, R). (11) 


4 
This is readily seen to be equivalent with 


ab+i* DE 


hk Oin B 


For a boundary point S it follows from the defini oy 
of ¢,(P.R) that 


€( S, R) 6(S, P) } 


if 6(S,R) denotes the function which is 0 for \ «<R 


and 1 for S=R. For n>0 we have e,(S.R) -9 
These facts and eq 12 lead immediately to the s 
ment of the theorem. 

Remark 1. The function £(P, R) can be considered 
the finite analog of the normal derivative of Greon’s 
function, since the solution that assumes the bound- 
ary values /(R) can be written in the form 
par PR) {(R) which recalls the formula 

y , , 
oat _R) f R)ds 
JC on 


from the theory of partial differential equations. 

Remark 2. \t is well known that the probability of 
a particle leaving the domain B at a preassigned point 
R of C satisfies the difference eq 12, if the particle is 
subject at each step to a probability of “dying” equal 
tol—k(P). The corresponding function g(P) is then 
every where nonpositive In the next section it will 
be shown that the procedure of the present paper is 
applicable under considerably milder restrictions on 
g(P). A further advantage of our method is that at 
each step only one random decision has to be made 
and not two as in a random walk involving the possi 
bility of diving. 


Ill. Validity of the Procedure 


The statement of theorem 1 includes the assump- 
tion that (P,P) is finite. In this section we shall 
prove that this is the case, if and only if the eigenral- 
ues of Au+g(P)u in B are all positive. Here and in 
the sequel we mean by the eigenvalues of a linear dif- 
ference operator L[{u] in the domain B those values of 
\ for which the problem 

Llu) +d\u=0 in B,u=0 on C 
has a nontrivial solution. The eigenvalues in this 
sense are the negative of the eigenvalues of the matrix 
formed by the left members of the system of linear 
equations represented by L[u]J=0 in B, u=0 on ¢ 
It will simplify our terminology if we call the negative 
of this matrix the matrix “belonging to L{u].”’ 

The proof of the theorem mentioned requires son 
preparations. 


wee 8 
Lemma \: The eigenvalues of —Au 4 ( I~} Ju are 
3 


the same as those of Au-+ gu. 
Proof: The matrix A belonging to 


Au+gu=0 in B, 


is of the form A=G+H. Here, @ is a diagonal 


4 
matrix formed with the values of — g(P) +; » (P in 
( 
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porn 
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atrix 77 has in the row belonging to any given 


5 ea L. 3 alas 
P. the element > In column corre 


h 
ling to an interior neighboring point of P, and 
everywhere else. The matrix A* of the 


every 


sSVs- 


S , _ ° 
0 tem —Au- ( I=} ) u=0 is, similarly, given by A” 
¥ 


(}—/1. We may interpret these symmetric matrices 
5 is belonging to quadratic forms in the variables u(?), 
(d- P in B). In the space of these variables consider 
« orthogonal transformation which consists in re- 

ng the value of u(?) by —u(P) in every second 
int P and leaving u(?) unchanged in all the other 

ts, in such a way that a point and its neighbors 
This can, e. g., be 


always treated differently 
ne by setting 


1 


u'(P 1)" u(P (14) 


where (7,y) are the coordinates of P This orthogo- 
nal transformation changes the quadratic form with 
“1 matrix G+-//7 into the one with matrix G—J//. These 
wo matrices have, therefore, the same eigenvalues. 


Now we define an analytic function @(P,R:r) of 7 


au (PR: pM ue? 
a 


r. R (15) 


4 : 
We assume throughout that g(P? j , and therefore 
: 


P)>0. To show that the power series in definition 
' has positive radius of convergence consider the 
all random walk problem in which g(?’) has been replaced 

rywhere by max g(/’), and denote by e)(P?,R) and 


PeB 
P) the quantities corresponding to e,(P,R) and 
P?) in this new problem. Then, by the definition 
” ; f ¢,(P,R), 
P,R)<k*"(P) p,(P, R)<k* 
; here pa(P.R) is the probability of moving along a 
; h that leads in exactly n steps from P to R. 
Hence, the radius of convergence of the 
( formula 15 is at least equal to 1/k* 
lhe formula 15 may, of course, be interpreted as 
finite analogue of a Laplace transform, applied to 
solution ¢,(P?,R) of the “parabolic” difference 
ition 9. Pursuing this analogy, we find, by com- 
ng formulas 9 and 15, for @ the elliptic difference 


tion 


series ih 


4 l 
k(\@ + go a) (1 )¢ 0. (16) 


() unction @ is that solution of the difference eq 16 
which assumes the values 6(?,2) on the bound- 


we denote by T° the diagonal matrix whose ele- 
s are the values of #(P) in B and by A the 
ix belonging to A¢-+-g¢@, and if we set, for abbre- 








viation, 75 (1 ) A, then T.l—A/ is the matrix 
t 
belonging to the difference eq 16. The solutions of 
eq 16 are rational functions of \ and therefore of 7 
It follows from definition 15 that r=0 is not a pole of 
this rational function 

For later use we mention that the poles of @ are 
simple. For if ) is a matrix such that /?=T, we 
have TA—AJ=D(DAD—rANDD The elements of 
the resolvent (DAD —XJ)~'! of the svmmetrie matrix 
DAD have only simple poles (see, @. g., {7}, p 26), and 
this property is not destroyed by the matrix trans- 
formation with the nonsingular matrix J. Hence, 
the elements of (.4—)J/)~‘ have only simple poles as 
functions of A, i. e., @ has only simple poles as 
functions of 7 


_ } ] ‘ 
Lemma 2 Che number j ( | ) is an elgen- 
ys 7 


l 
value of k(Ad+g¢@), if and only if r=r, is a pole of 
o(P,R:r) for at least one pair of points PR. 
Proof. If 7 risa pole of (P.R:r) for some P 


: } 
and RP then i2( 1 


is a pole of some element of the 


resolvent of the matrix belonging to k(A@4 

Hence, this number is an eigenvalue of that matrix 
4 
h? , 
k(A@+ge@) and assume that $(/,R;r,) exists for all 
P and R. In order to show that this implies a con- 
tradiction we make use of Green’s formulas for the 
difference operator Au. These formulas are (ef. [1] 


dQ) 


] 
Conversely, let (1 ) be an eigenvalue of 
i. / 


h _ fu j Wu? AS rAu S ri u 17) 
B+¢ B ( 
and 
hS ‘rau uAr S“Irl u) u(r)! i) (1S) 
B+ ( 


where u,, u,, ete., denote the first forward differences 
and where the expression I'(w) is defined as follows 
if w assumes the value w) at a given boundary point 
R and the values 1, .. uy(vS3) at the y neigh- 


boring points of R in B+ C, then 
I 19) 
u u ul vu : 
h 


Green's formulas 17 and 18 are valid for any functions 
that vanish outside of B+C. If we substitute in the 
identity 18 for ¢ the function ¢(P?,R:r,) and for « an 
{” corresponding to the 


eigenvector eigenvalue 


; ( | : ); it follows that ['(/ 
i~ ry 


all boundary points R. Substituting / 
identity 17 we see that the bilinear form 


E + sar ( | )}e > 0 20) 


vanishes identically for any choice of the function ¢ 
Taking for ¢ the function that is 1 at an arbitrary 


0, identically for 


for vw in the 


. 
itcl 
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interior point and 0 everywhere else, the identity 20 
is seen to imply for ¢° the first order difference 
equation 


4 l 
/ Uhl a+, (1 , ) Je 


Together with 


Oin BLU =—OonC. 


(21) 


k(AL’ +-gl +5 (1 ‘yu 0 


this \ ields 


hAU Om Bl Oon C, 


h)—2U (2.y)=0nm BU =O0o0nC., 
This last difference equation has the unique solution 
7’=0 as can be seen, e. g., by calculating its values 
from point to point from the boundary inwards. 
Since /” is, by assumption, not identically zero we 
have arrived at a contradiction, and the proof of the 
lemma is completed. 
and prove our theorem 
If gP in 


Theorem 2 j B, then E(P,R) 
lim £, (PR) is finite for all Pin Band all R on C, if 


” 


4 


and only if all eigenvalues of Au+gu are positive. 
Proof: The difference expression k(Au+gu) is, in 


general, not self-adjoint, i. e. the matrix belonging to | 


it is, in general, not symmetric. But its eigenvalues, 
i. e. the numbers A for which the problem 

Au+gu+dAk-'(P)ju=0 in B,u=O0onC (22) 
has a nontrivial solution are positive if and only if 
those of Au+-gu are In fact, let A be the me trix 
belonging to Au-+-gu and denote by DP) the diagonal 
matrix whose elements are the positive determina- 
tions of A1(P). Then eq 22 can be written in the 
form 


If we make the substitution 
u— Dr 
this is equivalent with 
DADer=wn. 


Hence \ is an eigenvalue of the symmetric matrix 
DAD and therefore real. If we write DAD” instead 
of DAD we see that DAD’ is the matrix of the 
quadratic form obtained from the one whose matrix 
is A by the transformation 23. Since this trans- 
formation not destroy the positive definite 
character of a quadratic form, and the same is true 
of its inverse, A is positive definite, if and only if 
all the eigenvalues of eq 22 are positive 


does 
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We are now ready to formulate | 


the eigenvalues of k(Au 


positive if and only if = ( 1— ) is positive fo 


By lemma 2, + gu) ar 


poles of ¢(P,R;r), i. e. if these poles are all gri 
than r=1, or negative. On the other hand, by 
definition of @(P,R;r) the expected mass E(P 
is finite, if and only if all poles of @ are numeri: 
greater than 1. Our proof will be completed, i! 
show that the occurrence of negative poles gre 
than or equal to —1 implies the existence of 1 
positive eigenvalues of Au+gu. To show this 


A* ; 
The 


short calculation using the expression 3 for k(P), 


substitute \ into eq 22 and obtain, aft: 


8 
k[Au-4 (5: q) uJ -A*u=0 in B, u=0 on C 
i 


In view of lemma 1 and the first part of this pro 

the existence of nonnegative eigenvalues \* 

equivalent with the existence of nonpositive eige! 

values of Au+gu. On the other hand, A*=0 mean 
5 8 

Ao h? 

k(Au+gu), if and only if ¢(P,R;r) has a pole 
1=r<0. This completes the proof. 
Corollary: If 


and by lemma 2 such a \ is an eigenvalue 


g(P) p V(P) 


where p is a parameter, lim /,(2, #) exists, if and 


only if p is less than the lowest eigenvalue of Au 
V(P)u. 

Proof: The eigenvalues of Au— V(P)u are obtained 
from those of Au+(p—\(P)u by subtracting » 
Therefore the eigenvalues of the latter expression aré 
positive, if and only if those of Au— Vu exceed p 


IV. Practical Bound for Validity of Random 
Walk Procedure 


If the method described in the preceding section is to 
be used, we have to be sure that the condition of 
theorem 2 is satisfied. Since the eigenvalues are, in 
general, not known, the following sufficient criterion 
may be useful. 

Theorem 3: The expected value E(P. R) exists. if in B 


W(P\< x ( 

a* 
where a and b are the sides of the smallest rectangl 
with sides parallel to the axes that contains B. 

Proof: We substitute for the problem defined by 
formulas 1 and 2 the similar problem obtained 
replacing g(P) everywhere by its maximum y in 
by replacing the domain B by the cireumscrib 
rectangle #2, with sides aand 6, and finally by 
posing the boundary condition: u-=1 on the bounda 
C, of B,. Consider the mass upon arrival anywhe 
on (C, for this new problem. Its expected val 





if it is finite —is at least equal to the expected 

ie(P, R) of the original problem, and it solves 
iodified difference equation problem. 

m theorem 2 we know that /\(/) is finite, if 
yu—0 has only positive eigenvalues correspond- 
o the domain #,, i. e. if the smallest eigenvalue 
Au 0 exceeds ¥ 

e eigenvalues p; of Au=0 in a rectangle can be 
ilated in literal analogy to the familiar procedure 
he Laplace equation. (See [2], p. 258). They 
out to be the values of 


wn rh / 


in ’ 
Ib \ 


ettin | and increasing the expression by 
substituting the first two terms of the power series 
for the sines the proof of our theorem is at hand 


V. Green's Function 


Using the procedure and the notation of the pre- 
vious sections let d, ar. (/) be the expected amount of 
mass that passes through an inner point 4 on the 

-th step of a walk that starts at P and ends on the 
boundary C. Clearly, 

Pd 


‘ 


bP. Q), g,(R, Y=0. (n>0), Ron C 


The difference equation 


k(P ; 
(26) 


, 
qd / .@ 


2 g(P,Q 


G,(P.Q 


GP, @)= lim G,(P, Q), 


G(P, @), if it is finite, solves the problem 


g(P)u (P,Q) in B, u 


Qon C' (29) 


t 
h k(P 


h corresponds to eq 12. When there is no indi- 
m to the contrary the symbol A is always meant 
rate on the point P. 

e solution of 

Oin Bu 


Au+gqu+f(P O0on C 


be written in the form 


u(P) 


h? — 
SO GP, ORO) 
4 «eB 


In analogy with the terminology for differential 
equations we shall call the function 


h 


} 


K(P, Q) -k(Q)G(P,Q) (32) 


Green’s function for our difference equation and the 
domain B. 

The following experiment leads to a random varia- 
ble whose expected value is u(P) provided G(P,Q) is 
finite. 

On every step of the random walk multiply the 
amount of mass at the instant by the value of 
k(Q) #(Q) at that particular point and add the prod- 
ucts thus obtained. If the cumulative sum after N 
random walks starting at P is multiplied by awe 
have an estimate of u(P). 

As to the existence of G(P,Q) we have a theorem 


analogous to theorem 2 


Theorem A: If g(P): in B, then lim G, (P,Q) is 


h? ae 
finite for all Pr and A] in B, if and only if all elge Ne 
values of Au qu in Bare positive : 


Proof: Define ¥(P,Q; r) by 


Sr"ga(P.Q 
n= 


ViP.Q;r) 


The convergence for sufficiently small jr. is shown 
exactly as in the case of the function @(P,Ry 
From formulas 25 and 26, ¥ is seen to solve the 
problem 

6 P.Q)in By—vone 


k(Ay-+ av) + 4( _ )y — (3 
34) 


It is therefore a rational function of 7 The poles of 
} l : 
r and the eigenvalues z3( l ) of k(Ay+-qgy) are 


shown to correspond to each other, just as in the 

proof of lemma 2. The reasoning is somewhat sim- 

pler here, since instead of I'(/ 0 we obtain here 

immediately SOU (P)4(P,Q)=—0. If this is true for 
PeR 

all @ in B, it follows that ((P)=0 and the contra- 

diction is at hand. The subsequent reasoning is 


precisely the same as for theorem 2. 


VI. Random Walks in Unbounded Domains 


If B is an unbounded domain other than the full 
plane, then the quantity £(P?,R) may still exist. The 
function G(P,Q), on the other hand, and its proba- 
bilistic interpretation may exist even in the full plane 
We limit ourselves, therefore, to a discussion of this 
latter quantity. 

Some facts from the theory of infinite matrices will 
have to be used. For convenience we summarize all 
the definitions and properties used, and add _ refer- 


ences to the corresponding pages of [7] 
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a) If A dy}, (, k=1,2,-++, is an infinite 
matrix, the finite matrices A, obtained from A by 
letting ¢ and & run from 1 to j only are called the 
seqment matrices of A. (p. 121). The set consisting 
of all eigenvalues of all segment matrices and their 
accumulation points is called the segment spectrum 


of A, (p. 124). 
(b An 
quadratically convergent, if >> «+, 7 Is finite, (p. 125). 


ce) An infinite matrix called bounded, if 
there is a constant M independent of n such that 


( >* Gut ys) =M-SO 2, ?-S>\y?, (p. 124). The 


_—— 
if l I 


infinite vector s,, (j—1, 2,---) is ealled 


is 


product of two bounded matrices exists and is 
bounded, (p. 131). For the multiplication — of 
bounded matrices the associative law is true, (p. 131). 
If the segment matrices are normal, i. e., if A, A 

A, A,, then A is bounded, if and only if its segment 
spectrum is bounded, (p. 124 


d) An infinite Hermitian matrix is called positive 
definite, if there is a positive constant yu, independ - 


ent of n, such that for all x 
A pF Ey Su p> r,*» (p. 124). 


It Is called nonnegative definite, if 


>) ayr,%>0, for all n and all s,, (p. 124). 
kel 


e) If there exists a constant C such that $3 a, 


k=1 

C, >> | (, (¢=1, 2,---), then the infinite 
matrix =(a,) is bounded (p. 153 

(f) If A is bounded and normal, then it possesses 


a bounded inverse if and only if A A’ is positive 


definite (p. 138). 
(g) If A=(a,) is positive definite, then it possesses 


a bounded inverse .17! a,) and, if A>! 
the inverse of the segment matrix A,, then 


Is 


(a; 


for all i and &, (p 


iy lim a 229). 


h) If A* is positive definite and A is nonnegative 
definite, then A is positive definite 

The definitions at the beginning of section 4 can 
be applied without significant change to an un- 
bounded domain. 

Let the points of the domain be numbered in such 
a way that the first 7 points form, for every j, a 
simply connected domain &, of lattice points. If 
the expected value G(P,Q) is finite for a given pair 
of points P and ¢/, then the expected value G’ (P,9) 
exists for the corresponding random walk problem 
in all those domains 2, which contain P and 4, i. e., 
for all sufficiently large 7. One has, furthermore, 


lim G°(P, Q)=G(P,Q). (35) 
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We shall call the expression Au-+gu posi 
definite, if the real and symmetric infinite matriy 4 
belonging to it is positive definite. In this ma‘. ix 
the rows and columns will always be assumed to je 
arranged in the order of the numbering of the poi. is 
of B which was described above. (But the posit. y, 
definite character of an infinite matrix is, of cou 
not destroyed if the rows and columns are subjec: od 
to the same rearrangement.) 

Theorem 5: If ASu+-qu is positive definite i) 


: - as 
unbounded domain B, and if MP) Sj . in B, then th 
2 


expected value G(P,Q) exists. 


Proof: Since Au 
a fortiori positive definite in all domains P, 
matrix belonging to Au+gu in Rf, is the segment 
matrix A, of A. By theorem 4, the expected value 
G”(P, @) is finite,—provided P, Q are points of 
B,—, and solves problem eq 29. Denote by A, (P, () 
the elements of the matrix A, in the row and column 
corresponding to P and (, respectively, and by 


+gu is positive definite in 2, it is 


| ilt 


AS” (P, () the corresponding element of the invers 
matrix. Then 
’ > t ‘ > v ’ . ’ 
GF () > Oe “(P.S)-k 'US)5CS, Y) 
h? SeR 
} es 
i (P, QOkE-"(Q) 


From this equality we conclude, by virtue of relation 


35, and of the property (g) stated above, that 
(? (P,Q) exists and is given by 
1p 4 ynyp 
G(P, Q) j A™”(P,Qk-(Q) 3 
¥ 
where AW" (P, Q) are the elements of the matrix 


Aa” 


Corollary: If, in addition to the assumption of 
theorem 5, it is known that g(?)>const.> ein / 
then @ (P, 4), as an infinite matrix, is bounded. 

Proof: For then, 0< k7' ((@)< const. - min B 


and the right member of eq 36 defines a bounded 
matrix, since A~' is bounded by virtue of (g). 
Theorem 6: If GP. (}) erists for all P and ), ther 
gu ix nonnegative definite 
Proof: If GOP, @) finite, so are all GY (P, | 
and therefore all segment matrices A, are positiv: 
definite by theorem 2, hence <A nonnegativ: 
definite, by definition (d). 

Theorem 6 cannot be strengthened to a full con- 
verse of theorem 5. This can easily be seen from 
the following example. Let g(P) be zero, so that 
k(P)=1 and eq 1 reduces to Au=0. Then @(P, 
is simply the expected number of visits at ( for a 
moving point starting from P and being absorbed 
as soon as it meets the boundary C for the firs 
time. This quantity is known to be finite for 
infinite quadrant but infinite for the whole pla 
(see [5] ), although Aw is semidefinite in both the 
domains, in consequence of the results of section 
But the following weaker statement can be prove 


Au 


Is 


Is 








eorem 7: Assume that G(P, () erists and is a 
led infinite matrir If then g(P) satis fies in 
hounded domain B the inequality 


q(P)==5—«, (€>0) (3 

g h? . 
Au qu is positive definite in BR. 

Proof: If GUP, () exists then it satisfies the dif- 
we equation 29 in B, i. e., if ACP, @) are the 

ents of A, 


, ’ y ’ 4 . } 
24 A(P, S)G(S, Q) re iP, Q)k-"(P) 
In other words, if T denotes the infinite diagonal 
matrix formed with the elements &#(P?) and © de- 
otes the infinite matrix with elements @ (P, (), then 


Pr is a right inverse of A The matrix A is 
t 


bounded, thanks to assumption 37 and property (é) 


since T is bounded, because of assumption 37, the 
h? . ° . 

matrix ‘ GY is bounded in consequence of proper- 

tv (¢ Thus, A is bounded and normal (even 


symmetric) and possesses a bounded inverse. There- 
lore A’ is positive definite by property (7), and 
property h) assures the positive definiteness of A 

elf 

Remarks: In view of the facts mentioned after the 
proof of theorem 6, it follows from theorem 7 that 
the expected number of visits at a point Y during a 
random walk starting from a point P of an infinite 
quadrant forms an infinite matrix that is unbounded 


VIl. Sampling Method for Calculation of 
Lowest Eigenvalue 


Che fact that the radius of convergence of the 
ctions @ of formula 15 and y of formula 33 de- 
mine the lowest eigenvalue of ki Au qu) can be 
sed for an approximate experimental determination 
this eigenvalue. This gives us an approximation 
the lowest eigenvalue of the differential expression 
) ou 
») oy" 
' nded by a@ plecewise smooth eurve ce 
ng in its interior the same lattice points as B. 
For k approaches 1 as fast as h?, as A —> 0, and the 
nvalues of Au gu approach those of the cor- 
sponding differential equation, (ef. [1] 
he following lemma is needed for the subsequent 


+ygu in the domain B or in any domain 3B’ 


and con- 


ission 

emma 3: The poles of the function ¢(P,R;r) (and 
those of Y(P?,Q:r)) lie symmetrically with respect 

0. This is equivalent to the statement the 

eigenvalues of k(Au+-gu) lie symmetrically with 
} 

h 

wf: The argument applied in the proof of lemma 
he matrix .1 belonging to du+gu can be literally 


ct tor 


extended to the matrix JAD of the proof of theorem 
2, with the result that the eigenvalues of k(Au+- qu) 


: Ss 
are the same as those of I | au 1 (9 he But 


in the proof of theorem 2 it was also shown that, 


“ 8 
if X is an eigenvalue of k(Au+-qu), then jz is an 
¥ 


I: [ a. \ (y 2 )u } Combining 


these two facts the proof of the second statement of 
our lemma is at hand. The first statement is 
equivalent to the second by virtue of lemma 2. 

We first discuss the use of the function @ for our 
purpose. The numerical computations will be 
simpler if, instead of operating with @ directly, the 
function 


eigenvalue of 


o*(P.r) S35 oP, R;r) 


Ret 


is used If we define ¢ nd Fig by 


Ret 
then formula 15 implies 
2 
o*(P)= do r"er(P 38) 


The function @*(/’) satisfies the difference equation 
16 in Band assumes the boundary value 


o*(R 1, for Ron C 


Let A, be the lowest eigenvalue of &(Au-4 qu By 
lemma 2 and lemma 3 the radius of convergence .of 


o(P.R:er) is then 
h? 
a (1 { AN ) 


The radius of convergence of @*(P.r) is, in general, 


also 7 Exceptionally, it may be larger, owing to 


cancellation of poles in the sum >) @(P,R3 In 
Pet 


this case the method to be described would lead 
to some higher, rather than to the lowest eigenvalue. 
If this is suspected, the result may be checked by 
repeating the computations with a different P? or 
with formula 15 directly We shall exclude this 
case from our considerations 

The quantity eF(2/) is the expected mass upon 
arrival anywhere on the boundary for random walks 
starting at P and consisting of exactly n steps. By 
performing a sufficient number of random walks 
starting at P those ¢*(P) for which n is not too large 
can be estimated experimentally 

In order to find 7, from these data we recall that 
@*(P) is a rational function whose poles are real 
and symmetric with respect to the origin. Hence, it 
is of the form 


> > [“ ax rd iP dP 34 
VQ } f’. / ry, 


* P:, 
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Here, N is the number of points in B. For even N 
the quantity d,(P) is zero. ¢,(P) is, by assumption, 
not zero for the point P chosen. From formulas 38 
and 39 we find 


e(P)= >. f; 


2 le.(P) +(— 1)"e_ (P)], 
<N 


If nm is not too small, we can write, with sufficient 
approximation, 


(P)), 


By using several values of n, the estimate of r can 
be improved, and at the same time the admissibility 
of the approximation used is tested. The lowest 
eigenvalue is then given by 


4 i 
- iA! “ty 


The function ¥ of section 5 can also be used to 
compute A,. One defines ¥*(P,r) and g*(P,r) by 


¥°(P.r)= SO WP,.Q;r), 
Veh 
and 
gi (P) So 9A(P.Q:r). 
Ve R 
Then v"(P,r) Dor" gP) g-(P is the expected 


mass anywhere in Bat the nth step of a random walk 
startingat P.  Thisquantity can be found experimen- 
tally by performing a sufficient number of random 
walks starting from P?, and recording the massat every 
step. The sum of all recorded masses for the nth 
step divided by the total number of walks performed 
is an approximate value for g*(P). From thereon 
the procedure is exactly like the one described for 
o*(P;r). 

If the domain B is unbounded, the concept of the 
lowest eigenvalue has to be replaced by that of the 
left endpoint of the spectrum, which may, of course, 
beat —o. We limit the discussion to cases in which 
the spectrum is bounded from below. 


Since no experimental method can take into ae- 
count random walks beyond a certain length, all 
experiments take place in some finite subdomain 
B, of B, so that we must answer the question whether 
the results thus obtained approximate the correct 
values for the given unbounded domain. This 
question is answered affirmatively by the theorem 
that the left endpoint of the segment spectrum 
coincides with the left endpoint of the spectrum, 
provided the matrix is halfbounded and symmetric 
(\7] p. 231). 
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VIII. Comparison with Method of Donsk 
and Kac 


In this section we shall assume that we are deal 
with a one dimensional problem. This will ena 
us to add a few heuristic remarks on the relations 
between our method and that of Donsker and K 
(3], without introducing notational complications 

Let x, be the position of a point at the jth step 
a random walk in a linear net of mesh length 4. T 
particle is supposed to start from s=0 and to 
equally likely to choose either of its two neighbors at 
each step. Then s; is the sum of 7 identically and 
independently distributed random variables with 
mean 0 and standard deviation h. Donsker and Kae 
are primarily concerned with the limiting distribution 
as h->0, ofa random variable which, in our notation, 
can be written (our g(r) is—2V (2) in their notation 


_h 


ath, t) 


5 p> g(x)). 
Here ¢ is an arbitrary positive parameter. 

The quantity g*(P?) of section 7 is, in one dimen- 
sion and for P=0, the expected value of the random 
variable 


y(h,n)= I k(x) 


n <n 
For small A, and with 
t=—h?n 


we have, approximately, 


a@tA.t) 


y(h,n)~e (40 


is the expected 


Our function ¥*(P,r) of section 7 
value of the random variable 


Dor" y(h,n). 


n=() 


With the substitution 


= ‘2 
and using relation 40 it is seen that, approximately 
for small h, 


l . 
So r*y(hn)~ J e~tt- ah Ode. 
n=0 h*, 0 


If ofh.at) is the distribution function of a(A,t), w 
have, therefore, 
* > l ~s em h, t) st 
vy (P,r)~75 ’ ie d,ath,t) (4 
Jo J0 


On the other hand,y*(P,r) satisfies a differen 
equation problem which differs from eq 34 only 





















he factor 6(P,4) in the right member is to be 
ed by 1 and the factors 4 by 2. This last change 
essary because we are dealing now with the one- 
sional problem. Because of eq 42 we may 
for small A, 


hermore, &~1, and the difference equation may 
eplaced by the corresponding differential equa- 
With these changes, we obtain from the differ- 
equation 34 the differential problem 


d*u 2 
g—2s)u4+ 


dx?" ‘s 


Oi B,u=OonC, (34*) 


vhose solution will be approximately y*. In the 
paper by Kac B is the whole line. In this case, the 
solution of eq 34* at s=0 can be written in the form 


u ; x(x) dz 44) 


wre x(r) is Green’s function on the whole line for 
he differential equation 


1 d*u q 
saatls s)u 0 


(45 


For small A, the left members of eq 43 and 44 are 
approximately equal. If we denote by o(a,t) the 
imit of the distribution function o(h, @, t), as h->0, 
t seems therefore plausible that the identity 
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é . “dala,tydt 


y(r)da 46) 


will hold exactly. This identity, proved by Kac in 
[4], is the theoretical basis of the method of Donsker 
and Kae in [3] 

In [4] the passage to the limit, as A->0, 
since that paper is concerned with the distribution 
of certain Wiener integrals. The sampling method 
of [3] operates, of course, with finite sums only. The 
methods of [3] and those of the present paper are 
therefore closely related. In order to be sure that 
the results obtained approximate those for the corre- 
sponding differential equation we refer to [1]. The 
analogous part of Kac’s theory in [4] also makes use 
of the methods of [1] 


is essential, 


Dated: March 21, 1950 


IX. References 


{1} R. Courant, K. Friedrichs, H. Lewy, Ueber die partieller 
Differenzengleichungen der mathematischen Physik 
Math. Ann. 100, 32 to 74 (1928 

[2] R. Courant and D. Hilbert, Methoden der mathematischer 
Physik, 1 (Berlin, 1931 

(3) M. D. Donsker and M. Kac, A sampling method for 
determining the lowest eigenvalue and the principal 


eigenfunction of Schroedinger’s equation, to be pub- 


lished in J tesearch NBS 


[4] M. Kac, On distributions of certain Wiener functionals 
Trans. Amer. Math. Soc. 65, 1 to 13 (1949 
[5] W. H. MeCrea and F. J. Whipple, Random paths in two 


and three dimensions, Proc. Roy. Soc. Edinburgh, 60 


281 to 298 (1930-40 


[6] J. Petrovsky, Ueber das Irrfahrtproblem, Mathematisclhe 


Ann. 109, 425 to 444 (1934 
[7] A. Wintner, Spektraltheorie der unendlichen Matrizen 
Leipzig, 1929 


Los AnGeLes, March 31, 1950 














rnal of Research of National Bureau otf Standards 


Vol. 46, No. 1, January 195] Research Paper : 


Checking and Interpolation of Functions Tabulated 
at Certain Irregular Logarithmic Intervals 


Herbert E. Salzer 


Many functions f(r)behave as polynomials in log + 


When tabulated for arguments 


in geometric progression, f(r) can be checked by ordinary differencing, and interpolation can 


be performed to a fine extent with existing tables of Lagrangian coefficients. 


But in practice, 


f(r) is often known or calculated at some or all of the points 1, 2, 5, 10, 20, 50, 100, 200, 500, 


and 1000 (same theory for the points 0.001, 0.002, 0.005, 0.01, ete., 


or any constant multiple of 1, 2, 5, 


etc., » « 
The present tables have a twofold use: 


laied at some of the more frequently occurring combinations of points 1, 2, 5, ete. 
1 \ I ; 


or 0.01, 0.02, 0.05, 0.1, 


Checking the correctness of f(r) when tabu- 
This also 


includes their use to estimate the least number of tabular entries for interpolation of given 


accuracy 
given by 


W.d 


Tavilor 


I. Introduction 


Functions which behave like polynomials in log x 
are encountered in numerous fields, such as statistics, 
actuarial studies, economics, biometrics, electronics, 
nuclear physics, biophysics, physical chemistry, ete. 
When a function obtained experimentally is suspected 
to have this form, it is usual to examine this by plot- 
ting the values on semilogarithmic graph paper and 
observing whether they lie on a smooth curve. When 
it has been decided that the function is approximately 
a polynomial in log x (or when it is known to be ex- 
actly a polynomial in log +), the question arises as to 
how this information can be used to facilitate check- 
ing the table and interpolating in it numerically. 
This is the question to be discussed here. 


II. Arguments in Geometric Progression 


If the tabular arguments are in geometric progres- 
sion, such as 1,2, 4,8, . or 1, 10, 100, 1000, 
the problems of checking and interpolation are quite 
simple, because the function is effectively tabulated at 
a constant interval in log z. Thus an examination of 
the ordinary differences of the tabular values will 
reveal any errors, as well as the lowest degree which 
an approximating polynomial in log « must have in 
order to vield a certain desired accuracy. Either the 
differences themselves can be used for interpolation, 
or if we wish to avoid using differences, any one of 
several well-known tables of Lagrangian interpola- 
tion coefficients may be used for interpolation to a 
very fine subdivision of the tabular interval. The 
most extensive of these tables of interpolation co- 
efficients are contained in [1].! 


III. Divided Differences for Checking 


In practice, however, it is often found that the 
function has been determined at some or all of the 
points 1, 2, 5, 10, 20, 50, 100, 1000 (or a con- 
stant multiple of those values, e. g., 0.1, 0.2, 0.5, 
1.0, .. .). Then it is no longer possible to use the 


Figures in brackets indicate the literature references at the end of this paper. 


Il. Facilitation of Lagrangian interpolation by 
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a generalization cf a scheme 


ordinary differences for checking purposes Instead 
we use a certain generalization, known as “divided 
differences,”’ a subject that is treated fully in most 
textbooks on finite differences (see [2] or [3]). The 
(n 1)th divided difference with respect to log r. of 
a function f(r) tabulated at the n points a, 2, 


n 
can be written in the form S*A, f(r,), where the 


Ins — 


! 

A, are certain numbers depending only on the 2, 
not on the function. It is convenient to tabulat: 
these quantities, once and for all, for the usual sets 
of values of s, and for the various values of 
so that the divided difference can be obtained by 
a single accumulation on a caleulating machin 
This has been done here for n=3(1)10 and x, ranging 
from 1, 2, 5, ., to 1000, in table 1, A; for n 
and s,; ranging from 1, 5, 10, to 1000 in table 
1, B; and for n=3(1)7 and vs, ranging from 1, 2, 10, 

to 1000 in table 1, C. All entries are given to 


n, 


3(] ‘ 


eight significant figures, correct to about a unit in 
the last place. 
The explicit expression for A, is 
A l/r’ (logyar;—logyr,), | 
where x’ indicates the product over all 747. It ean 
be shown that if f(z) is a polynomial of degree n—2 in 


log x, then the (n—1)th divided difference is identi 
cally zero. If the (n—1)th divided difference 
sufficiently small, then f(r) considered as a polyno- 
mial of degree n—2 in log ys, is in all likelihood, fre¢ 
from error. (However the user is cautioned that a 


Is 


n 
sufficiently small value of $).A;,/(z,) is only a strong 
I 


: 
indication, rather than conclusive evidence, of th 
correctness of the n entries, considered as values of « 
polynomial of degree n—2.) Now the tabular entries 
may be entirely correct and the user might still wan! 
to employ the most economical Lagrangian formula 
Again, this information is conveyed in the sufficient) 
small (n—1)th divided difference, which establishes 
the adequacy of an (n—1) point Lagrangian inte! 

polation formula. In brief, if the function ts check: 


ry 
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1. for checking and interpolation, based upon the points x, 











A,’s for a certain n, we use a different set of 
for interpolation. 


by the 
A,’s corresponding to n—1, 


IV. Interpolation 


Once we have ascertained, by use of table 1, that 
a function is adequately represented by a polynomial 
of degree n—1 in log x» (remembering that n—1 is 
here the »—2 of IIL.), we can interpolate by use of a 
rearrangement of Lagrange’s formula, which was 
suggested by a paper of W..J. Taylor [4]. The original 
method and notation of Taylor are not described in 
this article because he developed them only for the 
special case of equally spaced arguments, and here 


they would be superfluous. The interpolation 
formula is 
f(x) ~( Ya f(r) S\ a (2) 
—_—— — 
t= | l 
where 
a, A log a log 7 (33) 
and A, is given by (1 


V. Logarithms to Other Bases 


The coefficients A, have been computed from 
common (base 10) logarithms, according to (1). It 
would have been possible, instead, to have the A,’s 
calculated for any other logarithmic base. Further- 
more, it is permissible when computing the a,’s from 
(3), to use logarithms to a base different from that 
underlying the A,’s. To illustrate this point, natural 
logarithms (to the base ¢e=2.71828 183) have been 
used in the two examples in section VIII. Extensive 
tables of log.x are given in [5] 


VI. Conversion to Arguments in Geometric 
Progression 


If a large number of interpolations are required 
for a function given for any one set of values s,, it 
may be convenient to prepare, by the method in 
section IV, an auxiliary table giving the function at 
a new set of values Yi, where the y, are how mp geo- 
metric progression Then we can use ordinary in- 
terpolation formulae and coefficients, as suggested 
in section II. This point was called to the writer's 
attention by Churchill Eisenhart and Julius Lieblein, 
of the Statistical Engineering Laboratory. 


VII. Use of A; for Other Arguments 


If the given set of values x, does not begin with 1, 
the table of coefficients A, may sometimes still be 
used by a suitable change of the independent varia- 
ble x. Notice that a polynomial in log (az) or log 
(6/r) is still a polynomial in log rv. Thus, for instance, 
if a function is given at the points r= 10, 20, 50 and 
100, we may consider it as a function of 2’ = 100/z, 
given at the points 2’=1,2,5and 10, so that table 1,A, 
may be used. Again, if a function is given at the 
points r=20, 100, 200 and 1000, it may, instead, be 
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considered as a function of 2’=1000/r given at 
points 2’=1,5, 10 and 50, so that table 1,B, may 
used. Schedules A, B, and C list such transforn 
tions of the independent variable, giving both 
sets of x, for which the change of variable to | 


applicable, and the new arguments z;, which 
proportional to 1/z,. 
Transformation schedules 
ScHEDULE A, 
J I I I I I I 
2 5 2 50 2 500 5 20 5 () 
5 2 5 20 5 200 10 10 10 100 
10 l 10 10 10 100 20 5 20 ‘) 
20 5 20 50 50 2 50 %” 
50 2 ni 20 100 l 100 () 
100 l 100 10 200 
200 5 500 3 
500 2 1000 
1000 l 
ScHEDULE B 
J I r I I I 
2 50 2 500 5 10 5 100 
10 10 10 100 10 5 10 7h) 
20 5 20 0) 50 l 50 10 
100 l 100 10 100 > 
200 5 500 
1000 l 
SCHEDULE C 
J J i I I 
5 20 5 200 2 10 2 low 
10 10 10 100 10 2 10 20 
5O 2 50 20 20 l 20 10 
100 l 100 10 100 Z 
500 2 200 


1000 l 


VIII. Illustrations of Use of Tables 


The following two illustrative examples show how 
to use the tables of A 


Kram ple l. Give n 
r fir ; (ry, 
l 0. 52 50 6333. 56 
2 1. 11 100 15752. 58 
5 50. 30 200 34648. 66 
10 343. 74 500 85417. 78 
20 1447. 15 1000 155833. 72 








degree polynomial in log x (any base Hence only 
seven points are needed for the interpolation. The 
work may be arranged in the following way: (Note 
that an extra place is carried in some columns to 
reduce the error that would arise in rounding before 


leulate {(18 In the absence of knowledge of 
‘plicit formula for f(x), we find that its seventh 
led difference, employing the coefficients A; in 

|, A, for the eight points from 1 to 200, is less 
0.05 in absolute value. This indicates that the 


function behaves, up to the last place, like a sixth 


the final answer is reached. 


25 log. 18 log ry a, A log. 18 log.r, a.ftlz 
| 2. 89037 176 0. 371942 0. 193 
2 2. 19722 458 2. 288576 2. 540 
5 1. 28093 385 11. 903106 598. 726 
10 0. 58778 667 38. 427634 13209. 115 
20 0. 10536 051 144. 713527 209422. 181 
50 1. 02165 125 4. 921950 31173. 466 
100 1. 71479 843 0. 626925 9875. 686 
) > 168. 859664 Sa f(z 200737. 137 
1 fat 
Sa f(x Ja 1188. 78 
The answer found here, 1188.78, is correct to its ; To calculate {(160 To use the present tables of 
st significant figure. The function f(z) was chosen | coefficients A, it is necessary to change the variabli 
é . T to a’ LOOO 7 See Transformation schedule *B 
to be (log. 3(log,2)°+-—- » Fo _ 
; 6 If fla called f(a then 7(160 (6.25 As in 
Krample 2. Given the previous example, in the absence of knowledge of 
rn an explicit formula for f(r), we find that the absolute 
0 ~ value of its third divided difference, employing the 
coefficients A; in table 1, B, for four points, is less 
than 1‘ units in the fifth decimal place. This indi- 
20 50 15. 45081 — - 7 : a. ™ = * as 
100 0 10. 07131 cates that the function behaves, up to the last place, 
200 5 54. 18217 like a second degree polynomial in log z Hence 
1000 95. 09949 only three points are needed for the interpolation, 
Which is carried out as before 
log, 6.25 — log 1, log, 6.25 — log at 
i} 
‘) l 1. 83258 146 0. T8068 92 74. 24314 
) 7) 0. 22314 355 21. 20841 84 1153. 09453 
' 10 0. 47000 363 7. 06787 75 283. 21911 
Sa 27. 58560 67 Sa. f 1362. 97050 
> f | %%e 10. 40876 
This answer, 49.40876, is also correct to the last IX. References 
e given, the function f(r) having been chose 
2 9) (log. 1.43 log.z 1] Tables of Lagrangian interpolation coefficients Columbia 
™ Universitv Press, New York, 1944 
me 2} L. M. Milne-Thomson, Caleulus of finite differences 
ch. I, p. 1 (Maemillan and Co., London, 1933 
" 3] I T. Whittaker and G. Robinso Calculus of observa 
1¢ author expresses his appreciation to Churchill tions, ch. II, p. 20 (Blackie and Son, London, 1944 
I thart and Julius Lieblein for their thorough 4] W. J. Tavlor, A method of Lagrangian curvilinear inter 
y and constructive criticism of this paper in its polation, J. Research NBS 35, 151 (1945) RP1667 
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